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Section 1 


INTRODUCTION 


This report describes the initial attempts by 
the Radiation and Electromagnetics Division of SCIENCE 
APPLICATIONS, INC. to. develop a detailed model of the 
dynamics of arc discharges in electron-irradiated dielectrics 
This work, which was sponsored by the NASA Lewis Research 
Center, will be described in detail in the following 
sections. In this introduction, we note the requirements 
which we feel a detailed model of arc discharges must 
possess. These requirements and observations have motivated 
our attack on the problem — we have chosen a first principle 
approach to provide a basic physics model of the process. 

We summarize this model briefly after reviewing the general 
requirements. 

Any model of "the breakdown and discharge of elec- 
tron irradiated dielectrics should possess most of the 
following features — features which have been gleaned from 
the extensive literature on observations of such phenomena. 

(1) The model should admit permanent damage 
patterns (Lichtenberg figures) in the 
discharged sample. 

(2) It should also provide for the preferential 
channeling of current. 

(3) The propagation of the discharge with speeds 
of the order of 10 4 5 m/sec is also required. 

(4) The model should be capable of explaining 

the optical luminosity of discharges. 



(5) Discharge propagation and damage should 
preferentially occur in regions of large 
trapped charge density . 

(6) The model should contain intrinsically 
statistical features to account for the 
apparent statistical nature of discharges. 

(7) The overall features of the model should 
be generic and not material specific, as 
similar phenomena are observed in many 
different materials. 

(8) A discharge, once initiated, appears to 
continue to propagate independently of the 
field level at which initiation occurred, 
and into regions in which predischarge 
fields are substantially below breakdown 
levels. A discharge model should be capable 
of explaining this fact. 

(9) The model should provide quantitative informa- 

tion about most of the important physical 
quantities: propagation velocity, current 

density, current, temperature, channel size, 
and discharge time. These should be para- 
metrically explainable (in principle) in terms 
of material type, sample dimensions, and 
prebreakdown trapped charge conditions. 

(10) Quantitative predictions should be in agree- 
ment with experiment. 


In this report, we present a model which satisfies 
all the above requirements, with the possible exception of 
(10). Requirement (10) should not be applied too rigorously 
in testing the reported model. The reasons for this are 
simple: while the model is generic, its application to any 

specific material and material configuration requires the 
specification of a large number of materials specific para- 
meters — parameters which to a major degree are relatively 
unknown. Detailed numerical predictions of the model will 
reflect the uncertainty in these parameters. Detailed agree- 
ment cannot be expected when these uncertainties are large. 
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This is not to say that the predictions reported ■ 
herein are in disagreement with existing data — we do not 
believe they are — but rather that the model should be 
judged fibst on its generic features. Thus, the qualitative 
features should be judged, the various scaling features 
addressed, and the order of magnitude of the predictions 
should be consistent. When detailed numerical comparisons 
are to be performed, they should be performed on well 
categorized materials. This type of comparison remains for 
the future. 


With this caveat, the model conforms to the 
requirements. Our approach to the problem utilizes a first 
principles physical description of the underlying electron 
dynamics. The model inputs are descriptions of the basic 
electron scattering processes in the solid. With this informa- 
tion, the model consists of the following: 


(1) A description of the single electron distribu- 
tion function in the conduction band in the 

.presence of high fields. This includes a 
specification of drift velocity, diffusion, 
and avalanche length, as well as mean energy 
— all as a function of electric field. 

(2) A description of the primary spatially and 
.temporally coherent process of electron 
multiplication by avalanche. 

(3) A description of the self-consistent evolution 
of the electron avalanche into a primary nega- 
tive tip streamer which propagates with the 
field, leaving a high conductivity charge 
neutral region behind. 

(4) A description of the propagation of the posi- 
tive tip of the neutral region by successive 
electron avalanche/streamers, this phase oc- 
curring in three dimensions. 
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(5) A description of the current flow, channel 
evolution, charge release, and temperature 
rise associated with the above process. 


The above model, whose physical implementation is 
described below, provides for a complete description of the 
discharge process from single electron dynamics to Lichtenberg 
figure formation. It provides for a quantitative description 
of all relevant physical quantities. Code implementation of 
the descriptions of (1-3) are described in this report. Their 
names are, respectively, SEMC, CASCAD, and ACORN. A descrip- 
tion of the processes involved in (4-5) is also given, and 
a possible code implementation is described. A detailed 
application of the model is presented for a rough material 
model of Teflon. 
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Section 2 


SINGLE ELECTRON’ DYNAMICS 


In this section we give a discussion of our ap- 
proach to the motion of individual conduction band elec- 
trons of a dielectric in the presence of an impressed 
external electric field E. The discussion is organized 
by presenting the details of the theoretical structure 
which has been assumed, followed by a brief discussion 
of the limitations of the approach. This ordering of 
the discussion allows the limitations to be viewed 
against the concrete framework which is presented. Once 
the theoretical basis has been established, we present 
our computational approach to solving the problem which 
has been posed. A final section discusses computational 
results for two particular materials gaseous molecular 
Nitrogen and FEP Teflon. The material models employed 
in these calculations are discussed in Appendices. 
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2.1 


THEORY 


The fundamental assumption of our approach to 
breakdown dynamics is that the dynamics of the electrons 
of interest "may be followed with a linear Boltzman equa- 
tion. 


We assume that the behavior of a single elec- 
tron may -be described with a distribution function 
f(x,v,t), where x denotes the position of the particle, 
v its velocity and t is time. The distribution function 
f is interpreted as the probability of finding the elec- 
tron at position x, with velocity v, at. time t. The 
distribution is assumed to evolve in phase space (x,v) 
according to the Boltzman equation 


3f 

St 


-b 




\ ot / scatt 


( 2 . 1 ) 


The left hand side represents the evolution in 
phase space of the distribution function according to the 
ordinary phase space trajectories of an ensemble of non- 
interacting particles under the influence of the external 
force F = -eE. (The Louiville theorem.) The right hand 
side represents the changes in these trajectories due to 
various scattering processes. The linearity assumption 
requires that the interaction of various conduction band 
electrons may be neglected so that the scattering is 
specified by the interaction of the electron of interest 
with the remainder of the solid. Specifically, the as- 
sumption requires that the rate of change of f due to 
collisions be given by a linear collision integral. 


2^2 



( 2 . 2 ) 


(f)scatt * J d9 o *<*■*„•* > 

where K(v,v q ) is the scattering kernel. The scattering 
kernel K(v,v Q ) consists of two parts, that due to gain 

, and that due to loss K . Thus K. ( v , v ) is equal to the 

• """TO 

rate at which an electron of velocity v Q wil’l scatter 
into the new velocity v, while K_(v,v q ) is the rate at 
which electrons of velocity v are scattered into velocity 
v . The scattering rate is normally thought of in terms 
of the sum of several rates due to individual elementary 
scattering processes. We denote the differential rate 
for the i'th scattering process as 

K±(v,v o ) = ^4 (V— v) . (2.30 

dv 

where indicates partial differentiation with respect 
to all elements of the volume element dv. The notation 
implies that the rate in general depends not only on the 
initial state v Q but also the final state v. The total 
rate for the i^* 1 process v^(v Q ) is obtained by integrat- 
ing over final scattering states 

v jL (v Q ) = fdv dv l (v — -v) (2.4) 

J dv 

The rate for the i'th process is normally related to a 
mean free time t^(v ) for the i'th- process by the equation 
The total rate for all scattering v is the sum 
of all the individual rates. With these definitions, we 
have 

K + (v,v ) =X) (v — v) (2.5) 

i dv 

K_(v,v q ) = v(tf Q ) S 3 (v 0 - v) (2.6) 
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(2.7) 


and K(v,v q ) = K + (v,v o ) - K_(v,v q ) 


where 


6 3 


is the usual Dirac delta function. 


The interpretation of equation (2.1) together 
with eqs. (2.2 - 2.7) .requires some explanation. We are 
interested in processes in which the electron of interest 
(described by f) induces a transition from a bound state 
(trap or valence band) into the conduction band. If we 
were to attempt to describe this multiparticle situation 
with the Boltzman equation (multipart icle distribution 
function described by a mean single particle distribution 
function), then equation (2.1) would require the addition 
of a source term which takes into account the changing par- 
ticle number. Likewise, sink terms would need to be added 
for processes taking an electron out of the conduction band 
(trapping, recombination) . .We expressly ignore these 
processes in Eq. (2.1), and assume that we can meaningfully 
follow- the evolution of a singl-e distinguishable electron. 
Thus, the effect of the ionization process on the motion 
of this single electron is taken into account, but the 
multiparticle-, multi-band nature of the actual configura- 
tion is ignored,. The limitations of this assumption will 
be explored after further development of the theoretical 
framework. 


Our discussion indicates that we are consider- 
ing the Boltzman equation under the circumstances in 
which particle number is conserved. This follows from 
our description of the scattering kernel given by equation 
(2.7). With the usual interpretation, we may define t'he 
total number of particles N described by f to be 


J' 


N ( t ) ~ dxdv f ( x , v , t ) 


( 2 . 8 ) 
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Using the Boltzman equation (2.1) together with 
the scattering kernel property 


I- 


dv K(v,v o ) = 0 


(2.9) 


which follows from eqs . (2.5 - 2.7) 
dN 

that -xr = 0 . Thus, N is a constant 
dt 

veniently be chosen to be unity. 


it is easily verified 
in time which may con- 


Our use of the Boltzman equation to describe 
the motion of the electron is demanded not because there 
are many electrons to describe, but rather because we do 
not know the precise state of the solid during the electron 
interaction. Thus, the ensemble average implied in the 
use of the Boltzman equation is the average over the en- 
semble of solid configurations. The usual statistical 
mechanical interpretation permits the use of the Boltzman 
equation for the description of the single electron. 


We have assumed in eq. (2.2) that the scatter- 
ing kernel does not depend on position — the assumption 
of a homogeneous material. For the time and distance 
scales of interest, this is an excellent assumption. For 
kernels of this variety, the usual velocity space form of 
the Boltzman equation may be derived. Defining the veloc- 
ity space probability distribution g(v,t) by 

g(v,t) = fdx f ( x , v , t ) (2.10) 


we find 



Discussions of transport in solids (and gases) normally 
involve eq. (2.11) rather than eq. (2.1). In fact, most 
attacks on transport properties specialize even further 
by considering the equilibrium solution (v) (time inde- 
pendent) of eq. (2.11). 


The relation of these various functions and 
approaches is best seen by considering various moments 
of the distribution function. In the usual fashion, 
define the mean value of a function of position and 
velocity p (denoted <p>) by 

<p> = fdxdv p(x,v) f(x,v,t) (2.12) 


The following relations may then be easily established 
from eq. (2 . 1) ; 


d<x> 

dt 


<v> 


(2.13) 


-r-<X.X . > 
dt i j 


<x . v . + v . x . > 
1 J 1 3 


d<v> _ -eE £c 

dt m m 


(2.14) 

(2.15) 


~r<v.v .> = 
dt i j 


— (E . <v . > + E . <v . >) 
m i j J i 




whe r e 




dv v K( v 




g(v Q> t) 


and 


T. . 
11 


I JdvJdv o v iVj K(v,v o ) g(v o „t) 


(2.16) 


(2.17) 


(2.18) 
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The important point to note is that the evolu- 
tion of the various velocity moments {<v.>, <v.v.>} 

i x 3 

depends only' on the evolution of the velocity space dis- 
tribution function g. In particular, if equilibrium 

has been obtained = 0), the mean velocity and energy 

dt 

are constant, and the mean position increases linearly 
with <v>. Thus, the drift velocity v^ is easily identified 
from the equilibrium solution 


v_ = <v> 

D o 


= j*dv v g Q (v) 


(2.19) 


Define the displacement £ from the mean position 

<x> by 5 = x-<x>. What is not quite so apparent from eqs . 

(2.12 - 2.19) is that <£; . £ . > also evolves linearly in time 

x 1 

when g has reached equilibrium. The proof of this fact 
requires a somewhat delicate application of the Chapman- 
Enskog theory and is relegated to Appendix 9. Suffice it 
to say here that the coefficient of time in the evolution 
of is interpreted as a diffusion coefficient. Let 

e„ be a unit vector in the - direct ion of the electric field, 
and define the parallel and perpendicular components of x 

A A 

by x„ =(e 11 *x)e„, x x = x - x„. For diffusive motion, we 
have the definition of the parallel and perpendicular dif- 
fusion coefficients by the relations 


= 2D„t (2.20) 

<L *L> = 4D x t (2.21) 


These terms may be isolated from the solution f of the 
Boltzman equation using equation 2.12 - 2.19. Formulae 
for the D's in terms of the equilibrium solution. g Q may 
be found in Appendix 9 . 
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Another averaged quantity of primary concern 
in our discussion is the mean distance between ionizing 
collisions', variously called the avalanche length or 
the first Townsend coefficient a. Let v^ on (v Q ) 
ionization rate for electrons of speed v . Then the 
rate n(x,t) at which ionizing collisions are occurring 
at a particular point in space x at time t is given 
by (#/unit volume/unit time) 

n(x,t) = J d ^ 0 v i0n (^ 0 ) f (x,v o ,t) (2.22) 

The rate at which ionizing events are occurring at all 
locations N(t) is obtained b.y integrating n(x,t) over 
all space, 

N(t) =J*dx n(x,t) - j dv Q v ion (v Q ) g(v Q ,t) (2.23) 

The total number of ionizing collisions which have oc- 
curred to up time t, N(t) is obtained by integration 
N in time. When g has reached equilibrium, N increases 
linearly with time. Let tj be the time period required 
for N to change by unity. The inverse of t^ is of 
course the ionization rate 0. The mean distance between 
ionizations is related to t^ by a^v^t-j.. 

The spatial probability distribution n(x,t) 

of ionizing events is obtained by integrating n(x,t) 

in time. A convenient way to disentangle the many 

ionizing events is restrict attention to the first 

event. Let t be a time such that N(t-i) = 1. The spa- 
1 ^ ± 

tial distribution n(x,t=t ) =-n (x) may then be inter- 
preted as the probability that the first ionization 
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occurs at x. The mean location of the first ionization 
is then 

< ^I > = J n 1 (x) x (2.24) • 

The temporal evolution of the ionization may 
then be thought of in terms of a sequence of- events 
identical to the first- event, with spatial probability 
given by n^(x), and temporal probability of occurrence 
N(t) [cumulative]. In so far as t^tj, the definition 
of <x ]; > may be identified with a- This has been the 
standard procedure in the literature [Baraff , 1964]. 

Experimentally, the mean quantities a, , 

D„ and D a are the quantities which have been measured. 

We will refer to specifics of such measurements below 
in discussing our computational results. Quantities 
such as these are also- the natural quantities to use 
in proceeding to a macroscopic description of multi- 
electron dynamics, which will be discussed in Section 
4 below. Note that a specification of a material 
model requires a model for the various scattering 
processes which occur in the description of the scat- 
tering kernel K(v,v q ). 
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2 . 1.1 


LIMITATIONS 


'The description of single electron dynamics 
implied by the use of the Bolnzman equation is purely 
classical. That is, the probability distributions 
which are obtained are related 'to the classical aver- 
aging over ensembles which is customary in statistical 
mechanics. These probability distributions are in no 
way directly related to the probability distribution 
associated with the quantum mechanical description of 
the electron. The fundamental limitation of the method 
is therefore related to the fundamental domain of ap- 
plicability of classical statistical techniques. For 
a general problem involving electronic transport 
properties in solids, density matrix techniques of 
quantum statistical mechanics are required, leading 
the Kubo theory of transport phenomena (Kubo, 1959). 

We believe, however, that the classical de- 
scription which we have invoked for the calculation of 
the electron transport properties is adequate for the 
problem of interest. This technique is not unique to 
our approach and has been invoked by essentially all 
authors who have studied this problem at reasonably 
high temperatures O room temperature). A somewhat 
detailed justification may be found in Ziman's book 
(Ziman, 1960) and references cited therein, while a re- 
view of the literature may be found in 0 'Dwyer's book 
(O' Dwyer, 1973). Though the details of the argument 
can become complex, the underlying idea is quite simple. 
The correspondence principle demands that wave packet 
solutions of the Schroedinger equation obey the classical 
equations of motion. The fundamental limitation of the 
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classical methods comes only when attempting to provide 
more localization in phase space than is allowed by the 
Uncertainty Principle. If, for the energies- and times 
of interest, the classical averaging provides an uncer- 
tainty-greater than that demanded by the Uncertainty 
Principle, then the quantum uncertainty is masked by 
these effects, and classical techniques may be expected 
to be applicable. For example, if we wish to discuss 
(free) electrons of roughly thermal energy (a, 0.025 e.v.), 
then we are justified in applying classical techniques . 
if we do not try to localize the electron any more than 
Axo^-, or about 25 A. For higher energy electrons, this 
length becomes smaller. 

The problem may roughly be factored, with 
classical techniques being used for distances large 
compared to Ax'v|— , and quantum techniques being used 
over smaller distance scales (with corresponding state- 
ments about time scales). If the fundamental scattering 
processes are specified using quantum techniques (de- 
scription of the scattering kernel), the longer distance 
scale transport properties should then be describable by 
classical techniques. That this approach is adequate 
for the problem of concern is taken as a working hypoth- 
esis, a hypothesis consistent with those of other workers 
in the field. 

A second limitation of the approach is more 
severe. In applying a Boltzman type transport equation 
to the problem, we require that the interactions which 
are occurring in the quantum regime are describable in 
terms of a scattering formalism. This roughly requires 
that the electronic states of the conduction band be 
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describable with, quasi-free wave functions (i.e., that 
a quadratic relation exist • between electron energy and 
Bloch, wave ' number) . While this is probably a reasonable 
assumption for a perfect crystalline insulator (trap- 
free), it is decidedly not true for non-crystalline- 
polymeric trap-dominated insulators. Shown in Figure 

2.1 is a schematic representation of the electronic 
energy levels expected in solid insulators with a large 
degree of non-crystallinity (Mott and Davis, 1971). We 
can see that the crystalline concept of a forbidden gap 
has been replaced by a minimum in the density of states. 
Near the top and bottom of what would normally be thought 
of as the valence and conduction bands respectively are a 
group of states which are normally called "trap" states. 
These states correspond .to localized wavefunctions , and 
hence do not have a translational wave number associated 
with them. Localized states such as these (so-called 
Anderson states (Anderson, 1958) cannot be represented 
as quasi-free electrons, so that a scattering formalism 
is inapplicable for a description of their properties. 

In particular, the states just below the "conduction 
band," the shallow trapping states, are the lowest 
lying levels above the "valence band." Statistical ar- 
guments at ordinary temperatures suggest that most elec- 
trons which are not in the valence band in fact occupy 
these states. Because of this, the normal conduction 
properties of a material with an energy level structure 
of the type shown in Figure 2.1 are dominated by the 
behavior of electrons in these states. The normal de- 
scription of the conduction process is then by "hopping" 
conduction; that is, successive detrapping and retrap- 
ping. Techniques other than Boltzman equation techniques 
are required to describe this type of transport . The 
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FIGURE 2.1. SCHEMATIC ENERGY LEVEL DIAGRAM 
FOR NON-CRYSTALLINE INSULATOR 
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background on the required techniques may be found in 
the book by Mott and Davis (1971). 

Despite the above remarks, we persist in the 
use of the Boltzman equation for our model of avalanche 
induced breakdown. Our reasons are as follows. We are 
only interested in performing electron transport calcu- 
lations at or near breakdown fields. It is our working 
hypothesis that in this regime an electron which is 
initially in the conduction "band" (above the trap 
states of Figure 2.1) will undergo significant "heat- 
ing." That is, the rate of energy loss to the solid 
will not equal the gain from the field until the elec- 
tron has an energy substantially greater than typical 
thermal energies-. The presence of the field, therefore, 
induces an effective transition rate upwards in the con- 
duction "band." This effective shift of the energy of 
conducting electrons effectively removes electrons from 
those states which have the greatest probability of 
trapping (within kT of the trap states). Thus, if the 
heating rate is substantially faster than the trapping 
rate, we are well-justified in concentrating attention 
on the conduction "band" dynamics, and treating the 
trapping as a perturbation on these processes. 

The validity of the above assumption is 
clearly a function of field strength. It is our work- 
ing hypothesis that for the field strengths associated 
with breakdown this condition is satisfied. Once a 
calculation of the electron, distribution function has 
been performed, an a posteriori check of the assumption 
should serve to check its adequacy. With this assump- 
tion, it is reasonable to expect that a quasi-free 
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electron description of the conduction "band" electrons 
should be adequate, that a scattering formalism is ap- 
propriate , 'and that the Boltzman equation is applicable 
In fact, the results below suggest mean electron ener- 
gies in the eV range. Because we are well up into the 
"band" in this regime, we expect potential energy ef- 
fects to be less important than kinetic energy effects. 

We uniformly assume therefore that M* * M in relating 

© © 

energy and momentum. 

We have suggested above that we are well 
justified in treating trapping as a perturbation of 
the energy gain rate is substantially greater than the 
trapping rate. There are a whole host of other limita- 
tions of our approach which are based on a similar com- 
parison of rates (of energy gain and conduction band 
scattering). In particular, the deep trapping and re- 
combination rates must also be negligible compared to 
the energy gain rate. The electron-electron energy- 
momentum exchange rate must be negligible compared to 
the electron-lattice scattering rate. Finally, the 
electric field must be sufficiently small that the de- 
tails of the ionization process (energy/angle) do' not 
determine the evolution of the distribution function 
(scattering rate much greater than ionization rate). 

We believe all the above restrictions to be satisfied 
for the fields of interest. 
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2.2 


COMPUTATIONAL APPROACH 


A great deal of detailed work has been per- ■ 
formed on methods for computing, the equilibrium solution 
of the Boltzman equation for gases. In fact, the low 
energy cross-sections for gases have been determined 
from measurements of macroscopic parameters as a func- 
tion of electric field by inverting the process dis- 
cussed above, i.e., by adjusting cross-sections in the 
calculations until an adequate fit of the measured 
macroscopic data, is obtained. We refer the reader par- 
ticularly to the work of Phelps [see e.g. Hake and 
Phelps (1967) and Frost and Phelps (1962) and references 
cited therein] . Further references to the literature 
may be found in the compendium of data prepared by 
Dutton (1975) and, the recent book of Meek and Craggs 
(1978), while the older literature may be found, from 
Loeb's book [Loeb (1955)]-. The techniques developed by 
these workers are quite adequate for the proposed prob- 
lem. Unfortunately, the 3 ^ are rather specialized, in 
that various electric field regions are distinguished, 
various approximations employed, and specialized solu- 
tion techniques are employed in the several regions. 
Further, the application of any of these solutions to 
the spatial and temporal variation of the distribution 
function requires the application of the macroscopic 
diffusion theory. 

In solid dielectrics, much less work has been 
performed on developing numerical transport techniques. 
The state of the art appears to be represented by the 
work of Baraff (1962), who determined a universal set 
of curves for the avalanche length under a constant mean 
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free path assumption. Earlier work, and more complete 
descriptions of the various competing processes may be 
found in Heller (1951) ,. Wolff (1954) and Shockley (1951), 
while a summary may be found in O'Dwyer's book (1974). 
Again, the , techniques employed have either been very 
specialized, or very difficult to employ for general 
circumstances . 

Rather than rely on any of the specialized 
techniques which have been developed, we have chosen 
to attack the computational problem by a brute force 
Monte Carlo technique. This method (to be described 
below) 1 has several strong advantages: There are es- 

sentially no restrictions on the type, number or 
character of scattering processes which can. be con- 
sidered; as the technique involves path tracing, it 
can readily solve problems associated with both the real 
space and velocity space- positions of the distribution 
function simultaneously; the technique is conceptually 
simple to implement — resulting codes are quite small; 
the technique may be used in principle to determine 
any information desired about the evolution of the dis- 
tribution function; finally, SAI is well-versed in im- 
plementation of such codes. The technique has several 
disadvantages; Monte Carlo codes are inefficient — they 
are expensive to execute compared to comparable codes 
using other techniques; Monte Carlo codes provide essen- 
tially no basis for intuitive understanding of the re- 
sults computed; without special sampling techniques, 
quantities depending on low probability events are im- 
possible to compute; the statistical nature of the 
Monte Carlo output often requires that special fitting 
routines be developed for data smoothing. Our decision 
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to go with Monte Carlo techniques was driven by ease of 
implementation and generality. 

The basic concept of Monte Carlo transport is 
extremely simple. The left hand side of the Boltzman 
equation describes the evolution of the distribution 
function in the absence of collisions. It is well 
known that this evolution is described by individual 
particles following their mechanical trajectories, i.e., 
f is constant on phase space trajectories. The scatter- 
ing integral represents the scattering of individual 
particles from one velocity to another. This occurs at 
a rate of v(v o ). Thus, on the average, an individual 
particle will undergo free notion for a time 
t(v q ) =ly^ e (v 0 ). This mean free time may be associated 
with a mean free path A(v q ) by the standard relation 
X(v ) = |v I t(v ). Thus, the picture is that of indi- 
vidual particles undergoing free motion, together with 
successive scatters at a mean time increment of x(v ). 
The scattering events of the scattering integral may 
be described according to single particle scattering 
theory — which is precisely in accord with the method 
by which the scattering integral was developed. The 
scattering events themselves are random events chosen 
according to the differential rate (v q — - v) , again 

in accord with the meaning of the scattering integral. 
Averaging the results over many particle paths chosen 
as above yields the distribution function. That the 
above procedure indeed yields a solution to the 
Boltzman equation is intuitively apparent (the Boltzman 
equation was developed from this picture). The rigorous 
mathematical demonstration of the equivalence -is much 
more difficult. A discussion of the proof of this 
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equivalence would take us too far from our goals — we 
accept the equivalence and refer the reader to the lit- 
erature fo:d further discussion. 


We. have chosen to implement the above picture 
for the transport of electrons in the presence of an 
electric field. The method chosen is a single-scatter 
Monte Carlo method (to distinguish it from other Monte 
Carlo algorithms), i.e., each individual scattering 
event is modeled. The method as implemented is brute 
force — no special sampling techniques have been em- 
ployed. The restrictions resulting from this are dis- 
cussed after a description of the algorithm. The algo- 
rithm is as follows: 


dv i 

0 A set of differential rates — rr=? (v Q — *-v) 
is prescribed, i=l,. N, where N is 
the number of scattering processes to 
be included. 

• The total scattering rate v(v ) is de- 
termined, together with the mean free 
path X(v Q ) = !v o l/v(v o ). 

0 The differential scattering probability 

functions are determined according to 
the rule dP. dv (v — *-v) 

dv (v o“^ v) - dv 

. ■ . 

0 Cumulative probability functions C , v 
are determined by appropriately deiined 
integrals of dP^ 

“dv 

0 An initial velocity space distribution 

function f Q (v) is specified. The cumu- 
lative distribution function, C(v) is de 
termined from f (v) by appropriate inte- 
gration. 
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Particles are started from a preassigned 
position in space x . 

A single particle path is constructed by 
the following sequence. 


(1) An initial velocity v is chosen 
randomly according to°the distri- 
bution function C(v ). 

o J 

(2) The mean free path X(v ) is noted. 

(3) A free path Z is chosen according 
to an exponential distribution of 
paths with mean X. 


(4) 


The particle is allowed to traverse 
a path length l under the influence 
of the external field, with initial 
conditions x Q , v Q . This requires 
the determination of the time <5t 
required to travel a gath length Z 
(this depends on v , E and v Q *E). 
The particle position and velocity 
are updated according to the equa- 
tions of motion 


x = x Q +-v o 6t - l/2(^|)(6t) 2 

v = v - — St 
o m 


(5) A scattering event is described at 
the new location For velocity v', 
a new velocity v^ is randomly se- 
lected according to the probability 
distributions P^(v,Vq). The particle 
now has position x, and velocity v^ . 

(6) The steps (2) - (5) are repeated for 
the new position and velocity. This 
process is repeated until prespeci- 
fied conditions (e.g., fixed number 
of steps) are satisfied - . 


(7) The above sequence of path increments 
represents one Monte Carlo History. 

The sequence is repeated for a prespec 
ified number of histories'. 
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(8) Various quantities are scored during 
the course of each history. These 
quantities are averaged over the 
number of histories. 


The first four steps of the above algorithm have 
been programmed into a FORTRAN code called XSCPRP . The re- 
maining steps have been programmed into a FORTRAN code called 
SEMC (Single-scatter Electron Monte Carlo). The operational 
features of XSCPRP are described in Appendix 3. The opera- 
tional features of SEMC are described in Appendix 4. Some 
general features of the codes are as follows: 


• SEMC is designed to operate in the 
high electric field region where col- 
lision ionization is important. Many 
of the scoring and path stopping flags 
are keyed to these events. For low 
electric fields,, such events become very 
improbable. The code is not meant to 
operate in this regime without special 
direcxives . 

• The initial distribution of electrons 
for the computation is assumed to lie 
in the low energy domain (thermal ener- 
gies). Under high fields, the exact 
form of this initial distribution is 
immaterial to the ultimate development 

of the distribution function. It was con- 
sequently specified to be monoenergetic 
and unidirectional. 

• SEMC only follows each individual 
electron trajectory up until its first 
ionization event. It therefore assumes 
that the electron loses essentially all 
of its energy in this encounter, so 
that its further history is undistin- 
guished from starting with the initial 
distribution. This requires that the 
electric field not be so large that a 
substantial fraction of the ionization 
energy be gained in a mean free path. 
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0 


Sampling on free paths, rather than free 
times was chosen because the mean free 
path is more nearly constant in energy 
than mean free time. SEMC samples from 
the mean free path distribution of the 
electron prior to the time step. In 
principle, the distribution should be 
integrated along the trajectory, and an 
energy weighted path distribution sampled 
from. If the fields are not too large, 
this approximation should be adequate. 


The path length traversed £, is related 
to the time required to traverse that 
distance by the formula 


£" = 1/2 jcos 0 + (t "-cos 9) " 2 -2t "cos 0+1 

o 

+ sin 0 In 
where 


t"-cos 0 +vt 


vt. 


2t "cos 0+1 


(1-cos 0) 


} ( 2 ' 


24) 


£" 
t " 


e 

Eli 

mv^ 

o 

e 

E It 

m 

v | 
o 1 


(2.25) 

(2.26) 


and, 
cos 6 



(2.27) 


For a choice of path length i, eq. -(2.24) 
s inverted, and the new velocity and position are up- 
ated with the equations of motion for elapsed time t. 


# The location of the i'th. ionisation event 
x ± is scored together with the time of 
occurrence t^. These results are binned 
to provide the spatial distribution of 
ionization events. 
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• At the end. of each history, the following 
particle attributes are scored. The mean 
value of the component of x along the 
field provides the ionization length. 

The mean value of- (x-j_/ti) provides t'he 
drift velocity. The RMS speed of the- 
final posit ions - perpendicular to the 
field provides the diffusion co.efficient 
perpendicular to the field D A . Because 
of the strong correlation, between ioniza- 
tion time and distance traversed, the RMS 
speed of the ionization location about the 
mean ionization location does not provide 
a good measure -of the longitudinal diffu- 
sion coefficient D„. Instead -a computa- 
tion of mean dispersion is scored for a 
fixed time (to.be specified on input). 

This provides the longitudinal diffusion 
coefficient. 

* ■ The scattering distributions required 

for performing the random walk are re- 
quired input parameters for the code. 

o The angular dependence of the inelastic 
processes is currently being neglected- 
This is a reasonable • approximation so 
long as acoustic phonon scattering 
(elastic scattering.) dominates. 
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2.3 


RESULTS 


In this section we present the results of the 
single electron dynamics calculation for two materials, 
and CF 0 . One of the features which makes the Monte 
Carlo appealing i-s ease with which a running computer 
code may be developed. Unlike finite difference methods, 
for example, the Monte Carlo method does not have stabil- 
ity criteria and in general convergence is assured. For 
this reason it is imperative that. the Monte Carlo codes 
be checked, if at all possible, against experimental 
data. As our test case we have chosen N2, primarily 
because of the abundance of experimental data available. 
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2.3.1 N 2 RESULTS 

The model for the cross section set used for 
N 2 is presented in Appendix 1. For all results pre- 
sented for N 2 the solid curves are from an electron 
swarm data compilation (Dutton, 1975) and the points 
are data calculated by the SEMC code. All N 2 results 
are for a gas density of 2.7x10 molecule/m 3 . In 
Figure 2.2 the ionization length is plotted as a func- 
tion of electric field. The ionization length is the 
reciprocal of a, the first Townsend coefficient. In 
Figure 2 ..3 the electron mobility (ratio of drift velocity 
to electric field) is plotted as a function of electric 
field. In Figure 2.4, the longitudinal diffusion coef- 
ficients (the coefficient in a direction along the 
applied electric field) is plotted as a function of 
electric field. And finally, in Figure 2.5 the trans- 
verse diffusion coefficient (the coefficient in a direc- 
tion normal to the applied electric field) is plotted as 
a function of electric field. 

On two of the curves, data was not available 
for high values of the electric field. The curves have 
been extended (dashed position of the curve) for com- 
parison purposes. The generally excellent agreement of 
the code calculated results with experimental data vali- 
date the code. 
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2. 3. .2 


TEFLON 


SEMC has been exercised for a material model of 
Teflon which, is described in Appendix 2. The model is 
necessarily crude. Shown in figure 2.6. is mean distance 
for e-folding the number of electrons (inverse first 
Townsend coefficient). This quantity equals the mean 
ionization length divided by In (2.). For an applied field 

Q 

of 8 x 10 V/m, this corresponds to a distance of 
-8 o 

1.5 x 10 . -m (130 A). As. for Nitrogen, the quantity in- 
creases exponentially as the field is decreased. The 
extremely small distance scale associated with these 
high fields should be noted, as. it will determine the 
scale of later computations . These results are consis- 
tent with an analytic estimate of the, avalanche length 
presented in Appendix 7. Shown in figure 2.7 is the ' 

electronic mobility p. Its remarkably flat value of 

-4 2 . 

about 1.. 6 x 10 m /V-sec can be ascribed to the small 

constant acoustic phonon collision mean free path, and ■ 
is consistent with simple estimates. Shown in figures 
2.8 and 2.9 are the longitudinal and transverse diffu- 
sion coefficients. While the order of magnitude of 
these coefficient .is consistent with analytic estimates, 
the rather strong increase in Dip at low fields seems a 
bit peculiar. This behavior remains to be explained. 

This feature does coincide with the onset of signifi- 
cant ionization. The equations of Section 4 require 
the ionization rate 8 = aV D . This quantity is plotted 
in figure 2.10. Its variation with electric field is 
extremely sharp, and reflects the properties of p and 
a noted above. Finally, figure 2.11 shows the marginal 
distributions f.or the locations of ionization -events 
parallel and perpendicular to the field. It should be 
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noted that these distributions show a somewhat greater 
width than would be expected from- diffusion alone. 

This is a reflection of the correlation between ioniza- 
tion location and ionization time, and was discussed in 
Section 2.2 above. 
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Section 3 


ELECTRON AVALANCHE - 


In the previous section we developed, the compu- 
tational apparatus for calculating the production of a 
second conduction band electron due to impact ionization 
by an existing conduction band electron. The newly cre- 
ated electron is expected to undergo uhe same acceleration 
process due to the electric field, i.e., its dynamics 
should be identical to those of the first electron. An 
electron multiplication (avalanche) with drift is ex- 
pected (the Townsend mechanism). For an avalanche gen- 
erated by a single initial electron, we envision the 
evolution of the avalanche to be a spatially and tem- 
porally coherent process (within the confines of the 
Uncertainty Principle discussed in Section 2.1). In 
this section we discuss this phenomenon, our approach to 
computing its evolution, and some numerical results for 
Teflon. 
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3,.1 


THEORY 


Extensive cloud chamber and electrical obser- 
vations in gases have demonstrated that spatially and 
temporally coherent electron avalanches can be generated 
by a single electron. We refer the reader to the excel- 
lent book by Raether (1960) which describes the founda- 
tions of this research, and to the recent book by Meek 
and Craggs (1978) for a current survey of the literature. 
Shown in Figure 3.1 is a schematic representation of such 
an avalanche. A single electron was freed at the tail of 
the figure. As the el-ectron drifted with the field, it 
created a second electron — this pair of electrons in 
turn generated four electrons, and so on. As the elec- 
trons drift with the field, they also diffuse due to 
collisions with, the background gas, while the relatively 
immobile ions remain more or less stationary. 

This picture of single-electron induced ava- 
lanche is conceptually precise. Neglecting complications 
due to other rate processes (attachment, recombination,, 
etc.) it admits, on the average, a simple description. 

The total number of electrons N in the avalanche is 
given by 

N = e ax (3.1) 

where a is the inverse ionization length- (first Townsend 
coefficient) computed in Section 2 above, and x is the 
drift distance along the avalanche direction, given by 

x = v D t . (3.2) 
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FIGURE 3.1- CONCEPTUAL MODEL OF A SINGLE-ELECTRON 
INDUCED AVALANCHE 
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where v^ is the drift velocity and t is time. The volume 
occupied by this number of electrons is roughly spherical 
(centered at x) and grows in radius r by diffusion as 


r 2 = r 2 + 
o 


6Dt 


(3.3) 


where D is the diffusion coefficient. Thus, the electron 
density as a function of distance is given by 


n e (x) 


e «x(t ) 
4/3Trr 3 (t) 


(3.4) 


while the positive ion density n + may be obtained by in- 
tegrating the equation 


dn 


+ 


dt 


an. 


(3.5) 


Insofar as a continuum description of the ava- 
lanche development is adequate, the above description is 
essentially complete (of course we have only computed 
mean values, the actual description requires solution of 
the drift and diffusion equations). These formulae 
provide the basis with which the phenomena may be ana- 
lytically treated. Extensive analysis for gaseous ava- 
lanches may be found in Raether's book (1960) as well 
as many other places (e.g., Meek and Craggs , 1978). We 
merely note here the availability of this formulation 
which may be called upon as needed. 


The above process continues until either the 
avalanche encounters a conducting boundary, or secondary 
effects become important. The major secondary effects 
ivhich have been observed are those due to attachment and 
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recombination, photoionization and space charge. At- • 
tachment and recombination both tend to quench the 
avalanche. ' They must be accounted for if either effect 
is important at the field strength of interest . Photo- 
ionization of both the ambient gas and the electrodes 
can be caused by photons generated during radiative de- 
excitation of molecules excited by the' drifting electrons. 
While this process is extremely important in gas break- 
down, it has little impact on the evolution of the ava- 
lanche until space charge effects become important. 

The major effect at high amplification is to make the 
avalanche head somewhat more diffuse. 

The major secondary effect which alters the 
development of the avalanche from the above simple pic- 
ture is the space charge effect. The effective drift 
and multiplication of the electrons, together with the 
immobility of the ions leads to charge separation. 

This separated charge has a self -field which grows as 
the amount of free charge in the avalanche grows. When 
this self-field becomes comparable to the applied field, 
the further motion of the electrons is influenced by the 
total- field. Further development requires the self- 
consistent evolution of the field and drift motions. 

While it is relatively easy to give estimates of this 
self-field (Raether, 1960), we will content ourselves 
here with some general observations, as we present a 
numerical prescription and evaluation of the self-field 
in results below. For gas pressures typical of labora- 
tory experimentation (a, 1 Torr) , it has been observed that 
space-charge effects do not become important until the 
avalanche has reached an amplification of 10 D -10 . For 
parameters characteristic of solids, however, one expects 
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self-field effects to become important at much smaller 

o 4 

amplifications, of the order to 10 -10 . Some reflec- 
tion on Figure 3.1 together with inspection of figures 
given in the next section demonstrate that the space 
charge fields are rather complicated. In particular, 
a one dimensional description (as has often been used 
for estimates) is far from adequate to describe its 
features. Qualitatively, the field is that of a dipole, 
so that outside the charge distribution the electric 
field is enhanced over its ambient value, while internal 
to the distribution the field tends to oppose the applied 
field. These features of the space charge field (enhance 
ment outside, decrease inside) .are the characteristics 
'which determine the further evolution of the avalanche. 

It is our working hypothesis' that single- 
electron induced avalanches of the type described above 
•also take place in solids, and that these avalanches 
may ultimately evolve into breakdown channels. In the 
next section we describe our approach to computing the 
properties of these avalanches up until the time that 
space-charge effects become important. Further evolu- 
tion in the ' space-charge dominated regime is described 
in Section 4. 
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3.2 


COMPUTATIONAL APPROACH 


As with the solution- -of the Boltzman transport 
equation, we have opted to simulate the electron ava- 
lanche with a Monte Carlo method. Also as before, a 
'’crude" Monte Carlo has been used. The algorithm is 

(1) The set of distribution functions for the 
location of ionization sites is integrated 
giving cumulative distribution functions 
for ionization sites. 

(2) A 2-D cylindrical (r,z) grid is defined 
and a single electron initiated at 1 point 
in the volume. 

(3) The site of the next ionization is ran- 
domly selected from the cumulative func- 
tion defined in (1). 

(4) The electron is moved to the selected 
position, -the. location of the ionized 
atom and the incident and newly released 
electron are separately recorded. 

(5) The 'newly released electron is added to 
the aggregate of free electrons. 

(6) The process (steps 3-5) is repeated for 
each electron in the current generation. 

(7) After updating the electrons in the cur- 
rent generation, the charge density is . 
calculated. 

(8) The potential at the outer boundary of 
the grid is set via' a multipole expan- 
sion. 

(9) The potential within the volume is solved 
with a relaxation scheme. 

(10) The fields are calculated and if the 
fields are less than an input value, 
the next generation of electrons is 
followed. 
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(11) If the fields exceed the maximum input 
value, the positive ion density and 
electron density are printed and written 
to disk. The net charge density, poten- 
tial and electric fields are also printed. 


The implementation of the above algorithm is 
uncomplicated except for the Poisson solution. Poisson's 
equation in cylindrical coordinates (r,z) (with 0<r<a, 
0<z<L) is 


+ 

3r 2 


1 3(p 
r 3r 



P / £ 


(3.1) 


where p is the charge density and e is the dielectric 
constant of the material. 


A grid can be imposed on the region 0<r<a, 
0<z<L as follows. Let I , J be the number of grid points 
in the r-z directions, respectively, and define 

Az = j— . (3.2) 

The grid coordinates are 


Ar = 


2a 


21-1 


• r = (i-1/2) Ar Z = (j-l)Az ,(3.3) 

J 

for l_<i_<I , Note that r^ = a, Zj = L. Let the 

notation f^j = f(r^,z^) be used and let F = -p/e. The 
differential equation is anchored to the grid by 
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+ 


F. . 


(3.4) 




ij 



1J 


where l<i<I, l<j<J. The potential on the ' 2 = 0 , z=z 

max 

and r=r boundaries is set by a multipole expansion. 


The derivatives are approximated by 


= ^i+1,.1 *1-1,. i 



ij 


jrV 

Sr 2 


9 2 4> 
3z 2 


iJ 


iJ 


2 Ar 


'i+ij 


2d>. . + ®. . , 

y l.T ' 1 - 1,3 


Ar 


»i,J 1 " 2 *i.1, + 
Ar 2 • 


(•3.5) 


(3.6 ) 


(3.7) 


for l<i<I, l<j<J. The radial derivatives at i = 1 are 

derived as follows. Consider the regular grid points 

(r. ,z.) as lying in the 0 plane and let the grid points 
i 3 

(r Q ,z^) lie in the 0 + ir plane, where r Q = -Ar/2. The 
symmetry condition requires that » Kj.<J- If 

these conditions are imposed -for i-=--l the result is 
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3r 


■ 


ifi 

3r 2 . 


2 Ar 


* 2,1 - * 1.1 
Ar 2 


(3.3) 


<3. 9) 


lj 


The above difference framework can be anchored 
to ‘a simultaneous relaxation scheme. Let the residual 
R. . be given by 


R. 


1J 


vV 

3r 2 


J i 14' 

r 3r 


+ 


(3l£' 

l - 2 , 


F. . 
1J 


(3.10) 




1J 


-1J 


where l_<i<I, l<j<J. The relaxation scheme must converge 
which means the residuals should become smaller as the 
iteration process proceeds. Suppose the influence of 
.all poin'ts surrounding the ij th point be neglected and 
let k denote the coefficient of . . Then k is given by 


k = -2 


i + 


Ar 


Az 


(3.11 ) 


and 


R. . = krf> . . - F. . 


(3.12 ) 


Let be a small change in <ju ^ and SR^. the 1 


corresponding change in R . Then 


R. . + SR • • = kb. . + k5d> . . - F. . 
xj xj Y iJ ij 


(3.13) 
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Combining the above yields 


{R ij =■ kS hj 


(3.14) 


The new residual R. . + SR. . is zero provided R. . 

ij ij - ij 


Thus 




— R . . n , 
ij/k 


(3.15) 


Let d>'. . = d>. . + Scb... Then, <j> ! . = <f> . . + cR. where 

■ij • v ij T xj ij • lj 

c - -1/k. 


In order to put the preceding results into an 
iterative scheme let n denote the number of the itera- 
tion, l_<n<N. The simultaneous relaxation scheme can be 
written 


,n+l , n , - n 

♦id = hj + “ cE u 


(3.16) 


where a is a convergence parameter which under relaxes 

or over relaxes the scheme. Note that d>.A , l<i<I, l<j<J, 

' -KJ 

denotes a starting state that must be specified. 

The simultaneous relaxation scheme can be con- 
siderably improved relative to computer run time by 
using all new information immediately after it is com- 
puted. Thus, if the computations for all j are performed 

. . . . # 

for each x as l increases, the quantities 6 . <f>. 

n+l 

are available in the computation for <f> ^ . Such a scheme 
incorporating this new information is a sequential relax- 
ation iteration procedure, and is the scheme used in 
CASCAD . 
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The convergence criteria used for the sequential 

scheme is 


max 

ij 


n+1 

> ij 


, n 

- 4-ij 


V 1J 


< £ 


(3.17) 


Convergence is influenced by starting state-, convergence 
parameter and spatial grid size. 


The above procedure has been implemented in a 
FORTRAN code called CAS CAD . The operational features of 
CASCAD are described in Appendix 5. Some features of 
the code are: 


•- the distributions of first ionization 
sites, calculated by SEMC are required 
as input . 

• rhe size of the grid must be large enough 
so that the multipole expansion of the 
potential on the boundary will converge 
with a few terms (<10). 

• maximum number of generations is limited 
by central memory available (<2l- 4 ). 
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3.3 - RESULTS 

3.3.1 N 2 RESULTS 

Using the parameters calculated by the SEMC 

code, we have calculated the avalanche history through 

12 generations (4096 electrons). As expected., the electric 

field remained small (550 volts/meter). Estimates indicate 

that 20 to 30 generations would be necessary for field to 
7 

approach 10 volts/m. 

3.3.2 TEFLON RESULTS 


Employing the SEMC code results for Teflon, the 
electron avalanche time history in Teflon was computed using 
CASCAD. The computation was terminated when the space charge 
field reached 50% of the applied field. For an applied 

g 

field of 8 x 10 V/m, this occurred after 12 generations. 

This is entirely consistent with the following analytical 
estimate of the space-charge field: 


E 


E 

o 


e^ 

? 3/2 


where t, is the number of generations, and E is related to 
the problem parameters by: 


E 

o 


e 

4-rrs 

o I 



Xj is the ionization length, and r Q is the elementary 
diffusion radius given by: 
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Figure 3.2 shows the net charge density for the 
avalanche after twelve generations. Note its essentially 
dipole character. The positive and negative charge densities 
are shown in Figures 3.3 and 3.4. -Note their very similar 
nature, and the delicate cancellation which gives rise to the 
net charge density of figure 3.2. These charge densities 
are consistent with the estimates given by equations 3.4 
and 3.5. Shown in figure 3.5 is the electric field after 
12 generations in Teflon. Note that it is rather complicated 
but that external to the charge distribution, it is dipole 
in character. These results form the initial conditions for 
the computations of Section 4. No surprises are apparent. 




ALL CHARGE CONTAINED WITH DASHED 30UNDARY . 


REGION 

DENSITY Ccoul/m 3 ) ■ 

SIGN OF CHARGE 

1 

> 5 x 10 5 

+ 

2 ' 

> 1 x IQ 5 

+ 

3 

> 5 x 10 4 

+ , 

4 

> 5 x IQ 5 

- 

5 

> 1 x io 5 

- 

6 

> 5 x 10 4 

- 


FIGURE 3.2. CHARGE DENSITY AS A FUNCTION OF SPACE FOR AN APPLIED 
ELECTRIC FIELD OF 8 x 10 ° volts/m AFTER 12 GENERATIONS 
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Section 4 


NEGATIVE STREAMER DEVELOPMENT 


In the previous section we have developed com- 
putations of the evolution of single electron avalanches 
up until the time when space-charge fields become com- 
parable to the applied field. Further development re- 
quires a description of the self-consistent evolution 
of both the fields and the electronic motions. We 
describe this phenomenon, our computational approach to 
calculating its evolution, and some model results for 
Teflon . 
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4.1 


THEORY 


In Figures 3. 2-3. 5 we have shown the electron, 
hole, total charge and electric field distributions due 
to a single electron induced avalanche in Teflon. We 
see that a negatively charged region appears near the 
front of the avalanche, while a positively charged re- 
gion appears at the tail. In front of the negative 
tip, the field is enhanced over its ambient value. 

Behind the positive tip a similar situation prevails. 
Between the positive and negative charge regions, the 
field is diminished from its impressed value. Further, 
the field begins to develop a significant radial com- 
ponent . 


.In gases [Raether (1960), Meek and Craggs 
(1978)] it has been observed that when the conditions 
of strong space-charge prevail, a different set of 
physical phenomena occur. At both ends of the charge 
distribution the configuration distorts from the 
simple picture of Section 3 above, and a rapid propaga- 
tion phase sets in. One end of the distribution propa- 
gates toward the anode, the other toward the cathode. 
Though the terminology is not completely uniform, these 
are usually called the "anode-directed streamer," and the 
"cathode-directed streamer." The anode-directed streamer 
evolves from the negatively charged tip, while the 
cathode-directed streamer evolves from the positively 
charged tip. It is our working hypothesis that similar 
phenomena can occur in solids. In this section we de- 
velop the theory of negative streamer development. A 
discussion of positive streamer development may be 
found in Section 5 below. 
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For the purposes of this discussion, we will 
assume a single trapping level. We will further assume 
that only electrons can be trapped — that hole trapping 
is not of importance,. The generalization to include 
more complicated carrier kinetics will be apparent. • 

We will further assume that the time scale of interest 
is long enough that all microscopic parameters are in 
local equilibrium, i.e., that macroscopic parameters 
may be used to describe the relevant physics. This 
should be an excellent assumption. For the conceptual 
purposes of isolation of effects, we will also assume 
the dielectric to be locally neutral; with the impressed 
electric field given by some distant sources. 

Let us define the following symbols: 

n^ - number density of conduction "band" electrons 

n + - number density of holes 

n^ - number density of trapped electrons 

Nj_- - number density of trap sites 

p - electronic mobility of conduction band 
electrons 

p + - hole mobility 

- "effective" mobility of trapped electrons 

v e - drift velocity of conduction band electrons 

v + - drift velocity of holes 

v. - "effective" drift velocity of trapped 
electrons 

p - net charge density 
E - Electric field 
$ - Electrostatic Potential 
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With these definitions, let us consider the 
various processes which can alter the basic number - 
densities (n g , n + , n^) [we label this as respec- 
tively] . Each species may be transported by drift 
with the electric field, and diffusion due to number 
density gradients. The basic transport equation (de- 
rivable from Boltzman equation arguments) is 

9n i 

__ + v*J. - sources (4.1) 

d X X 


where 


h = n i f i - D i 5n i < 4 - 2 > 

The quantity J\ is the flux of particles of 
type i, n.v. is the flux due to field drift, -D . vn . is 
the flux due to diffusion, and D, is rhe diffusion co- 
efficient for the i'th species. Each of .the species 
number densities may also change due to sources, i.e., 
due to transitions from one type to the other. The 
number of electrons in the conduction- band may increase 
due to creation of electron-hole pairs (valence band 
ionization), and thermal and collision induced detrap- 
ping, while it may decrease due to trapping, and 
electron-hole recombination. Assuming that all number 
densities are small compared to the valence band number 
density (an excellent assumption), these terms contribute 
sources gn , 6 n, , <$ 0 n n,, -an (N.-n. ), and -bn n,, where 
g is the avalanche rate, 5^. the thermal detrapping rate, 

6 2 the collision induced detrapping rate, a the trapping 
rate, and b the recombination rate. Similarly, the num- 
ber of trapped electrons may 'increase by trapping 
[an (N -n.)l and decrease by detrapping [-6 1 n . -6„n n. ] , 

©XX X X ^ 0 X 
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while the number of holes may increase by pair creation 

[fin ] and decrease by recombination [-bn n ,] . 
e __ e + 

The coefficients in these kinetic source 

terms are in general a function of the electric field 

strength [E[, and the lattice temperature T, as are the 

mobilities and diffusion coefficients. The charge 

density p = e(n + -n -n t ) determines the. space charge 

field through Poisson's equation, while the impressed 

field E is assumed invariant. The lattice temperature 
o 

is determined by the ambient conditions together with 
the joule heating of the lattice J*E, where J is the 
net electrical current . It is assumed that the time 
scale is sufficiently fast that all energy deposition 
is instantaneous, that is, lattice thermal diffusion, 
and mechanical pressure relief are assumed much too 
slow to allow for energy transport away from the depo- 
sition site during the times of interest. The heating 
is therefore a strictly local process, the lattice in- 
ternal energy increasing with the Joule heating, and 
the temperature increasing as a constant volume process . 


With this background, we have the following 
equations describing the evolution of the system: 


3n 

__® + y.J e = £^+(<$.^6^) 


- an e (N t -n t ) 


- bn e n + 


(4.3) 


J = n v -D Vn (4.4) 

e e e e e 
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v = -y E 
e e 


3n, 


— + V- J + = -8n e -bn e n + 


- J + n + v + D + 7n + 


v, = y + E 


3n, 


-3¥ + v ' J t 


an e ( W 


-< 5 l +{i 2 n e )n t 


J t “ n t v t " D t Vn t 


5 t * -’ , t E 


P = e(n + - n e - n t ) 


J - e(J + - J e - J t > 


3_p 

3t 


+ V*J = 0 


V 2 $ = | 


E = -V® + E, 


3T — 1 — * 

C v f = J ' E 


(4.5) 

(4.6) 

(4.7) 

(4.8) 

(4.9) 

(4.10) 

(4.11) 

(4.12) 

( 4 . 13 ) 

(4.14) 

(4.15) 

(4. 16) 

(4.17) 
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is the specific heat at constant .volume for the sub- 
stance . 


A substantial simplification of these equations 
is possible for the conditions of interest. Trap densi- 
ties in most solids are usually at most of the order 
of lol^/cm 3 . As we saw in Section 3 above, -electron 

avalanches quickly reach densities of the order of 
18 3 

10 /cm . The effect of the trapping is therefore ex- 
pected to be of the order of 1% at most. We may there- 
fore neglect trapping kinetics entirely, unless a 
posteriori inspection of the solutions suggests otherwise. 
Formally, we require n t =N t =6 1 =6 2 =a =° ■ ?or the highly 
mobile electrons well up in the conduction band, we also 
have the inequality y + .<<y e . Thus, for the time scales 
of interest (10“-^-10-H secs) we may neglect hole drift. 
This same inequality, together with the Einstein rela- 
tion implies that we may neglect hole diffusion as well. 


In the results presented in Section .4.3 below,- 
the computations also assume that b=0, that is, recom- 
bination has been ignored. This is justified a posteriori 
by noting that the maximum 'hole densities reached are of 
the order of 10 19 /cm 3 . The recombination coefficient for 
Teflon is thought to be ~ 10 -9 cm 3 /sec [Gross (1978)], 

so that the maximum value of this term is . For 

e 

the time scales of the calculations ( *> 10 -1 ^) this term 
may be neglected. For the larger times associated with 
the discussion of Section 5, this neglect is no longer 
justified. The subject will be readdressed during that 
discussion. 
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The above set of equations, together with the 
simplifications provide the basic model for the descrip- 
tion of the development of the negative streamer. These 
equations should be good if: 

a) there are sufficient numbers of par- 
ticles available that a continuum de- 
scription of the particle properties 
is justified 

b) the time scale is sufficiently long 
(>10~13 sec) that microscopic equi- 
librium has been achieved, and 

c) the time scale is sufficiently short 
(>10-10 sec) that the conditions on 
trapping and recombination noted 
above are valid. 

All the above conditions are thought to be satisfied for 
computations given below. 

What is the behavior of such a streamer ex- 
pected to be? The literature associated with the compu- 
tational modeling of this phenomena is fairly limited, and 
has been limited to investigations of gaseous discharges. 
Further, most of the work has been either one-dimensional, 
or quasi-two dimensional. Only the recent work of Davies, 
Evans and Woodson (1978) and Geary and Penney (197-8) [see 
also Geary (1973)], appear to treat the fully two- 
dimensional problem. Earlier work may be traced from these 
references, as well as Bayle and Bayle (1974), Kline (1975), 
Ward (1965), Ward (1976), Kline and Sianbis (1971), 
Reininghaus (1973), Yoshida and Tagashira (1976) and 
Meek and Cragg (1978). These results are a bit diffi- 
cult- to interpret because of complications due to photo- 
ionization, parallel plate circuitry, secondary cathode 
ionization processes and an emphasis on the cathode 
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directed streamer. The results suggest that a thin 
shell of space charge develops which propagates at 
essentially a uniform velocity, and that an essentially' 
neutral conducting region is left behind. The space 
charge distortion of the field reaches an equilibrium 
at a field substantially higher than the impressed 
field, and the effective velocity of propagation is con- 
sequently larger than would be expected from drift in 
the ambient field. Features such as these are best 
seen in the paper of Geary and Penney (1978). We will 
not discuss these qualitative features further here, 
as numerical results will be presented in Section 4.3 
below. We merely note that it is the propagation of 
an ionizing space charge wave which we are attempting 
to model. We concentrate in this section on the nega- 
tive tip alone. 
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COMPUTATIONAL APPROACH 


4.2 


In this section, we write down the specific 
equations for streamer propagation and describe the 
numerical methods used in the solutions. The system under 
consideration is axially symmetric and the equations are 
solved in 2-D cylindrical coordinates. A description of the 
operational feature of the code is given in Appendix 6. 

The code is an adaptation of the code SPARK2 which is 
publicly available from the CPC Program Library of the 
Queen's University of Belfast. This code has been documented 
by Davis, et al . (1978) in Computer Physics Communications . 
The reader is referred to this primary reference for a 
detailed discussion of several of the major algorithms. 

Our version of the code is called ACORN. The equation of 
basic interest is that for electron charge continuity and 
is obtained by combination of equations (4.3) and (4.4): 



V 


n v 
e e 


Bn 


+ D V 2 n o 
e e 


(4.18) 


To complete the set of equations, we require equation (4.5) 
and specification u , 8, and D , The space charge part of 
the electric field is formally given in terms of a potential 
(equation 4.16), which, in turn, is obtained through 
solution of equation (4.15), the Poisson equation. The net 
charge density p in this equation can be obtained equivalently 
in two ways. We may solve eq. (4.6), along with eq. (4.3), 
to obtain the net density given by equation (4.12). Alternately, 
we can solve the net charge continuity equation. (4 . 14) directly. 
The latter is the procedure we have taken in the numerical 
algorithm. The expression for the positive density can then 
be found from equation (4.12). 
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In summary, the equations ancillary to equation 
(4.18) are:. 


V = ■ -U E 
e *e 


p = e(n + - n e ) 


J ~ -eJ = en V + eD Vn 
e e e • e e 


l£ 

3t 


+ V * J = 0 


■V 2 f = £ 

£ 


E — — V<f> t* E 


(4.19) 

(4.20) 
(4. .21) 

(4.22) 

(4.23) 

(4.24) 


The derivations for y 0, and D have been 

0 © 

discussed previously. In principle, the equations can 
admit arbitrary functions of the electric field. We have 
found for the Teflon case that y g and D g are nearly inde- 
pendent of the field, and we have programmed simple constants 
into the code. On the other hand, 0 is found to take on the 
form 


S = a 0 I^el exp( ~ a E / l E l ) 

where a Q and ctg are constants. 
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Equation (4.18) is solved numerically by a time- 
splitting process. The diffusion part of the equation is 
first stepped forward over a trial time step. The diffusion 
equation is solved implicitly so that the result is stable. 

The Peacemen-Eachf ord ADI method was used to solve the 
equation. The advantage of the scheme is that the diffusion 
equation can be written in tridiagonal form, thus simplifying 
the matrix inversion required for the solution. The remaining 
part of equation (4.18) is solved in the last frame of reference 
moving at the electron drift velocity, that is along the 
characteristic curve. This’ concept implies that the velocity 
dtiring the time interval for the calculation must be known, 
and this is determined by an iterative scheme. At the beginning 
of the iteration, the drift velocity corresponding to the 
electric field found at the end of the previous time interval 
is used to estimate the characteristic derivative, during the 
trial time interval. Then the continuity equation is inte- 
grated along the characteristic. The net change density- 
equation is then found by an implicit solution of equation 
(4.21). Once the charge density is determined, Poisson's 
equation is solved to find the potential from which the elec- 
tric field and corresponding electron drift velocities are 
obtained. These new values of the electron drift velocities 
are fed back into the beginning of the iteration sequence to 
find better estimates for the characteric curves. The solution 
along characteristic curves, as- well as the diffusion solution, 
are calculated at least twice. Convergence tests are employed 
after solution of the potentials and after determination of 
the electric fields. A failure in these tests means that the 
solution along characteristics, as well as the time split 
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diffusion solution ,- must undergo at least two more itera- 
tions. When appropriate, a time step reduction is made before 
a new set of iterations on the characteristics is carried out. 

The 'Poisson equation is solved in the z-direction 
by a method similar to Hockney's fast Fourier transform 
method. The remaining ordinary differential equations in 
the r-direction depend on the Fourier variable K and are 
solved by standard methods. The prescribed boundary conditions 
.are introduced into the equations as pseudo charge sources. 

The boundary potentials for the streamer problem 
correspond to free boundaries enclosing a charge distribution. 
With each iteration step where a new 'Charge distribution is 
computed, a new set of boundary potentials is also determined. 
The boundaries are placed far away from the region of charge, 
so that a multipole expansion is appropriate. This technique 
is especially valuable for axisymmetric charge distributions 
because the azimuthal angle integration can be carried out 
analytically. There remain a set of two-dimensional charge 
moment integrals to solve. In order to make the resulting 
series converge as fast as possible, the origin for the moment 
calculations- is taken to be midway between the positions of 
positive and negative charge centers, respectively. 

Convergence of the iterative scheme was determined 
through two major tests. The first test was that no "fluid" 
element would drift more than halfway across a cell in one 
time step. This is the standard Courant check required for 
the stability of explicit solutions of the continuity equation. 
This’ technique is used to ensure that the trial time interval 
at the beginning of the iteration is indeed sufficiently small'. 
Another check requires that the change of electric field over 


4-13 



the time interval did not exceed a prescribed value. This 
criterion is a measure of the conductivity of the system, 
for an electric field relaxes very quickly in a highly 
conductive medium. Since the conductivity is easily 
shown to be proportional to. the electron density, the 
highest electron density in the system determines the 
smallest time step according to this criterion. In prac- 
tice, this test was conservative and gave a smaller time 
step then was actually necessary, and a- different stabil- 
ity criterion described below has been implemented. The. 
reason that this criterion is overly strict is that the 
highest electron number densities occur behind the head 
of the anode' streamer in a "passive''' region where electric 
fields are small and near charge neutral conditions pre- 
vail. Under these conditions, .some overshoot of the 
electric field relaxation can be 'permitted without harm- 
ing the overall streamer solution. It is. important, 
however, that only a fractional change of the electric 
field be allowed during a time step in the streamer head 
where the net charge density is negative. Since the 
electron density is usually smaller here than in the 
neutral region., the time increment will be larger when 
the relaxation criterion is applied only to the head of 
the streamers. 
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4.3 


COMPUTATIONAL RESULTS 


We have exercised ACORN to compute the evolution 
of the electron avalanche in Teflon when self-electric 
fields are important. The' computation was performed for an 

g 

applied field of 8 x 10 V/m. The parameters used for the 
electron drift quantities were those reported in Section 

2.3.2. The initial conditions were those taken from Section 

3.3.2. The computation was performed on a 21 x 64 grid with 

-9 ° 

both radial and axial grid sizes of about 7.5 x 10 m (75 A). 

The computations were -carried out for a total evolved time of 
0.73 ps. By this time in the calculation, the initial distribu- 
tion has reached very near the boundary of the grid, and the 
required time steps have become quite small. Longer time 
computations have been performed (^2 ps) using slightly different 
material parameters on a grid with twice the above spacing. 

Though these results (meant to represent Lucite) are for a 
different material, they indicate that the results given above 
are indicative of the general features of the evolution.. We 
will present these results elsewhere, and restrict our attention 
to the above noted Teflon computation. 

Shown in Figure. 4.1 and 4.2 are the time evolution 
5 3 

of the high density ( >10 C/m ) regions of electron and hole 
concentrations.- Note that the radial scale on these (and 
many subsequent) figures has been expanded to better display 
this variation. This region can be seen to be increasing 
in size by drift with the field and diffusion radially. A 
somewhat better feel for the spatial dependence of these 
densities may be seen by considering the contour plots of 
these quantities at the end of the calculation shown in 
figures 4.3 and 4, .4. Note the sharp variation in profile near 
the axis in the direction of propagation. A better feel for 
the propagation can be obtained by considering the regions of 
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high net charge density. Shown in figure 4.6 is the region 

5 3 

of Targe net negative charge density (>10 -C/m ) , as a function 
of time. This region shows clearly the propagation of this 
"active" region to the right. This region propagates to the 
right and grows radially by diffusion and field distortion. 

The growth of the positive charge region which is left 
behind is shown in figure 4.7. Note its growth in both the 
radial and axial directions. Between the two regions of 
positive and negative .charge, an essentially neutral region 
develops which grows as the active space charge head propa- 
gates to the right. This may be seen in the net charge 
density contour plot of figure 4.5, Note particularly the 
sharp variations in net charge density in the axial direction, 
in rough conformance with the results' given in Appendix 8. An 
intuitive feel for the evolution of the field and charge 
density may be obtained from figure 4.8 which shows both the 
field and high charge density regions at the end of the 
calculation.. Note the enhancement of the field over the 
ambient field both near the negative charge region and the 
positive charge region. In the internal essentially neutral 
region between the two charge regions, the field has been 
relaxed to a small value. External to this region, the field 
points either toward or away from the region, i.e. , the fields 
are like those around 1 a conductor. This type of feature is 
continuously evolving in t.he calculation. 

A quantitative feel for these results will be 
obtained from figure 4.9, which shows the axial electric 
field and net charge density as a function of position 
along the axis at the end of the calculation. The enhanced 
fields at either end are clearly associated with the exposed 
charge regions. Note a factor of two increase in field near 
•the positive tip, with not quite so large an increase near 
the negative tip. The time evolution of the field along the 
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axis is shown in figure 4.10. Note the clear propagation 
of the negative head, and the continuing increase in enhanced 
field. 


Some feeling for the variation in the radial 
field may be obtained from figure 4.11, which shows this 
variation on two fixed radii as a function of axial distance. 
Note the correlation with the exposed charge regions. The ra- 
dial variation of the radial field may be- seen in figure 4.12, 
where several, fixed Z profiles are displayed. The profile 
at Z=0.1 y coincides with the positive charge region, that at 
0.3 y with the negative charge region,’ while the Z=0.2 y 
profile -falls in between. 

The peak value of the axial electric field on axis 
is shown in figure 4.13 as a function of time. It is appar- 
ent that this is not nearly saturated. Ohmic heating of 
the lattice gives rise to a temperature change of the solid 
material in the avalanching regions. Shown in figure 4.14 
is the maximum temperature a-s a function of time. While this 
does not represent a temperature rise nearly sufficient 
to cause major alteration of the. solid, it does represent a 
significant temperature rise for so small a period of time. 

A profile of the temperature change is shown in figure -4.15. 
Note the sharp variation near the avalanche head. 

Shown in figure 4.16 is the axial centroid of 
electron density as a function of time. While the data is not 
accurate enough to determine adequately the acceleration of 
the head, the mean velocity from the slope of 1.65 x 10^ m/s 
is larger than the drift velocity in the ambient field. 

Another useful diagnostic on the propagation is given in 
figure 4.17, which gives the length of the effective dipole 



as a function of time. The calculated results indicate 
an expanding dipole with acceleration. The strength of 
that same dipole can be seen as a function of time in 
figure 4.18. Note its rapid growth. Finally, the total 
number of electrons involved in the entire process is 
shown in figure 4.19. Note that this number is signifi- 
cantly smaller than that which would be obtained by un- 
retarded avalanche. 
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FIGURE 4.1. REGION FOR WHICH ELECTRON DENSITY > IQ 5 coul/m 3 IN TEFLON 
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ELECTRON NUMBER DENSITY CONTOURS AT .73 ps IN TEFLON' 


FIGURE 4.3. 
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FIGURE 4.5. 


NET CHARGE DENSITY CONTOURS AT .73 ps IN TEFLON 
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FIGURE 4.7. 


PROPAGATION OF THE REGION OF- POSITIVE 
CHARGE DENSITY > 10 5 coul /m^ in TEFLON 
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FIGURE 4.8. VECTOR FIELD 
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FIGURE 4.9. AXIAL FIELD AND NET CHARGE DENSITY 
ON AXIS TEFLON AT t = .73 ps 
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FIGURE 4.15. T 



FIGURE 4.16. POSITION OF ELECTRON CLOUD 
CENTER VS. TIME-TEFLON . 













Section 5 


SYNTHESIS 


In the previous three sections, we .have described 
the computational features of the basic elements of our 
model of dielectric breakdown in electron-irradiated solids. 
Basically, we have provided a mathematical description of 
what we perceive as the initiation stage of the breakdown 
development. It is the purpose of this section to attempt 
to interpret these results and to provide a conceptual model 
of the continued development and propagation of the discharge 
through tree formation. 


We consider first the development of the negative 
streamer. This phenomenon has the following features: 


(1) Drifting, avalanching electrons give rise 
to a space charge separation which leads to 
enhanced electric fields, at the head of the 
avalanching region. 

(2) The enhanced electric field gives rise to an 
enhanced avalanche rate through the field 
dependence of avalanche coefficients. This 
provides further enhancement of the field. 

(3) The avalanche production of mobile electrons 
provides su-fficient carrier number density 
to permit the relaxation of the electric 
fields internal to the streamer channel in 
very short times (S.10~ 13 sec). This region 

of the solid locally appears to be a conductor. 

(4) The acceleration process at the' head of the 
streamer saturates. A space-charge ionization 
wave evolves which propagates with constant 
velocity in an essentially self-similar 
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fashion. This saturation velocity, together 
with the spatial character of the field and 
electron density in the front appear to be 
roughly independent of the initial conditions 
of the streamer development. That is, once a 
critical number density is reached in an .ava- 
lanche, such that self-fields are important, 
the further development is driven by these 
self-fields . 

(5) Very large fields are developed near the 
positive end of the streamer. 


Strictly speaking, not all of the above items may 
be taken as conclusions from the. results which have been 
computed to date. Rather, they represent extrapolations of 
the results based on reasonable theoretical -expectations . 
Further work will be required to computationally verify these 
assertions. For now, we content ourselves with various 
theoretical arguments in support of these expectations. 

The assumption that the electric field- is derivable 
from an electrostatic potential means that V x E = 0. This, 
in turn, means that the magnetic field is negligible. From 
Ampere's law, with V x B. = 0., we find: 



where J is the current density. CRelation (5.1.) should not be 
taken too literally, since it would imply that the field can 
only change in the conducting region — It is useful for an 
estimate, howeverl . The current density J is given by 
Ohm's law, J = aE, with the conductivity, o, given by: 

a = n e ev e (5.2) 
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where n g is the electron nvimber density, e the electronic 
charge, and the -electronic mobility. Relation (5.1) then 
represents the usual relaxation equation with relation time 
t given by: 
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r 


£ 

a 


n e e ^e 


(5.3) 


Choosing the mobility to be roughly constant 
—4 2 

(V e ~ 1.6 x 10 m /v-sec), as computed for Teflon in 

Section 2 above, we can see that this relaxation time 

depends only on the electron carrier density. A relaxation 
-13 

time of 10 secs is achieved for an electronic charge 

ft o 1 R R 

density of 1.1 x 10 C/m (n g =-6.9 x 10 electrons/cm ) . 
This relaxation time is the dominant time of the problem, 
and should be kept in mind’ during the following discussion. 
Of course, this time depends critically on the density of 
electrons which are produced.' 

A reasonable estimate of the maximum value for 
this quantity may be developed. Denoting the convective 
time derivative by D/Dt (D/Dt = 3/3t + v • V), eq. (4.18) 
for the electron dynamics may be written: 

* (D7n e ),(5.4) 

Ignoring 1 the- diffusion term, we expect the electron density 
to level out when: 

g - V ♦ v e ss 0 (5.5) 
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Asstiming a roughly constant mobility, using the relation 
£ = a|v D |, where a is the first Townsend coefficient, and 
invoking Gauss's law, eq. ('5 .5) may be rewritten as-: 


P MAX 


eg 

U 


eEa 


(5.6) ■ 


Referring to Figure 2.6, we can see that this is a 
rather sharp function of the field strength. Taking 

the values indicated for the ambient field, this relation 

. * 6 
.gives a value of p 1.1 x 10 . 

This number may be expected to be within a factor 
of two to three of that of the individual particle densities 
(i.e., p is achieved in a region where p + ~0). These estimates 
are entirely consistent with the numerical results. 

The enhancement of the electric field ar the • - 

negative streamer tip is apparent from the figures of Section 
4. This enhancement is precisely that which may be expected 
from a dipole charge distribution and is similar to that seen 
in the primary avalanche of Section 3, The results of 
Section 4 also indicate that this enhanced electric field is 
accelerating the avalanching process' that the region of net 
negative charge is moving with a velocity faster than would 
be expected from the ambient electric field. What we see 
from Figure 4.6 is that the region of net negative electron 
density is propagating in an -essentially self -similar (in the ■ 
dZ direction) fashion, with .radial growth due to diffusion. 
While some acceleration of this region is apparent from our 
results, the numerical computations have not been carried far 
enough to see an expected saturation of the tip velocity. We 
expect saturation near the tip for the following reasons. 
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The inverse ionization length (first Townsend coefficient) 

takes on a minimum value. The results of Section 2 for 

Teflon may be. fit in the usual form a - a exp(-£ /E), with 
o_l o o 

a ~1.85 x 10° m . This corresponds to the fact that at 

very high fields, the distance to achieve ionization must 

be at least (Ncr^) where N is an effective valence band 

electron -number density, and Oj is an effective ionization 

cross section. Thus,, the distance between ionizations tends’ 

to a minimum value as the field increases. Similarly, the 

electron drift velocity tends to reach a maximum. This is 

given roughly by assuming no collisions other than ionization. 

'The mean velocity is then roughly 1/2 that required to give 

an energy equal to the ionization threshold. For Teflon, we 

therefore expect the drift velocity to saturate at a value of 

5 

about 7 x 10 m/sec. 

The above two considerations imply that the 

ionization rate 8 = ctjv^l also saturates. There is therefore 

a maximum 'electric field which could be sustained by this 

avalanching process, and is given roughly by the field at 

which the above quantities saturate. For Teflon, we estimate 

9 

this field to -be about 6 x 10 v/m. In practice, the effects 
of diffusion and multiple dimensions may be expected to 
lower this by a's much as a factor of two. Nonetheless, the 
basic idea for the expected saturation is clear. We do not 
expect a completely runaway process. We do not have a good 
estimate at this time for the amount of propagation time 
required to achieve the expected saturation. 

The important aspects of this negative streamer are 
that it’ propagates and is self-sustaining. That is, if a low 
probability event gives rise to an avalanche of sufficient 
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length that space charge effects become important, these 
effects will tend to maintain the effect through further 
generations. Further, the effects lead to an enhanced field 
which makes the process self-sustaining . Assuming the break- 
down is initiated by the above streamer (to .be discussed 
below), this phenomenon would explain why breakdown- which 
is initiated at edges and weak spots at fields that are 
below bulk breakdown fields nonetheless can propagate into 
the bulk. That is, the low field initiated weak spot break- 
down does not remain confined, but, in fact, can propagate 
throughout the region which is presumably not weak. We will 

return to this point below. Note the propagation velocity 

5 

of the negative streamer (^2 x 10- m/sec) is sufficiently 
large that it can cover small macroscopic distances 
(1 mil 'v 25p) in very short times (125 ps ) . 

This brings us to the next important feature of 
the above process. We are not suggesting that the propaga- 
tion of the above negative streamer ij; the breakdown process. 
In fact, the rather small changes in lattice temperature 
rioted in Section 4 above indicate that this is distinctly 
not the case. If the above developed streamer 
propagated to- a conducting- boundary, and no further source 
of current were made available, then the net effect of the 
process would be merely a rearrangement of electric field. 
Internal to the channel, recombination of the mobile carriers 
would quickly decrease the conductivity, and the process would 
be at end. For breakdown to occur, further processes must 
be at work on a timescale more rapid than the recombination 
•process. To understand what these processes might be, we 
must shift our attention to the field rearrangement nature 
of the process, and concentrate our attention 'on the immobile 
positive tip of the streamer. 
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The region of net positive charge left behind occurs 
roughly on the outside and rear of the conducting central 
region. The field near this region is substantially per- 
turbed from its ambient value. In particular, the field is 
enhanced and is directed toward the conducting region. The 
positive charge density falls off roughly as given by the 
considerations of Appendix 9. It is our hypothesis that the 
negative streamer process which, has been modeled above now 
repeats, this time in the presence of the distorted positive 
tip field. -To account for this occurrence, we must now 
take into account the discrete probabilistic nature of the 
process. In order for an avalanche and negative streamer 
to form, not only must the field be sufficiently large, but 
an initiating electron must be available in the conduction 
band. Let us estimate this probability. Consider a roughly 
spherical volume around the positive tip with radius of the 

_7 

order of 10 m (region of field enhancement). The volume 

—21 3 

of this region of space is about 4 x 10 m . In regions 

16 3 

of dielectric which are uncharged, at most c,10 e/m are 

expected to be found in the conduction band EGross (1978)1, 

-5 

so that on the average, only 4 x 10 electrons are expected 

to be found in this volume. In other words, the probability 

of finding an electron in this volume is' quite small. However, 

in regions of the dielectric which have accumulated excess 

charge due to electron bombardment, the excess charge number 

21 22 3 

density can be as high as 10 --10 e/m CBeers, et al . 

(1979)1. In regions such as this, the number of excess 
electrons to be found in the positive tip field region is on 
the average, from 4 to 40. Assuming that, on the time 
scales of interest (1-10 ps), one of these trapped 
electrons makes an excursion into the conduction -band either 
via thermal processes, field assisted thermal processes 
(Poole-Frenkel) , or direct beam processes, then the enhanced 
field in this region essentially guarantees that this, elec- 
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electron will undergo avalanche in the direction of the local 
field. This direction is no longer determined by the ambient 
field. 


The appearance of mobile electrons due to this 
second avalanche allows the charge to rearrange precisely as 
given in our above picture of the primary avalanche and 
streamer development. The net effect of the second process 
is to transport the positive tip in the opposite direction 
of the avalanche. That is, the positive tip appears to 
propagate. (It should be noted that probability arguments 
of the type invoked above have led researchers in gas breakdown 
to include photoionization processes in their description of 
positive point Ccathode-directed streamer! propagation, this 
being necessary to account for the carriers necessary to 
permit the propagation. This photoionization process appears 
to be unnecessary to account for the behavior of positive 
point propagation in electron- irradiated solid dielectrics.) 

The propagation velocity is just that of the electron drift 
in the enhanced field. Repeated occurrence of the above 
process leads to a large scale propagation of the positive 
tip disturbance. Assuming that the controlling process is 
not the availability of avalanching electrons (the above 
estimates suggest it is not), then the positive tip also 
propagates at a velocity similar to the negative tip. 

We believe that this process is intrinsically 
discrete and probabilistic. Thus, the direction that 
secondary avalanches take is determined by the probability 
of finding an electron at the location associated with the 
given field direction. Thus, the direction of the propagation 
no longer needs to follow the applied field, but is determined 
both by the field direction near the positive tip, and the 
density of available carriers in different directions. 
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Needless to say, this will provide for the preferential 
propagation of the positive tip into regions of high 
electron density. It is also quite clear that the proba- 
bilistic nature of the process is conducive to bifurcation 
i.e, treeing. In regions of high electron density, independent 
avalanches may be expected to develop at nearly concurrent 
times. These processes can continue without interfering 
until their space charge fields overlap and influence one 
another. Branching, and direction change, are thus intrinsic 
feat\ires of the model. 

Insofar as the tree propagation also coincides 
with the regions of net excess charge density, each new 
avalanche/streamer contributes a net charge transfer out of 
the region. This net charge will be preferentially trans- 
ferred 'along the entire conducting channel to the negative 
tip. This net charge transfer is important to the continued 
development of the breakdown, as it corresponds to the genera- 
tion of net current to the system, thus extracting energy 
from the stored charge field. 

It is clear how a computer model of the above 
process can be generated. A model of -the trapped charge • 
distribution is first developed. Then a probability function 
is ascertained for determining the probability of finding an 
electron in the conduction band during a sampling interval. 

This probability distribution is used to develop a sampling 
algorithm for choosing the location of the appearance - of 
the next electron during the sampling interval. Having made 
the choice of the location of .-the next available electron, a 
stored description of the positive tip field developed in 
Section 4 above is consulted. The direction and magnitude 
of the ensuing avalanche/streamer from initiation to the 
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existing channel is determined. The positive tip field is 
transferred -to the new initiation location. This process 
is repeated sequentially to generate many channel segments 
of a tree. The radius of each tree segment "is updated at 
each time step to account for radial diffusion. The net 
charge transfer through each of the branches and trunk is • 
scored. Ohmic heating of the laxtice is determined in each 
channel segment. The overall driving field is relaxed 
globally according to external current transfer and circuit 
equation restrictions. Such a model would provide a 
reasonably complete description of the entire discharge 
process . 


Before describing the predictions of such a model, 

we make one further note about Ohmic heating. On the time- 

scales of concern to this problem, lattice thermo-mechanical 

energy transfer is unimportant — these processes are far too 

slow. With the heating rates noted below, if the primary 

discharge streamer propagates to a conducting boundary, then 

subsequent positive poinx propagation will result in net 

current flow through the channel. For the field treated 

20 3 

herein, the Ohmic heating rate is about 10 watts/m . This 
corresponds to a 1000° K temperature rise in 20-30 ps. Clearly, 
this process cannot go on very long before our assumptions 
about the electronic mobility are no longer valid, i.e., 
the electronic transport model given above depends on the 
lattice temperature and structure. We believe that a rough 
estimate of the final temperature may be obtained by assuming 
that the electronic and lattice temperatures come to rough 
equilibrium. Assuming the mean energy of the electrons for 
this lattice temperature is not severely different from the 
mean electron energy in the cold solid, we may equate this 
final temperature to this preheating mean electron energy. 
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For the Teflon calculations of .Section 2, the mean electron 
energy is approximately 0.8 eV. This corresponds to a tempera- 
ture of about 6000° K. ffe believe this to be an upper bound . 
on the temperature, and should be correct within a factor of 
two . 

• Based on this model, we make the following predic- 
tions of the model for the case of electron-beam irradiated 
Teflon. The sample is assumed to be circular, of radius R, 
and the beam deposition is assumed to be thin compared to the 
sample thickness, or R. 

(1) Peak mobile electron density: 
p e s 10 6 C/m 3 

(2) Peak current density: 

J = 10 x A/.rn 

(3) Channel temperature: 3000 - 6000° K. 

(4) Roughly linear current rise in time with 
slope of 6 ttDJ , or 

I = '8 x 10 8 t. 

(5) Peak current proportional to sample radius 
•with slope of 12frDJ/v D or 

J peak ^ 8 x X ° 3 R( "° 

(6) Discharge pulse rise time given by the 
discharge length over the drift velocity or 

x ~ 10" 5 R(m) 
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(7) Discharge tree in the plane of charge 
deposition. 

(8) Maximum channel radius r proportional 
to the square root of th§ sample radius 

r = 2 x 10" 4 R 1 ^ 2 (m) 
o 
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Appendix 1 


N 2 MODEL 


The scattering and energy loss parameters 
needed by SEMC will be presented in this Appendix. 

SEMC is a single scattering Monte Carlo code and thus 
requires information on individual electron interaction 
processes in contrast to most Monte Carlo electron trans- 
port codes which use the continuous slowing down and 
small angle scattering approximations and thus require 
only a stopping power and multiple scattering formula. 

The interaction parameters needed by SEMC are the elastic 
scattering cross section, total excitation cross section, 
mean excitation energy and total ionization cross sec- 
tion. These will follow for molecular nitrogen (Ng) 
which was chosen for validating the newly developed 
code . 


The elastic scattering cross section appears 
in Figure (A. 1.1) versus the incident electron energy 
(see Strickland, et al. (1976) and Jasperse (1976) and 
references therein). Its differential form is assumed 
isotropic in angle which should be acceptable in the 
energy range of interest (< 10 eV). 

The excitation cross section is given in 
Figure (A. 1.2). It contains contributions for several 
excitation processes starting with vibration excitation 
(which produces the feature peaking at 2.5 eV) and 
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extending through excitation of Rydberg states whose 
excitation energies are above 15 eV. The results shown 
above 8 eV are based on the individual cross sections 
given by Strickland et al. (1976), The vibrational ex- 
citation cross section was taken from Jasperse (1976). 

The differential form of the cross section is given by a 
delta function in angle and thus excitation does not 
lead to a change in direction of the electron. 

The mean excitation energy is given in 
Figure (A. 1.3). It is obtained from the formula 

W(E) = E W, a. (E) / E a. (E) eV 
j 3 J ' 3 J 

where the sums are over excitation processes (ionization 
excluded) , W. is the j th excitation energy, and cr. is the 

J J 

excitation cross section. The denominator gives the 
cross section appearing in Figure (A.1.2). 

The total ionization cross section is given in 
Figure (A. 1.4) and was taken from Strickland et al . (1976). 
This cross section serves to remove the electron from fur- 
ther transport. Thus, considerations regarding its dif- 
ferential form i'n energy loss (equivantly secondary 
electron energy) and angle are not required. The sum of 
all cross sections discussed above appears in 
Figure (A. 1.5). This concludes the presentation of the 
atomic information needed by SEMC. 
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Appendix 2 


TEFLON MODEL 


The material model for Teflon (CF 0 , with 
3 

density 2.15 g/em and dieletric Constant 2.0) is de- 
scribed in this Appendix. In treating the electron 
multiplication in a material, the basic parameters con- 
sist of (a) the mean free path for acoustic phonon 
scattering, (b) the mean free path for optical phonon 
emission, (c) the quantum energy of the optical phonon, 
and (d) the threshold energy of pair production and the 
corresponding cross sections. We will describe briefly 
in the following the simple approaches we took in 
achieving this Teflon model. 

A. 2.1 ELECTRON MEAN FREE PATH FOR ACOUSTIC PHONON 
SCATTERING 


We studied this quantity from two approaches: 

Solid State Approach — The inverse mean free 
path of slow electrons for acoustic phonon scattering 
is given (Seitz, 1948; Ziman , I960; O' Dwyer, 1973) by 




kTm* 2 

?l 4 s 2 MN 


(A. 2.1) 


where C represents the coupling constant, in units of 
eV, between the electron and the acoustic phonon, k is 
the Boltzman constant, T temperature, m* the effective 
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mass of electron, % Planck's constant, s sound speed of 
Teflon, M Is the reduced mass of carbon and fluorine 
and N the number of unit cells per unit volume. The 
values used in Eq. (A. 2.1) for Teflon are as follows: 

C = 14 eV 

T = 300°K 

m* = free electron mass 
5 

s = 2x10 cm/sec 

M =7.2 times mass of proton 

_3 

N = (3 times Bohr radius) 

Eq. (.A. 2.1) then gives X = 68. 

The choice of C=14 eV for Teflon is based on 
the general trend of C reported in literature (Seitz, 
1948). C=31 eV for diamond, 10 eV for silicon and 

5 eV for germanium. Kittel (1963) also used C 30. eV 
and s=5xl0 cm/sec in discussing the phonon cloud 
around a slow electron. 

The value for N (corresponding to a bond 
length of 1.58) is derived from a study of the Raman 
spectroscopy of polymer DUDF-PTFE (Latour, 1977). 

Atomic Approach — If we consider Teflon as a 
gas with density the same as that of the solid, then the 
elastic scattering in the gaseous state will correspond tc 
acoustic phonon scattering in the solid state. V/e have 
carried out an analysis for elastic scattering mean 
free path based on a screened Rutherford scattering 
model. The elastic cross section for an electron with 
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V(r) 


-Ze 

T 


(A. 2. 2) 


2 


e 


-br 


is given, in atomic units, by 


da 



__1 

C(l-cos 6) + 2n Q J 2 


where 



(A. 2. 3) 


(A. 2. 4) 


For the Thomas-Fermi model of atomic systems, 

b = Z 1 / 3 /0.885a o (A. 2. 5) 

where a is the Bohr radius, 
o 

The solid state effect on this atomic approacn 

has been discussed (Nigam et al . , 1959) and its -result is 

2 

effectively to replace n Q by n = n Q U ’where p=1.80. The di 
mensionless r) therefore becomes, with E in atomic units 
(27.21 eV is one atomic unit of energy), 

2/3 

n = 0.517 (A. 2. 6) 

A slow electron of a few eV in a moderate Z 
material will have a very large n compared to the 
(1-cos 6) factor i-n Eq. (A. 2. 3). The elastic scattering 
is therefore quite isotropic for slow electrons. Ne- 
glecting the (1-cos 6), and integrating over the angles, 
one reaches the following formula for the total elastic 
scattering cross section 
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a = 2.94 Z 


(A. 2. 7) 


2/3 2 

A teflon gas (CF„) of density 2.15 g/cm 3 con- 

09 " 0 9 

tains 2.57x10 carbon atoms and 5.15x10 ^ flourine 
atoms per cm . The mean free path of electron for 
elastic scattering in this gas is computed by Eq . (A. 2. 7) 
to be 48 compared to 68 from the solid state approach. 

The agreement of the results from these two 
approaches is considered very satisfactory in the sense 
that many uncertainties are involved in either approach 
and yet the main physics in each approach is accurately 
maintained. In the calculation, we use 68. 
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A. 2.2 


ELECTRON MEAN FREE PATH FOR OPTICAL PHONON 
EMISSION 


A detailed treatment for the mobility of elec- 
trons in pure non-polar insulators has been given by 
Seitz (1948). Besides the acoustic phonon, the non-acoustic 
phonons (optical phonons) were discussed according to 
the degree of symmetry in the insulator structure. 

For teflon, which we assume to possess minimal symmetry, 
the ratio of cross sections of the optical phonon emis- 
sion to the acoustic phonon scattering is given by 
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(A. 2. 6) 


In this equation, E is the energy of the slow electron 
capable of emitting an optical phonon of quantum energy 
/zo) (E>7iw-)* D and C, being almost equal (Seitz 1 , 1948; 
Shockley, 1951), are the coupling constants of electrons 
with optical phonons and acoustic phonons respectively. 

K is the .maximum wave number in the Brillouin zone. 


For Teflon, the energy of the optical phonon, 

as will be seen in the next section, is quite large 

compared to that in semiconductor or diamond. It is 

therefore not suitable to put TzKs'vpRo as Seitz (1948) did 

5 

for diamond. In fact, using s=2xl0 cm/sec and K=7r/3a o , 

we obtain ^Ks=0.026 eV, almost identical to the thermal 

9 

•energy of 300 K. By putting TiKs^^T and -g- — g- 1 , 

Eq. (A. 2. 6) becomes 


R - 


k T 
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Tim 
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(A. 2. 7) 
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The mean free path for optical photon emission 
is then given by 



(A.2.S) 
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A. 2.3 


QUANTUM ENERGY OF OPTICAL PHONON 


Based on experimental Raman Spectroscopy re- 
sults (Latour, 1977) of PVDF-PTFE we have determined 
the main vibrational frequency of the normal modes in 
Teflon to be %0.11 eV. There are three vibrational 
modes with energies close to this value: 
antisymmetric vibrations' of CF 2 and CC, rocking mode 
of CF 2 . 


We also use another approach to check 
the reliability of the above value for optical phonon 
energy. From the polymer structure of Teflon, we as- 
sume the bond strength ’of CF is in between that of 
dimolecules HC1 and CC1 . We are then able to locate 

the symmetrical vibration frequency of CF 0 by the re- 

2 ^ 

quirement that p^w =constanx for a given spring constant, 

where p, is the reduced mass and to is the vibration fre- 
A 2 ? 

quency . p.m =.141 and .099 m (eV) respectively for HC& 

A P 

and CCA (Handbook of Chemistry and Physics, 1969). For 
symmetric vibration of CF 2 , the mass of carbon can be set as 
infinite. The reduced mass of the system is then 19 m 

P 

which in turn leads to 0.086 eV and 0.073 eV for CFg 
vibration. The use of 0.11 eV therefore seems reasonable. 
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A. 2.4 


IONIZATION THRESHOLD AND IONIZATION CROSS 
SECTION 


The optical properties of TEFZEL (1973) indi- 
cates an energy gap between valance and conduction bands 
of about 7 eV. The sharp rise of the absorbence for 
smaller photon wavelengths corresponds to the fact that 

the edge of valance band is reached. For teflon, we 
adopt 6.5 eV as the energy gap. 

We also roughly check, using the idea of bond 
strength, the order of magnitude of this 6.5 eV energy 
gap. From experimental data (Handbook of Chemistry and 
Physics, 1969), ‘the bond strengths of HF and CF are re- 
spectively 153 and 107 Kcal/mole. From the ionization 
potential of HF, 17.7 eV, we infer that the ionization 
potential for CF is roughly 12 eV which is about a 
factor of two larger than the value we' used. 

The ionization cross section for electrons 
with energies above the threshold is not known. For- 
tunately, theories (Wolff, 1954; Baraff, 1962; DiStefano 
and Shatzkes, 1976) indicate this quantity is not impor- 
tant in the calculation of charge multiplications. We 
— 15 2 

have chosen 10 cm , a large enough number, as the 
cross section for ionization. 
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A. 2. 5 


SUMMARY OF PROPERTIES OF FEP TEFLON 


The following summary of the physical proper- 
ties of FEP Teflon have been taken in part from Dupont 
Technical Bulletin T-1C. 


Density 

Dielectric Constant 
Volume .Resistivity 


2 . 15 gm/cm^ 
2.0 

10"^ ohm-m 
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Appendix 3 


DESCRIPTION OF XSCPRP 


I. PURPOSE 

XSCPRP is a computer code to prepare the cross 
section set for SEMC. 

II. 1 START-UP REQUIREMENT - INPUTS 

The following inputs are needed for an XSCPRP 

run: 

CARD TYPE 1 C12A6) COMMENT (I), 1=1, 12 
A comment card defining the cross section set, 

e.g., 

N2 CROSS SECTION SET 
CARD TYPE 2 C 13) 

NDISC — ntimber of discrete excitation levels. 
CARD TYPE 3 CI3*F5.1) 

NION — number of energy points in ionization 
cross section set. 

THI ionization threshold (eV). 

CARD TYPE 4 C3E9.2) 


EION (I) rl=1 /NION — energy points for ionization 

cross section set. 



CARD TYPE 5 (8E9.2) 

SIGIONCl)rI = 1rNION - cross section- for 

ionization at the 
corresponding energies 
EION(I) (cm 2 ).- 

For each excitation level, K=1,NDISC, read 
CARD TYPE 6 (I3,F5.1) 

KDISC(K) — number of energy points in the 

excitation level. 

' THDISCCK) — threshold for the K excitation 
level (eV). 

CARD TYPE 7 C8E9..2) 

EDISC ( I /K) *1=1 ^KDISC (K) —energy points for 

the excitation 
level (eV) . 

CARD TYPE 8 C8E9.2) 

SIGD (I.*K) r 1 = 1 *KDI SC (K) — cross section for the 

excitation level, 
at (cm?) 

energies EDISC(I,IC). 

CARD TYPE, 9 (13) 

NELAS — member of energy points in 

the elastic cross section set. 

CARD TYPE 10 (3E9.21 

EEL (I)' '*1=1 *NELAS —energy points for the elastic 

cross section set (eV). 

CARD TYPE 11 (8E9 . 6) 

SIGEL(I) *1=1 *NELAS — cross section for elastic 

scattering at each of the 
energies EEL(I) (cm 2 ). 

II. 2 START-UP REQUIREMENTS r- DISK FILES READ 

No disk files are read. 
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Ill . 1 


CODE OUTPUTS 


PRINTED OUTPUT 


The code prints out the input data as read. The 
number of points in the energy grid for the total cross 
section and the master energy grid for the total cross 
section are printed. Finally, a tabulation of the total 
cross section, probability of ionization, probability of 
elastic scattering, probability of excitation, and the average 
energy loss per excitation event are printed at each energy 
point in the master energy grid. 

III. 2 CODE OUTPUTS - DISK FILES WRITTEN 


The following information is written to file XSEC 
in binary format: 


RECORD 1 HEADERCI) 
HEADER (1 ) - 

HEADER (2) - 
HEADER (3) - 
NEMG - 


✓1=1 /3^NEMG 

material for which the cross 
section set is written. 

date of preparation . 

time of preparation. 

number of points in the 
master grid. 


RECORD 2 EM (I) ✓1=1 ✓NEMG 

master energy grid (eV). 

RECORD 3 SIGTOT(I) ,1=1 ✓NEMG 

2 

total cross section (cm ). 

RECORD 4 PELECI) ,1=1 ✓NEMG 

probability of elastic scattering. 
RECORD 5 PEX Cl) ,1=1 ✓NEMG 

probability of excitation. 

RECORD 6 AELECI) ,1=1 ✓NEMG 

average loss per excitation event (eV). 
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V. 


PROGRAM LOGIC 


The logic flow in XSCPRP proceeds as follows : 


1) Input data read. 


2) A master energy grid is defined to include 
all energies in all cross section sets, 
elastic, excitation, and ionization. 


3) The total cross section at each point in 
the grid is calculated. 


4) The probability of each event type, elastic, 
excitation, and ionization is determined. 


5) The above results are printed and written to 
disk file XSEC. 
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Appendix 4 


DESCRIPTION OF SEMC 


I. PURPOSE 


SEMC is a single scatter electron Monte Carlo 
code which calculates the drift velocity, ionization length, 
the diffusion coefficient parallel to the electric field, and 
the diffusion coefficient normal to the electric field. The 
spatial distributions of the ionization sites are also, 
scored. 

II. 1 START-UP REQUIREMENTS - INPUTS 


The inputs for the SEMC code are in namelist 
format. Two namelists SIV and $0F are needed for each data 
set. The run is terminated with the following card: 

SIV HXHIST=-1$ 

The namelist format was chosen for ease of use, as well as 
for the ability to define default values for the variables 
in the namelist. The variables in namelist SIV are: 
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VARIABLE 

DEFINITION 

UNITS 

DEFAULT VALUE 

EZ 

' 

electric field 

volts/m 

10 7 

R VMD 

density of atomic 
bands 

/ 3 
/m 

5.14 x 10 27 - 

ESTRT 

■ 

initial energy of 
electron 

eV 

0.02 

MXHIST 

number of electron 
trajectories to be 
followed 


200 

•DR BN 

bin width for 
scoring radial 
location of ioni- 
zation sites 

m 

1.0 

DZBN 

bin width for 
scoring axial 
location of ioni- 
zation sites 

m 

1.0 

DTBN 

bin width for 
scoring temporal 
location of 
ionization site 

sec 

1.0 

TSCORE 

time to score 
the diffusion 
coefficients 

sec 

1.0 









In order to get the values for variables DRBN, 

DZBN, and DTBN, 'it is usually necessary to make a preliminary 
run with a small number of electron trajectories ^ 200 . The 
information printed on the page of output labeled DIAGNOSTICS 
is used t.o determine the bin sizes. The values are chosen 
as : 


DRBN > 0.1 * MAXIMUM RHO 
DZBN >0.1 * MAXIMUM Z 
DTBN > 0.1 * MAXIMUM TIME 

The value of TSCORE is determined from the value of ioniza- 
tion time printed on the first page of the output labeled 
RESULTS. 


. 2 5* 10 N I ZAT I ON TIM E<TS COR E< . 5*I0N IZATI0N TIME 


The variables in namelist .$OF .are: 


VARIABLE 

DEFINITION 

DEFAULT 

PRNTXS 

logical variable to request, 

. FALSE. 


print of cross section set. 


PRNTMD 

logical variable to request, 
print of the spatial and 
.temporal distributions of 
ionization sites. 

. FALSE. 


II. 2 START-UP REQUIREMENT - DISK FILES READ 

SEMC reads the XSEC file prepared by XSCPRP. 
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III.1 


CODE OUTPUTS 


PRINTED OUTPUT 


The code prints out the input data as read. The 
cross section identifiers as read from file XSEC are 
prinxed next. If PRNTXS is set to -TRUE, in namelist $OF, 
the cross section is printed. On the next page, labeled 
DIAGNOSTICS, only those values previously described are of 
general interest. On the following page, the drift velocity, 
ionization length, ionization time, mobility, the diffusion 
coefficient parallel to the electric field, the diffusion 
coefficient normal to the electric field, and the ratios of 
the diffusion coefficients to the mobility are printed. 
Finally, the spatial and temporal distributions of ioniza- 
tion sites is printed if PRNTMD is set to .TRUE, in $OF . 

III. 2 CODE OUTPUTS - DISK FILES 

No output disk files are written. 

IV. PROGRAM LOGIC 


The logic flow in SEMC is as follows: 

1) Input data read and printed. 

2) Cross section set read and printed 
(if requested). 

3) Random walk variables -initiated. 

4) Electron is started at energy ESTRT. 

5) The distance and time to next collision 
site is determined. 
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6) The electron is moved to the site of the 
next collision and the electron energy 
updated. 


7) The collision type is selected. 

8) If collision is ionization, the ionization 
parameters are scored and a new electron 
begins at 4) if number of trajectories is 
less than MXHIST. If number of histories 

is equal to MXHIST, the results are printed. 


9) If collision is elastic or excitation and 

the number of elastic scatters is less than 
2000, then the next step of the random walk 
is taken at 5) . 


SAMPLE PROBLEM INPUT PARAMETERS 


The sample problem was for one thousand elec- 
tron histories incident on Teflon with an applied 

9 

electric field of 10 volts/m. The inputs for this 
problem are 


SIV 

EZ = 1.0E+9, 

RVMD = 5.14E+28, 
MXHIST = 1000, 

DRBN = 3.0E-9, 
DZBN = 3 . 0E-9 , 
DTBN = 2.0E-14, 
TSCORE = 7.5E-14 
#END 

$OF 

PRNTI.ID = .T. 

$EITD 

$IV 

MXHIST = -1 
$END 
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VI. SAMPLE PROBLEM OUTPUT'S 

Sample output from SELIC for the inputs listed 
in Section V are given below. 
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S E M C 

SI NGLE -SCA HER F LFC TRO N MONT E C AR.L Q CO DC 
(032279 ,h37q2 ) 

INPUT S 








S E M c 

SINGLE-SCATTER ELECTRON HOME CARLO CODE 


TOTAL OF 1DDD HISTORIES 

TOTAL OF 995 IONIZATIONS 

TOTAL OF 0 RUNAWAY ELECTRONS 

TOTAL OF 0 HISTORIES EXCEEDING 2000 ELASTIC SCATTERS 

AVERA6E NUMBER OF ELASTIC SCATTERS PER HISTORY 108.059 

AVERAGE NUMBER OF EXCITATIONS 23.68 1 

MINIMUM Z = 5.549-Q09 

MAXIMUM Z - 2.337-008 

MINIMUM RHO = 5.455-Q11 

MAXIMUM RHO - 2.858-008 

MINIMUM TIME = 7. 600-Q1S 

MAXIMUM TIME - 4.960-013 

XBAR - 1.422-010 

X2BAR z 2.34?-ni7 

YBAR - -1. 006-010 

Y2BAR = 2. 630-017 

ZB A R = 8.750-009 

Z2BAR = 8.196-017 ' 

RBAR = 5.810-009 

R2BAR ~ 4.972-017 

RZBAR = 5.629-017 ’ 


(032279 ,113702 ) 

DIAGNOST ICS 






ORIFT VELOCITY 


IONIZATION LEN6TH 


IONIZATION TIME 


MOBILITY 


DIFFUSION COEFF (PARALLEL) 

DIFFUSION CQEFF (NORMAL ) 

0 (PARALLEL IMMOBILITY 


D 


{N0RMAL)/M0BILITY 


S E M C 

SI NGLE-SCATTER ELECTRON MONTE CARLO C ODE 
(032279 ,113702 > 

R E SULT S 


1.63 + L05 (M/SEC) 


1,26-008 <M) 


1.15-01 3 

(SEC1 



1 . 6 3 - CO 4 

( M* * 2 ) / (VOLT-SEC » 




(M**2/SEC> 



1 . 5 6 - CO 9 

<M**2/SEC> 



2 . 4 H - CO 1 

(EV > 



9 . 5 5 - CO 1 

(E V ) 






MARGINAL DISTRIBUTION IN 2 


2 PROB 


( M > 


] . 50-009 

.000 

<<.50-009 

1 .005-002 

7.50-009 

6.372-001 

1.05-008 

2.673-001 

1.35-008 

6.332-002 

1,65-008 

1.407-002 

1.95-008 

7.035-003 

2.25-008 

1.005-003 

2.55-008 

.000 


2.85-QQ8 .000 



S E M C 

SINGLE-SCATTER ELECTRON HONTE CARLO CODE 
<032279 ,113702 1 

R E S U L T S 



MARGINAL DISTRIBUTION IN RHQ 

RHO PROB 

(M) 

6.00- 009 2. 643-001 

9.00- 009 

1 . 20-008 
1.50-008 
I. 80-008 
2.10-008 
2.40-008 
2.70-008 ' 5.025-003 

3.00- 008 1.005-003 

3.30-008 1. 005-003 



3.377-001 
2 . 101-001 
1. 156-Q01 
4.523-002 
1.307-002 
7.035-003 






MARGINAL DISTRIBUTION I N TIME 


TIME 

PROB 

(SEC) 


3.00-014 

7.638-002 

5.00-014 

2 . Q80-Q01 

7.00-014 

2.050-001 

9.00-014 

1.417-001 

1.10-013 

9.950-002 

1.30-013 

6.935-002 

1.50-013 

6.231-0D2 

1.70-013 

3.618-002 

1.90-013 

3. 11 6-002 

2.10-013 

7. 035-002 


STOP 


IQ 


S E M C 

SINGLE- SC ATTER ELECTRON MONTE CAR L 0 CO DE 

(0 32279 , 1 1 37 02 ) 

RESULTS 
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DESCRIPTION OF CASCAD 


X. PURPOSE 

CASCAD is a Monte Carlo code which calculates 
the electron avalanche process. A single electron is 
initiated; the location of the next ionization site is 
randomly sampled from the distribution of ionization sites 
calculated by SEMC. The process is repeated for all of the 
electrons in each generation. After a predetermined number 
of generations, the charge density, potential, and electric 
field are calculated. If the maximum electric field' exceeds 
a pre-defined fraction of the electric field, the problem 
is terminated and the positive charge density and the 
electron density are separately written to disk. 

II. START-UP REQUIREMENTS - INPUTS 

The inputs for CASCAD are in namelist format . 

Three namelists, SCNTRL, $GEOM, and SDISTR, are needed for 
each run. The variables in namelist $CNTRL are: 
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VARIABLE 

DEFINITION 

DEFAULT 

EPREL 

relative dielectric 
permittivity . 

1.0 

ERRMAX 

relative error for the 
iterative Poisson solution. 

0.001 

ITERflAX 

maximum number of iterations 
for the iterative Poisson 
solution . 

500 

EFAC 

fraction of the applied 
electric field for run 
termination . 

0.5 

NGEN 

maximum number of generations . 

13 

MCHEK 

generation to begin field 
calculation . 

S 


The variables in namelist SGEOM are: 


VARIABLE 

DEFINITION 


DEFAULT 

If-IAX 

number of grid points 
in the radial direction. 


' 

20 

J MAX 

number of grid points 
in the axial direction. 


21 

L 

axial length. 

meters 

— 

A 

maximum radius. 

meters 

— 
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The variables in namelist $DISTR are: 



DEFINITION 


' EAPP 

value of the applied 
electric field for 
this data set. 

volts /m 

N1 

number of points in 
the radial distribu- 
tion of ionization 
sites. 


N2 

number of points in 
the axial distribu- 
tion of ionization 
sites. 


RMAX 

maximum radial 
displacement for the 
distribution of ioni- 
zation sites. 

meters 

ZMAX 

maximum axial 
displacement for the 
distribution of ioni- 
zation sites. 

meters 

PR Cl ) 

radial distribution 
of ionization sites. 


PZ Cl ) 

axial distribution 
of ionization sites. 



II. 2 START-UP REQUIREMENTS - DISK FILES READ 


No disk files are read. 
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III. 1 


CODE OUTPUTS 


PRINTED OUTPUT 


The code prints out the input data as read. 

Then the cumulative distribution of ionization sites are 
printed. The maximum values of the radial field, axial 
field, and maximum magnitude are printed for all genera- 
tions greater than NCHEK. 

When the maximum number of generations is 

reached or the maximum electric field exceeds EFRAC*EAPP, 

the spatial distribution of charge density RHO (in units 
3 

of coul/m ) is printed. Next, the electrostatic potential 
PHI (in units of volts) is printed. Then, the radial com- 
ponent of the electric field ER and the axial component of 
the electric field EZ (in units of volxs/m) are printed. 
Finally, the number of positive and negative charges in 
each cell are printed. 

III. 2 CODE OUTPUTS - DISK FILES WRITTEN 


The following three binary records are written 
to file COUT: - 


RECORD 1 IRSvIZS/'DRSrDZS 


IRS — number of cells in the radial direction. 
IZS — number of cells in the axial direction. 
DRS — radial cell size in meters. 

DZS — axial cell size in meters. 


RECORD 2 
RHON 


( (RHON (J ✓ I ) ✓ J =1 ✓ I RS ) ✓ 1 = 1 ✓IZS) 

— negative charge density in each of the 
cells (coul/m^). 
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RECORD 3 


( (RHOPCJ/I) ,J=1 ,IRS) ,1 = 1 ,IZS) 


RHOP — positive charge density in each of the 
cells (coul/m^). 


IV. PROGRAM LOGIC 


The logic flow in CASCAD is as follows: 


1) Input variables are read and printed. 


2) - Random walk variables are initialized. 


3) Each electron in the present generation is 
moved to the next ionization site. 


4) If the number of generations is less than 
NCHEK, go to 3). 


5) Calculate the charge density, potential, and 
electric field. If the fields are less than 
EFRAC*EAPP, go to 3). 


6) Print the charge density, potential, and 
electric field. 


7) Write the positive and negative charge 
densities separately to disk. 


V. SAMPLE PROBLEM INPUT PARAMETERS 


The sample problem was for thirteen electron gen- 
erations in Teflon for an applied electric field of 
9 

10 volts/m. The inputs for this problem are ■ 
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$CNTEL 

EPREL 

_ 

2.0, 



EFRAC 

= 

•75, 



NGEN 

§END 


13 



$GEOM 

SEND 





SDISTR 

EAPP 


1 . OE-s-9 , 



N1 

= 

10, 



N2 

= 

10, 



RMAX 

= 

3. OE-8, 



ZMAX 

= 

3.0E-8, 



PR( 1) 

- 

.264, .338, .210 

/ .116, .0452, 

.0131 



.00704, .00525, 

.001, .001, 


PZ(1) 

= 

0.0, ,00i; .637, 

.267, .0633, 

.01407 



.00703, .001, 0. 

0, 0.6 


SEND 






VI. SAMPLE PROBLEM OUTPUTS 

Sample out from CASCAD for the inputs listed, 
in Section V are given below. 
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> 

Ol 

I 


•t. . 1 


1 ' L 

- 

, ft >•! i , 

t AX 


. It -OC, 

ITtKf-X 

- 

500, 

If^AC 

= 

, 75L * 00 

flGU, 

= 

u, 

JJC.HL* 

= 

d, 

SE(nO 



Mil. 1 



\ *•■- 1 

= 

RO, 

J.MAX 

= 

21, 

L 

= 

■ oC-G o , 

I. 

_ 

. ibc-Oo. 


SELD 



RECCRO .'<0, 11 P2(l)«0,,.00j,.o37».2b7»,063i#.01fll».«070i# ,001*0. ,0. 

12iBSi7890121<l5t7e‘)ul23ijSo/b9UU3^bo7 890l2Jillso/e-)0.........' i 2S«io7bb‘0l 23 Abe /KOI) 1 2i«b6 /«<,(; l RlAlb/Bvl, 

CC^A RlitlM, AT £A0 OK RECORD - CCRMA A&SLRLD 

ERROR fiORE’c-R E<3 DETECTED bV IflCLMs Af ADDKESb 000207 

CALLED f ROM NA«lf.s A1 ADDRESS 0 0 0 1 41. 

called from cascad a? line 65 

.llrl’i, 

• 10, 

= 10, 

> .it-, I, 

• ,5i-07, • 

.2o<l£ + 00, .Jiei+00, .21EKOO, .IlcL+OO, .«i2E-01, ,131L*01, .70 hE- 0?, .b2St-0 2 , .lt-02, .1L-02, 0.0, 0.0, 

0*0/ 0*0/ O a 0/ 0,0/ 0*0/ 0*0/ 0*0/ 0,0/ 


0.0, .11-02/ .6i7L*0O, ,2c)7fcrCiO/ ,i>4iE~Vi, ,7U3fc-0<> 

0*0/ 0 , 1 / / 0,0/ 0,0/ 0.0# (1.0/ 


1 f, - 0 2 / i'*0/ 1,0/ 0*0/ 0,0/ I' . U / U*U/ 



-sv 


i'l V !•« 1 1 of IUI»i 


I IJf- / 


l 


1 

0. 

0. 

2 

b,0Gt,“<)9 

0. 

J 

1 .Ovt-Gfe 

1.00E-03 

if 

1 , DOE "00 

o . A Be -01 

b 

2.0AF.-0B 

9 , Obb ”0 1 

b 

2 . bOt -OB 

9.or,e-0i 

7 

3.001 -OB 

9,b,>f.-oi 

e 

3. bit -Lfl 

9 . B 9 L - 0 1 

9 

ft .OCb-Ufl 

9 . 9 0 b - 0 1 

10 

4. b0b-06 

9.90F-01 

1 1 

S.00L-0B 

1.00F+00 

J 

RnC 

FOFK 

1 

0. 

0. 

2 

3.00E-09 

2.B4E-01 

A 

o.rvt-09 

o.Ci'E-lf I 

R 

9 . Dor -09 

.i. 12E-01 

c 

1 ,2ut-C& 

9, 2l*E -0 1 

0 

1 .bOt-US 

9.73I-0I 

7 

1 . fiOE -03 

9.8BE-0I 

3 

2. lOE-03 

9.93E-01 

9 

2.90b -OB 

9. 992-01 

10 

2. /Cb-08 

l.OoE+OO 

1 1 

3.00b-CB 

1 . OGE+UO 


Gt-'if HiTlUi. 0 

MAXIM U>' nADIAL Fli.ll) i.« , W?t»07 

MAXIMUM AXIAL f II LD /f.OOBt + O/ 

MA < UX'M MAI1.NI TUDF 0.1‘»7t + 07 

OtM-PAlIul. 9 

MAXIMUM RADIAL MfcLO H,e03E.*n7 

MAXIMUM AXIAL rltLO Q.cJottU/ 

‘•jAXIMLF KA6NIIUUL l.H,?t + 0ti. 

CEM'RAJIOU 10 

Maximum radial FIELD 1.071EF0B 

MAXIMUM AXIAL I'lsii) 1 , 7 A i £ t U c 

MAXIMUM MAUN1 TU.'fc l.90!L<Do 

DLobRA | lUU 11 

MAXIMUM RADIAL FIELD 2.349E+0B 

maximum axial FIELD 3. lAOttOB 

".AxJMlM MAONI JUDE i,fc70b+OU 

CEf.f -fATIOIx 12 

MAXIMUM radial FILLD 3 . I r 3 1 V o m 

MAXIMUM AXIAL r J 1 lO o . i 1 - b t G il 

MAXIMUM •■■■AtM iUOr. t>.«94b-t&B 

obMhhAI K-ii 1J 

. f- AxiMUM RADIAL fit ID /.,:HiE-t01i 

■MAXIMUM AxIAL FIELD I.II'.EtOO 

MAXIMUM MAGNITUDE 1.22<>Lt09 




Ol 

I 

CO 



«2 

f* 

5<s 




, 


i 

t , 


i 

l\ 1 1 

, 

/. 





0. 

-7 

.OOc-b.J 

o.DOl -Ob 

0 

, Ul't-Oft 

l.iJOt-C/ 

1 , ") <1 I. - 0 t 

1 ,C oi ' *• / . *• V*. 

-V 7 

f 1 . L » J i 

I 

t 

!•( 1 7 

3, fit.- 0>3 

0. 

0. 


0. 

B, 

3 0 4 E 1 0 4 

e.no/t+o'i 

7.o'.‘9) lub 

“c , Sb Ji Of» 

-b 4 4bl f *ub 


i.ist-oa 

0. 

“7. 

bb3E+03 

4 , 30 OF + 04 

7 , 

■ / Li +00 

2.391c tOb 

h. 1 7 /ttfi 

\ , V/ •- S: ▼ i? b c » ] 

T«*<J 

- * . . » 7 , 1 f y 

3 

1 .12E -0# 

0. 

0. 


0.. 

n . 

023- +04 

2 . JS3t f 0‘i 

7.ilc'7tiOS 

1 , *; t* * * w > i » l w i - c. * v b * I , s i • : 1 

' *'%’ 

- 3 k :»J- 

— 

2.)>‘-' 

c . 

0. 


0. 

b . 

im+03 

o . Vb ot,+ 04 

i .milt t Ol 

l.l/cM'.o i'f /: >V0 M.I'U 

> 1+ *1 

" > * \f Z 0 • . *J 

b 

i , ll-ll -!,<) 

•0. 

0, 


0. 

0. 


c. / 10L *04 

1 , -*300 + 03 

4 &ui t - 1 „ i/ •*,<?!. - «;b -b./fts 

t0‘) 

“ l • * > 1 ♦ »/ 5 

b 

K23J -08 

0 . 

0. 


0. 

0. 


.-.60UE + 03 

0. t 7 if +00 

1, /oll-ibb -o . i-l+.t * to -2.fc‘*bl+0b 


7 

5. 01.1- -Cd 

0. 

0. 


0. 

0. 


0. 

2.207E+03 

2,b3hM04 -b.blOt+Oii -K2’j3t + 0b 

-a. 31 vttoa 

& 

b.77c-08 

0. 

0. 


0. 

0. 


0. 

0. 

4. 7821 +03 -Kloof +0b — 1 ,bl 7 ). + 0 4 

-UViiLfOi 

9 

o.s«r -ob 

0. 

0. 


0. 

0. 


0. 

0. 

0. - 1 . i-bt-t +03 -2,b_>lc 

*yS 

~1 , obt'&tOi 

[0 

7.312-05 

0. 

0. 


0. 

0. 


0. 

0. 

u. -7 , bboE ) 72 -7.3307 

* 

U, 

1 1 

.l.Otf "-Oil 

0. 

0. 


0. 

0. 


y. 

0. 

0 , » , 0 , 


u k 

12 

<>. bBr -Oil 

0, 

0. 


0. 

0. 


0, 

0. 

0 . 0 . 1- . 


t; fc 

13 

r / . 1. 1? L — Ob 

o . 

0, 


,o. 

i>. 


0. 

0. 

0 . 0 . 0 . 


0, 

14 

1 t fib£-u7 

0. 

0. 


0. 

o, 


0. 

0. 

o . o . o , 


u. 

13 

l,12t-u/ 

0. 

0. 


0. 

0. 


o. 

0. 

0. 0. 0. 


o. 

1 O 

1.1-71-07 

0. 

1). 


0. 

0. 


0. 

0. 

0 , 0 , 0 . 


o. 

17 

1 ,2/E-O? 

0. 

0. 


0. 

0. 


0. 

0. 

0, 0. 0. 


0. 

lb 

l.iSE-O/ 

0. 

0. 


0. 

0. 


0, 

0. 

o . • o . o • 


0, 

i V 

1.42E-07 

0. 

0. 


0. 

0 . 


0. 

0. 

0 . 0 . 0 . 



20 

1 .bOt-07 

0. 

0. 


0. 

0. 


0. 

0. 

0 , 0 . U . 


0* 




2 

7 


L 

L 


/ 

L 


/ 

/ 



3.O0E-07 

3. 30E-07 

' 3. 

bOE-07 

ip70i“07 

^,20f-07 

O.bOi -07 

o (>t -i)7 

b. 10J.-0/ 

b. 7 

b * 7 vi / 

] 

i.(l) 












i 

3. obt -09 

-b .730C+04 

-2.2691. +04 

(1. 


^ • 

0* 

0. 

ft 1 


t . 

«l p 


1 , 1 5i--0o 

- 5 . 736ct04 

-9 .S7.3-- + 03 

0. 


0 , 

0. 

0. 

y * 

J. 

b . 

o. 

> 

1 ,‘.2f “OlJ 

-2,c9bE-04 

0. 

0, 


\J * 


0. 

0> 

0 4 

0. 

li, 

y 

2,1 9 1 -Oil 

"2,«'J*lLt04 

0. 

0, 


o » 

V , 

0, 

li. 

0 . 

0. 

0 • ’ 

V, 

3 ,4ot -OB 

« l . 59HE+04 

0. 

0. 


0 . 


0. 

0 1 

0, 

1 , 

y , 

o 

4 . 23t -CO 

-b.^lbC-03 

0. 

0. 


0. 

o. 

0. 

IJ , 

0 , 

0. 

*j . 

7 

b . (.01 -OB 

-l.103t.0i 

0. 

0. 


y * 

o. 

0. 

0. 

0. 

V • 

0 . 

o 

b, 7 7t -Ob 

0. 

0. 

0. 


c. 


0. 

0. 

a , 

•/ , 

y ♦ 

9 

6 . b 4 r -Ob 

-H.436L+02 

0. 

0. 


0* 

0. 

0. 

0. 

o t 



10 

7 . 3 1 r -3b 

0. 

0. 

0. 


0, 


0. 

0 * 

V. 

0 . 

y i 

i t 

‘n. f*‘. --J b 

o . 

0. 

0 • 


Ok 

o. 

3. 

Vi • 

y+* 

li 4 

u. % 

li 

t'.oii^Oe 

0. 

0. 

0 . 


li t 

y . 

0, 

u < 

c , 

, 

i>. 

n 

<( . o2“. - 06 

0. 

0. 

0. 


u , 

<i* 

0. 

‘ 0. 

0* 

0* 

■ii. 

: i. 

1 .(.at -07 

o. 

0, 

u. 


0. 

o, 

0. 

V-. 

0 * 

a* 

0. 

i b 

1 . 12C.-07 

0. 

0. 

0. 


y. 

0* 

0. 

0, 

y* 

Vi • 

U 4 

\ o 

1,191-0 7 

0. 

0. 

0. 


. u. 


0. 

y. 

0* 

0. 

0, 

*7 

1 .2/£-07 

0. 

0. . 

0. 


0, 

»)* 

0. 

l>. 

s> * 

V ♦ 

0. 

U 

l.Ibt.-07 

0. 

0. 

0. 


0* 

o» 

0. 


0. 

li . 

V* 

Pi 

1 . i2i-0 / 

0, 

0. 

0. 


u. 

Ol 

0. 

0, 

y. 


V, 1 « 

20 

1 .bOc-97 

0 „ 

0. 

0. 


0. 

y f 

0 . 

0, 

0, 

1? , 

0 , 



A5-10 


c. ,0<JL -07 


t 


* . 

1 

1.15'-. ' 

0. 

1 

1 .‘.YL <1 

0. 


2, O'*. “00 

u. 


A, . ■« ft 1 "CD 

a. 


•1 . ^ 3fc— y fa 

o. 

7 

S.(M.L-Ot) 

v • 

e 

b./’l-OC 

0. 

<> 

t. ,‘j <h -On 

c , 

i V 

/.AIL. .J « 

0. 

. i 

i! . -1/11 

0. 

12 

B„ eb£ -oo 

0. 

12 

<?.62fc-0U 

0. 

14 

1 .0ilt-07 

0. 

15 

1 ,l2i.-07 

0. 

10 

1.19L-07 

0. 

17 

1.2 71-0/ 

0. 

1 *> 

1 . J'jl-07 

0. 

1« 

1.421 -07 

0. 

20 

1 ,50t-07 

0. 



A5-11 


i : ; l 


I 

Ml) 

0 • 

^ . (J m i " J ■ 

6 4 o j ; »vo 

•/ » U 0 1 * ** 0 

1 / 

1 ,‘> 0 i - 0 / 

: , t • .* / 

<r 4 > “ j f 

4.1 / 

1 

i 

3.060-O9 

1 ,5306*00 

2 , 6 9 9 1: + 0 0 

0.1536*00 

/ , <)39k* OO 

1 .2901 +01 

/.Jill 101 

> . i 7 i t t *J 1 

-1 ,/Olf +00 

-.3 . 1 o"> 6 + 0 1 

-i ,uiLi +ii : 

1 

1 . t3f -Ob 

1 ,V6o6+06 

2.6996*00 

0.1636*00 

/. 0396+00 

1 .2006*01 

2 , 3191 *01 

3.1716*01 

- 1 . /U 2 (.*(J 8 

• 5 * ] tO 1 

-1.9s-. *.!■ 

5 

1 . 020-015 

1.9906*00 

2 • 6 jot *00 

3.9/66*00 

8 . 7 ’ 1 6 * 0 0 

1 , 2 l 2 r +01 

2. 1 1 /! +01 

7.6?3t *01 

- i f * t* b c. * U l> 

-8 . (‘526*0 1 

- 1 • ftv 1 } tl'l 

n 

2,091 -0 fi 

1 .9176*00 

2.5596*01) 

3. 70OF * (10 

0.17 (11. + 0M 

l.M/ 4k *01 

1 .ili/Of -01 

2.2 71-1 *01 

- : * ? / *uu 

- i , t t »> i 

“*•<»'? f *• i 

< 

i * c f - 1> c 

1 4 d*'/r + <jo 

2.970t»'jw 

3 . 5 0 0 6 r 0 0 

* 0 (» 

9.32 If ,0 0 

1 « 4 / o', y '• 1 

1,/,'H , V, 

- 1 . ^ \ i/ r ' / *7 

- 1 , *» S 

“2 ..lO. • . i 

0 

‘..iTJv'Ob 

1 * t b U \ * if o 

li 4 3 o “ I. * 0 0 

, 3 31'f *00 

b * *00 

£■ ,0 0i ICO 

1 . ! 7‘‘; -0 1 

1 , i ■> *» r » »* 1 

-l+i'i/U * u u 

- 1 4«nn{.i If 1 

-..ii.'.i i 

i 

b • \> 0 - “On 

1 . 6oi!l'*00 

2.2566 +«0 

3. 1 006*0(1 

9,60 1i_t00 

f> . fcobfr 1 0 U 

9.0106 *00 

9. lOol, too 

“ 1 . btot M)U 

- 1 , l !. r L**»l 

-•* , j‘: -t - - v 

i 

3 . 7 / 1 -Of) 

1 . /6I6*0U 

2.1366*00 

2.83 7F *06 

0.0916+00 

6 , t* t- t ■* 0 0 

7. 5501 *00 

0,07 7c *00 

-1 .2301. *00 

•f* , /b'l» t Mil 

“ / . • r'Vi * vi 

9 

o „ b 4 i - u 3 

1 .7096*00 

2.0086+00 

2.6136*00 

3.0226*00 

0 ,9931 ♦ 00 

*U(J 

0.9/O6*OU 

-9.50/6-01 

"C « ♦ 'JU 

*■ U 4 <V V » 0 

10 

7,311-06 

1.6696*00 

1 . b /o6+ 00 

2 . 3 7 1 6 1 0 u 

3 . 1006+00 

0 .23,-6*00 

0.9(17f +00 

3. 7511 *00 

-7 , 68 3b -W l 

-'J, /cut* UU 


1 1 

8.0(16-06 

1 ,59t;£*0O 

1 ,79 lE+Oo 

2.1326+00 

2. /53C+00 

3.67/f + OO 


2.1)036*00 

-5, .0/91 -01 

-3, (,9,.( * ou 

-a, l * | r tuu 

12 

A , o5L -0 6 

1 .i3>r*oc 

1 . oO 26*00 

1 .aObt+OO 

2.0106*00 

3.0036-00 

i * U ■* «*u 

2 . 2u 36 - i0 

-0,1) )(•( -0 1 

tsl J 

i * i Ci’*! ▼ ' v 

\ i 

-.621 -Oe 

1 . 9 the. *00 

I .JfcOf *00 

1,6/06*00 

2,0/26*00 

2 , 6021 .* 00 

2.5/ 1 6 *-0 

1 „ / 0 i : i Oil 

- ).l 6l 1 -!• 1 

-i.f'O! • W 

-r * ^ 

1 

1 .Oil -0 / 

1.-15: too 

! . Sl ! ,t ‘00 

1 , 0 9 / L * u o 

1.7606+00 

£.0811+0(1 

2.0001 -60 

1 . 3 1 /l * 00 

•* t * , * / 1 * \i 1 

"1 . / >• ' t T »I (J 

♦* ; 4 ' i < *i . »\\ 

1 > 

! . l2r.-0 7 

i , 3b i£* </0 

1 . 1 Ut: +0u 

1 . 2? ? C l £ T 0 0 

1 .0616*00 

1 .i/lt -00 

1 . o 3 2 r. - 0 0 

1 4 0 1 it ■* (. 0 

*■ d 4 .'loi -u i 

* ! . i'i'-t ♦uo 

•- i . / *' > < » */ , 

iu 

1 .195.-07 

1 .2936 + 00 

1.0116*00 

1 , 0 l cfc 1 00 

1 . 17 1l*00 

1 ,320( +0o 

1 , 2 0 / ( -uu 

/. (.856-01 

-1 . /)9[ -01 

- 1 . 0 03. + 00 

•• 1 . 337, * 0, 

17 

1.276-0/ 

1 .233E + 00 

6, 9 b 7 £ — 0 1 

8.0696-01 

9.0956-0 ! 

1 ,012t +00 

9.62 36-01 

5,e7 3! -01 

-1 , 2801-01 

- / 4 t> V# 0 f - V [ 

“ i • O C /' C V V t 

16 

1.356-07 

1 . 1 7 9 6 r 0 0 

o. 7006-01 

5,9996-01 

0. 8326-01 

7,21106-01 

0.71)06-01 

3.5796-01 

**9 4 i>J7r -D? 

-6,31)1 ( -01 

-/, 2/1*1-0 1 

19 

1.026-07 

1.11 /£ ♦ Oo 

0.763i:-Ol 

3.9566-01 

0,3061.-0 1 

0.00)6-01 

0 , 3 1 0 f - 0 1 

2.3126-01 


-3,^*/7r-oi 

*-tj 1 

20 

1 .506-07 

1.0626*00 

2 ,6056-0 1 

1.9626-01 

2,0981.-01 

2,2/26-01 

2 , 07SL-0 1 

1,2006-01 

-2. /09l -02 

- 1 . hi5t-« 1 

-0 1 




2 

2 

L 

i 

L 

/ 


/ 


L 

/ 




3,006-0/ 

3.306-07 

3.606-07 

3.901.-07 

0 . 20f -07 

0,506-07 

0 , Hul -07 

0 , \ i> fc - 'J / 

V 



- 0 / 

J 

*'U) 













* 

S , 6 *. 1 " 8 0 

-9.3306*00 

-<4, 8/9L-U0 

* 2 , 7 0 c 6 + Ci 0 

-1 .(01-E + OO 

-1 .u32t. , 0 o 

-7.9236-91 

- 3 .loil.-u: 

{ l »" 1 


'. / 1 "Vi 

- 3 , U'ii 

I 

d 

1 , lif -08 

-9, J356+0S 

-0..*/'/ ftJU 

- ? . / 43 9 L *00 

- 1 . 1 0 96 1 0 0 

-1 .0 031 * 00 

- / .0 Silt -91 

-5.1 (4 ("1 - 1 

- - , “ \ } u - L« 1 

- . 

V j 

- b . 1 0 J > 

- J i 

j 

! .926 -0? 

-0,5026-00 

-0, ,00^*00 

-7. / 2 o 6 + 0 C 

— 1 .1,4 0; +00 

“ 1 . 0 0 o ( * 0 0 

-O , *,33f -01 

-5.11 •( - in 

• ■ v l 

“ + , 


-*'*J • 1 t ; t 

* v 1 

H 

2 -Cr <3 

-O, 3116*0!) 

-4. 63 (.5 *00 

-2.6S5r.tOO 

- 1 , 44 031. - Ou 

- 1 .0 ! (• C + 

— 0 . (' : * 4 1 — -1 1 

-3 . 07 /( - 1 * 1 

- 0 . 04*00 - 8 1 


i' 

- >* \ftl 

* : 

b 

3 ,08( -o'- 

- 7 . 5(,2II* <n. 

- 0 , 4*8 36 + 0 0 

-2 . 64 7 )6< 4/U 

-) .34.5* *06 

-9 ,0CJ4l! - 0] 

-0 . 4,4'Ot -01 

“<* , (**4 9 t - 4. 1 

— 6.2108 ** vl 

~ ' . 

- 1 1 1 - 0 l 

"*» « 1 J 

- »• \ 

o 

0.231 -00 

-O./51C+00 

-3.5686*00 

-2 . io 76 tOO 

-1 .00 01 +00 

-9,0796-01 

- 0 . 3 9 U 0 - 0 1 

“8.7% 4-441 

“0 ,l|6fcl -0 ) 


4 1 V. - Ul 

- 1 , V V 1 f 

- 1* i 

I 

6.001 -0,1 

-5 . 5*496 vOO 

-3. s3o£+00 

-2.2 S6L + 0U 

-1 .401..+ CO 

-9.0246-01 

-6 , ] if (. c -0! 

-0.34)1-01 

-i,9.)7.--ul 

. 

; ‘ist'-sn 

«u,*/uu 

- 0 1 

3 

5. / 76-08 

-6,! 926*00 

-3,2930*00 

-2.il/0f *00 

-1,3126+00 

“8.5086-1,1 

-5.7811-01 

-0.il/0-ul 

-3, /« /t- - (i 1 

“'■5 , 

S‘J/C -ui 

- fl . / V V 7 

-w 1 

9 

8,506-00 

- (i , 0 9 /.E * 0 0 

-47.9376*00 

-1.8976+00 

-1.2)71, *00 

-7. 90.‘t -01 

-6 .021 f -ul 

”u, o6 3 6- Hi 

- 4.671,0-. 1 

“3 . 

l 

- U « 0 n V < 


JO 

/. 31c«0b 

-3 . P63C+00 

-2 ■ is206tOO 

-i . oo 

-1.11b! *00 

“7.3.4 It -01 

-3,0290-01 

-.> , 7 art 1 

- . 1 . sl2r -4, 1 

-i. 

‘•n 1 *: “U 1 


* v 1 

ti 

8.C-31 -90 

-5.2976*00 

•2. 312c * DO 

-1 .30 U + 09 

-1.0121 -Of, 

.69^ -01 


-4.6* '*4 41-51 

• >.*! 4: -.1 

, 

*f ' f ‘.-"V* t 

*' ♦ , 2 f * . 

"Oi 

! *r 

8 . «3r 11 (1 

/906 + 00 

-2,0! 16*00 

- 1 , 3o*i ^ * 0 u 

- V . 3 6 0 f - 0 1 

“o, 08 , ‘ 6 “.il 

-4 , ! 1 1 

-3, 1 If. -ul 

-<•, 6.*.“.; 


1 / J.-w, 

** 1 /;' 

« .» * 


*, * c 21 -*C9 

-fi.iioftUO 

-1 .72 76-00 

-1. l-.fiL + Oo 

- / . 5 7 66 - 0 1 

“ 6 , 5060 - 0 1 

“3. / 168.-0 1 

-2.1. -.91 -01 

If 4. 1 

“ ii « 

<1 S r ' - y 1 

•• 

- . ; 

l m 

1 , 0 o t - 0 7 

-1 ,5366*00 

“1 .«oO0tOO 

-1 , J216 4 (JO 

-6.99 : 0 - 0 ( 

“O .oof! “0 1 

- 3 1 2 0 ’ t - 0 1 

-7.65 -,t 1 

-2.7 0>+.-0 1 

P L , 

c i'yi. -u 1 

->* 7 1m 

" U ‘ 

lb 

i ,i2t -:-/ 

-1.5/06.00 

- 1 .2 1 o6* 00 

-8. 3056-01 

-6.337E-J] 

“3, 930C-O1 

-2.//U -01 

-2,13/0-11 

" 1 , 9/( 1 -V ! 



- J , - i J c. 

-vi 

J ^ 

1 . 19: -07 

- 1 , c 6 0 6 ♦ 0 0 

- 9 , 7 6 4 r - 0 1 

-4; ,9801 -1,1 

-0. 75«e-o | 

-3,2810-01 

-2.29r’6~-i 

- 1 . / 1st- \> 1 

-1 ,4,68V -.1 

“ £ • 

l»9«t "‘il 


“w 1 

* » 

1 . 2 7 1 - 0 7 

-5 , 5706 -0 1 

- 7 . 5 o 1 6 - 1, : 

»5...5c£-ul 

-3.706f.-ol 

- 4 ;.. 5706-01 

- 1 . 81 ot - 0 1 

-1.0110-01 

- 1 . 3 44,6-oI 

** 1 

c 4|(> - u 1 


"‘-‘I 

J c 

1 . >5t »t, 7 

-c. ,6566-01 

-6 , 5026-01 

“0.OOot -01 

-2.7/50-01 

-1 .9(1 It-!,) 

- 1 . 30 31 -()> 

”1 . U04 t -s 1 

- V , 9 ! 71 ■ 1 7 

-i ! 

1 

- • 0 1 i i h 

**^ i 

i9 

1 . «26-i> 7 

-- .0396-01 

*• 3, 6c3t-0 1 

-2.8006-01 

-1 ,bl7f -01 

-:.?m56-01 

-6,8; 10-02 

-6. (%«oL -u* 

-8,6 U; -v2 

*»/; , 

b 1 L *■ 1 / * 

- . - 1 



1 ,306 *0 / 

-2.1 526-01 

••1 .7366-91 

- 1 . 2 7 IE -0 1 

-o ,6oi6“02 

- 6 , 0 ‘l/L -02 

- 4| , 0 2 1 L - C 2 

-3 , 34426-02 

« J . 78 3 1 - wt 


.‘ir -0 J 

- 1 / ,«/«. ] ^ 

- Jf' 



A5-12 


Wi ■* 0 7 


1 

2 .C^L-O') 

0. 

2 

1 . l bf.-od 

-o.io7E-01 

•3 

1 . 9 2 t - 0 .i 

-t> . joOE-O 1 

‘ ! 

?.6 , ;i-o« 

-t.i2i-.-01 

5 

J,uCt-o’> 

-o , 2 ft, t-0 1 

6 

<•.201-00 

-a .7 <lci "0 1 

7 

b , o 0 1 “ <i ti 

-o.l 9‘3t-0 1 

<5 

3.7 7 L - 0 fl 

-6. llPt-01 

9 

o , 5 0 i -03 

-t. 0751-0 1 

10 

/.3lfc-0d 

-6. 00 it -01 

11 

a.oet-oa 

-5 . 930E-0 1 

12 

a.DSt'-oa 

-b .3501-01 

13 

9,622-08 

-S.7O5E-01 

it 

1 . 0 .! t - 0 / 

-5.6761-01 

15 

1 , 1 

-b .bail -o 1 

U> 

1.1V1.-0/ 

-b , Ohoc-u i 

1 > 

1.271.-07 

-b.3«9E-ai 

10 

1 ,35t -07 

-5 , 2871-01 

19 

1.021-07 

-b. 18 ;t-0l 

20 

l. bCt-07 

-5 . 0 73E-0 1 



RAGE IS 


> 

m 

i 

w 





0, j.OOL-Oi' r,.08i -Oil 

vy:. 


. .’d* -<J 7 

1 . ‘ u ! - a / 

1 

<r. l <>i - •;/ 


f' ^ / ‘*t *•» i 

t 

* 

\ 

8(1) 

3.dbl'-09 

0. I,168c*0a l,072!.*ua 

2 , 1 9 9£ * Oa 

1 . 

6 3 9 £ 1 0 3 

a ,61 btrOj 

-;.l /oLtoa 

•Mill 

•5*lc9r. -fwtj 

“ 0 * b *f C i, t vj »< 

a 

1 , 1 b£ -08 

-i.OaSEtOfl - 7 . 7331.15b -2.3!bl.tO? 

-a .80b. toy 

-1. 

Obbf H/h 

-2, bb l r. * of, 

-6 , /.till t,jb 

3 , u / u r. * u i 

^ , j *> 1 v ( ll 

1 , ?*J 1 1 


1 , "2t -He 

-3.013£t0o -9 , /bbl.-ye -2.b371.rO/ 

-o.98;t 107 

-1. 

b 1 UL 1 O.i 

-a . i ’ ji * a.j 

-7.1 I-,- TOO 

- 'J v fc U / t ’Of 

*i k h‘/s» » t#f» 

d t *>/ t .rt* 

>i 

2 .t.9r -06 

~3,924E*0o -J .Inat.O/ -2.7o81t07 

- 7. 1 lot 1 67 

- * 

l* 691. tub 

- a , 8 9 3 1 * 0 a 

-7,23b t i Ob 

• 1 • n*MC t « • 

r r ^ <i,t 

^ * I i L r »j O 

i 

3.-*uL-06 

-a. 7b9£t0o - 1 . 5261 1 0 7 -J.9/ot+07 

-o.9<,ot * C 7 

-1! 

6‘7'7t to.l 

-3,/clfto.i 

- o . 0 1 b £ + 0 6 

1 *?!>.< »‘j7 

'j , / , f i -r UH 

8 % J.j ^ . 

a 

•;,2>£-oa 

-j.i'.’bt + Ob -I.ab0t*u7 -3 . 092£t07 

-0,098c +07 

-i . 

<I‘70£tOI) 

• 3 . ooal tOo 

- o , b 8 « i I 0 b 

i . lc<>* / 

c t 0 o 


7 

S . 0 GE-u d 

-6,lb7£*0b -l , iboE 1 07 -3.1‘jiiLt07 

-6.37 o£i 07 

-i. 

'30il*0« 

-2. a8b£ *08 

-3, lbf£*0» 

+u / 

i . S'/#* l T 0 c 

d*C 

7) 

8. 77E-08 

-o,7C7t*07. -1 . 6991 +07 -'3 . 1 7 1 t'.-i 07 

-6.01 IE i(!7 

-i . 

1 32c t Ob 

- 1 .9!te tilo 

-J.22otr «6 


^,»Ur 

1 , 

9 

b.ba£-08 

-7.1b3t*06 -1 .7 1 2E + 0? -3.1991*07 

-b.oEbl: 07 

-9. 

bail 1 0 7 

-l,b2br 1 0 1> 

- 1 .btiaLl (lb 

<? . 7 0 o t 1 0 / 

3 /AK'ctwt 


JO 

7.311-06 

-7,997£*06 -1.7o2E+07 -i.l02Ft07 

-5.237C-07 

-a. 

bb'/l *07 

- 1 ,220E tlib 

-1 ,5 'ml * Oil 

c'.fcWfc tO / 

! * a \ i (. *> ji) 

1 ilib 

1 ! 

fl.06t.-0a 

-7.7aa£*Oo -l.dOSEtO? -j.039F.t07 

-6*6611.0/ 

-7. 

abbt *0 7 

-9 . 9 1 btt (1 7 

-b.b/tl r 0 7 

I . o 1 a t 1 u / 

1 * I* ot.'. ▼ Ob 

1 . OViL tt)b 

12 

d .6b£-oB 

- 7 .r99E*0o - 1 . o921 rO 7 •■2.9o < <i T 07 

-9 .bual - o 7 

-6. 

b >M 1 1 0 7 

-8, )(>7' * 0 7 

- e , aOOt * 07 

1 .?i\ fc •<)/ 

t 1 ♦ T T v / 

<5 , i l « (• r i; ; 

li 

9 .,j2i - Oe 

-7.9?lct06 -1.0661-0/ -2 . dOOL+O 7 

-a. 1 tbc vu 7 


/ivt* 0 7 

n . 7 li ! r. t ij 7 

-b . j I t! r 0 7 

1 » (| i t 1 U’/ 

< .*>/-*» * o / 

o .'• .i/r 

14 

i.vat-j/ 

-7,S68C*06 -1.9alc*07 -s.djittO? 


-b , 

0 6 o £ - 0 7 

-b , 60;it mV 

-s.Obif *0 7 

r v)o 

b . 1 0 / 

j • f / J L ♦ l* 7 

IS 

1 . l2t-07 

-7.70ii.t0tj -2 . 0 1 9E -07 -2.77bE.T07 

-3 , o !2E tO 7 

-a. 

b 1 3t t 0 7 

-a , 7b2£+0 7 

-3.170t to/ 

o»hu5t tOo 

i,l / / tt07 

3 * v d 1 i i v / 

Jo 

1. I9E-07 

- 7 ,7 79t + 06 -2.13! E+07 -2.72(£+u7 

»3,aOi£t07 

-a. 

0 b / E * 0 7 

-a,090£t07 

-J.blbf *0/ 

5 . 7 1 Ot + fcn 

3 .G/^t+y/ 

M . i'iOl i\j l 

17 

l, 27E-07 

-7.bnf + 0c -2 . 20 c.Et 07 -8.6bbl.t07 

-3.262C+C7 

-3. 

6!iaEt07 

-3 , b b 0 1 *07 

-2,2031 r 0 7 

a. 90 l£. + fc’o 

«iS5lr i 0 / 

i.Ct f t / 

18 

1.351-07 

-7.K07t.t06 -2,b8at+0? -2,0961 rO / 

-3.025£t07 

-3. 

3bl r 1 0 7 

-3.102f.r0/ 

-1 .0071* 07 

a too 

^ » */ / 

>1/7 

19 

1 ,a2£-Q 7 

-7 , 1 79E *06 -2.U57E*07 -2.60bt + 07 

-2. 6c9£ * 07 

-3. 

1 a 0 E + 0 / 

-2, 9Obh + 0 7 

-1 . 700Ltu7 

i,b6 f ifct0o 

^ , 31/ ‘*c. ti) / 

d » l £ «l ? 




L 

L 

Z 

2 

2 

l 


i 


L 




3.00E-07 

3.30E-07 

3 , 6 0 £ - 0 7 

3,9ot-0; 

a .2UE-07 

-07 


b , 1 0 1 - 0 ? 

'j , J / 

/On 

-y / 

I 

1 

«tn 

J.nbf -09 

-6, 382h+0a 

-9 . 3 2 a £ + 0 a 

-o. 1 1 3i't(ia 

-9 , bbf)F vfla 

-8,9171 tua 

••/.a 7 nt t«n 


~ r j t i ('4 

-i, i'-.-L tv- 

- J , *j c 7 1 

* +-j 

2 

1 . Ibt-Od 

b, i2.'l£*07 

1 .l-.Obf * C / 

6 . 38 e t tOfi 

i.COoi lOo 

i ,o 2 too 

4. ; i /t tob 

J3 , ti r > S l tltS 

<J . / VO! " *«S 

i.i'i y e /e» i 

is? , tfr • 

•/ 

2 

1 . 02 : -06 

7.9381*07 

2 . 7 a 1 £ + 0 7 

1 , 06 '.i to 7 

babOtC rOtl 

2.9t)'»£ + C ;6 

1 » / o i £ 4 0 » 

1 # 1 b •< t h 0 f) 

6 1 L * u f i 

E a '< »' * L. » fc 5 

- . J > * . 


Cj 

2 , 6 0 f - U a 

9 , 7ui£+07 

.5 , 6 'iot * 0 7 

1 , 8 <ijl.t«/ 

/ . 2 bbf * Oc 

a. ii lit * a o 

i-fio 

i * f m? / 1 >'5+0 

1 . ^4 1 41.- 

, 4 . fc f i V> i 

/ + fc i f i 

»>+ *t 

b 

3.0i'i -06 

UOSaEtOb 

i. 9 .jo >.*07 

1 , 7-.9L t fl 7 

6 . Oibt. t Ot. 

5,0 3/1 * 6 +* 

3 . iiMJf » Ob 

.Mll'hJb 

1 . l>’f c i •* * l> 

o 

'< . ^ - ■> ; 

* • > 

o 

a . 2 lr. -60 

1 . 0 *. 8 f. * 0 n 

a.^rj.'l t .) ! 

1 .Or.atto/ 

1 .nib! *07 

b , V 2vt r I/O 

i , f J S 7^3 i |>C> 

l l* fc 

1 , » i>n 

tyo 


< . t 

7 

b . i/Ol-Ob 

a , d39£ + o7 

a , 1881 *0 7 

2 . 1 a a 1 1 0 7 

1 . 1 b 1 1 1 j 7 

6 , 7 0 0 E r C 6 

*1 • d £ * L < Oo 

M+b 


l • ** ^ 1 f, t \l o 

. * iS 1, 

. t’fc 

o 

b . 7 7 E - 0 d 

9,1, S9t 1 07 

u.3Mf.*0/ 

2.2a7£t07 

1 1 a 2 h 1 0 / 

7.3=7tt66 

t \ , 4 On 

3 ¥ <?/S£ + 0o 

d tOp 

d. iiffcc t(n 

i . o 

+ tf o 

9 

o.but-OC 

8,2 10£ *07 

a. 278tt07 

2. 3 ObttO/ 

1 .310E *07 

7,9001*06 

i.OV5t4<Jo 

tOfc 

♦Uet 

* (jo 

l fc OOPf 

T‘'T 

10 

7.311-08 

7. J/bt+07 

a , 1 1 1 1 * 0 7 

8 . 3 1 b £ 1 0 7 

l.ib/EtCV 

3,3311 rOb 

5 • ^iSMOb 

3 • cJ7 Vt4 Oo 

i,iv7L^vt‘ 

f.&VH 

1 iOh 

o 

1 1 

‘il.Obt-Od 

6 . bOOt t 0 / 

3.912t+07 

2.3081*07 

1 , 3boi, 1 0 / 

b , 6 S c F. t 0 6 

*2 * 7 1 5 h t 0 o 

^4 1 > 1 ? S t i 0 b 

i * 4b 1 £ + Ot. 

2 • ** *7 > c + o fc 

li ‘f 

c,Mj 

1 M 

o. 6 ‘,t -0 8 

b.t /ot * 0 / 

3.o9hEj 0 7 

8 . 8 o ? f. * 0 7 

1 .3981. + C7 

O,c7‘/k. + 0c 

5*43 .£.1 {,*> 

« # 3 ^ 5 1 ♦ ;i cj 

i * r >'/5' ▼Ot 

5 I Of i » i* o 

.'J.'ll 

*» 

1 i 

0 ,o 2 t -Ob 

b , c 92c* 0 7 

3,8711 r 0 7 

8 , 2 0 b E * 0 7 

1 . 1 Ool. t C 7 

9.010 r, + 6 6 

b » I'tisiE. ijo 

• 1 , ** t < Utj 

i # c y r \ 4 > 

>t »<.'.• 

ji . j u : i 

• s, n 

1 K 

1 . U a F - 0 7 

«. 6692*07 

3 . ob Jt * r 7 

2.1 3 « L 1 0 7 

1 Obiil+f, 7 

O.l'bfcetf.t. 

b. IbUz fc Oo 

aJ , fc 5 i t’ 't v fc 

^ „ ^‘/ciL f V 1 . 

3 fc 0 J.'l T V t 

A *<//?! 

f , i , (j 

15 

1 ,121.-07 

a ,« iafc+07 

3, Os.’.h t ( / 

2 , Ob VL t 07 

l . ioOfc + C/ 

0 , 02 bL *06 

o , tflr/ 

4. ?4tLi (rt» 

'+ , * b J 2 i o o 

fc* , ; 'i jL *\»it 

4 , M 71 

VlfO 

l b 

1 . 1 Ot -0 7 

3.8 111. *07 

2 . 8 S 2 L U17 

1.9/96+0/ 

1 . 33u£t07 

f! ,98bt + 06 

o # ? 0 5 £ » Ob 

^ . 7 3Vi * <;o 

« , c * } b fc ■» 4 + 1> 

•4 , .»V-t t On 

b. Hot. 

rfco 

* 7 

1 .271.-07 

3.a7aE*07 

2. o 7of. * 0 7 

S.OOoE+07 

1 ,201 C t07 

b , 7b6t+0o 

b* 1 liLtOb 

H « / / 7 t if 0 b 

* . it + un 

A * *i 4 o 4 t ■„ 

b, Ifcfcl 

T fc O 

lc 

1 .35r.-07 

3 . 1 9 /£* 07 

2.5lol+07 

1.81 a8*07 

1 , 2L Ot 1 0 7 

8 , 53ife * 06 

6,C10fcfCb 

4 , t /.oe * Ob 

L.*n t*uo 

tvo 

7 , S b o t 


! 9 

1 ,a8E~07 

c. 975£t07 

2, 3/9E+07 

1 , 73 5Lt07 

l,2O3t+07 

B,2b3i.*C6 

5,t/i /t + 0 o 

il , 5 fc t i 0 6 

9 6 iioc * yfc 

5 4 fco 

s , 43^h i’Uo 



A5-14 


1 

i.ebi -00 

• . 0 1, c - 0 7 
-i' , e 1 7{ <0/ 


■ . iSi' -<11> 

i . 3 1 1 E-vilb 

3 

! ,'<£c-0B 


•1 

d , oOc. -06 

0.bl2f*0S 

5 

5 , *1 ci t - u i 

b.ba SfctOS 

o 

>i , r ' <; ti 

o.b! 71 tbo 

? 

3, .(i- ••!,!} 

'/ , •l./l* / V O'! 

c 

i.77 l -0<1 

b^uitOb 

■) 

b ,iuL-Gb 

*) , QbbitbS 

10 

7 ,ilt-06 

9.77U + 0S 

1 1 

o,v6t-0B 

i .Obae^oo 

12 

b.GDC* O'! 

l . K)2Eti)D 

1 3 

9,< Si-Jb 

1 , 1 ‘5b r rOo 

l- 

i . CiiL-'lV 

: , i/Oit < Co 

• «. 
« -> 

1 . \il-bl 

1 , t 9 / t. i Oo 

1 O 

1 .Hi-07 

. 1.2«5£+0o 

17 

1 .27f “07 

1.31«!£*06 

Id 

1.3bt-07 

1 ,39/it Oo 

10 

1 .H2E.-U7 

l , 3 7 1 it Oo 





*/ 

l 

J 

/ 

l 

/ r 

7 

.• 

. 

/ 




0. 

1 , Oof - )i: 

l+.OOt-O+t 

9 , on -oa 

1 , t '[-07 

1 , SO * -9 7 

1 .211 -0 7 

2. ! in -i. 7 

(' - - 'l *-> 1 


H * 

1 

t 

«<n 

i.ftbC-09 

a.osotto? 

4. Sob! *07 

9. 6209+0 7 

l .9(, 56 + 06 

3, 40 1 6 + OH 

a . Itbbf.i 06 

-1,11 4L+ 09 

- 1 . 00o6 + 09 

4 , 1 / r t l > yd 

3 ,ii',f 

t •iti 

a 

1 . lSt-08 

2.4*>9t.07 

4 . HbSFi *07 

9.8201+07 

1 .9356 + 02, 

1. 00 1 £ + 06 

2.I1SHL + 0O 

-1,1 146+09 

- 1 .dJut+09 

*1. 1 /9L fOb 

1.3341 

iOn 

3 

i ,'iai -oa 

a, 31 3L + 07 

4.470E+07 

9.iaof >07 

1 * ftltab + 08 

3,0196,06 

a . 3 r at <oa 

-9.90 76 >06 

-4, Old +06 

Kl O 

3 . U c. 4 L 1 0 1" 

A 

a » o J u -orf 

a.l»l!)Ltu7 

a , 0 7»t * 07 

7.99ii <07 

1 .‘>1 7t < 08 

2.0.*SI.< 01' 

1 .‘,921 +06 

-8, ‘.‘SMt + t'O 

-/."'HI *i.c 

^ V »t i (M 

/ 1 *» ► 

*i*n 

s 

J.-ftt-OJ 

1 . 9421+07 

i.ebbt *07 

6 . 1+ 1> 1 C < 0 7 

1,2306 + 01. 

1 . 7‘<3:. + .''i. 

li , 3 7 ,1 + !> 7 

"o,37ii.+ Ob 

- S , 7 b - 1 + •• 11 

l , / ' *1 ' 

f'.Vt'Mt 


0 

4 . ait -06 

l , / 24 E + 0 / 

3.2336*0/ 

S.ooZL+O 7 

9,77 ol +07 

1 ,2‘>9tt0t> 

a. 1,91,6 + 07 

- 4 , 7 +1 7 1 +16 

-4.2001 + '+8 

1 , ') 7 o t mm 

i , *><* i i 

» .iu 

• 7 

b. 006-08 

i ,<m£*o/ 

2,8156+07 

a.9iet*07 

7 a o 1 3t ♦ 0 7 

11,6106 tC 7 

-1.0936+07 

-i.b'j/t +06 

-3, lb It + 06 

J , - 1 t 7 V / 

l . ) *!>=. 


a 

b,/7t-0U 

1 .24/E+97 

2.40/Eto/ 

a . 1 1 at<07 

b,9jyt+07 

b.oiS6+ 6 7 

-a. 9246 tit/ 

-2.o3bL+oe 

-a. 3J06+08 

<i, OtVVt T^/ 

o * ; f i « t 

to / 

9 

6,b-)F -06 

9.9676+Ob 

2.017tt07 

3.3849+0/ 

'1 .b4o6+07 

3,62/6+07 

-3.7036+07 

-1 ,4751 + 08 

- 1 « 7o2t + 08 

b.cdM rVb 


tO/ 

1 0 

7 ,3lt-08 

7,4106+06 

l.n49E+07 

2.7091+07 

3 , 4oo£ <07 

2.a39E+07 

-3.8b46+07 

-1 .bfl06 + 0fa 

-1,3441 +615 

- 9 * c t A. o L. i 0 o 

4 * ydot tO / 

1 1 

t.Oi.l -OS 

0 ,r 1 46 t Ob 

t.SGSl <07 

2,2026 <o / 

t.biaii+ov 

1.2996 + 1.7 

- 3 , 6 1 4 f + 07 

- 1 , l S06 < 6b 

- 1 . UJbEtOb 

- 1 y M VJl. + U / 

el * b i c’r. + y +' 

I -■ 

6. t-'+b-'jO 

2,1 7tlT0b 

O.rthCfc .Ob 

1 . 7 v S 1 1 0 7 

1.9996(0/ 

(i ,0901. + On 

-3,33. ‘ +!'/ 

-r. . 61,9c + 0 7 

- 5 , 0 3 bf + 0 7 

- i . ob y| ) u 7 

1 .■> Sv + 

♦ o; 

! 3 

o.tat-oo 

-5 ,cibt <06 

o . 90 06 < 0 b 

1 . 3aolt07 

1 , 4 3 c t < 0 7 

2,6726(00 

-2.919c + i)7 

-c. .bVbt <67 

- 5 ft +. + U 7 

- 1 • <■•<' If * V 1 

, . ,4 -f 

<;•/ 

, 

l.0*.t-07 

-3«31‘+6 + y6 

•3 , 3926 *0» 

1.0106+07 

1.0376 + 0 7 

8,0936+04 

-2,4bd, + 0/ 

-b,3<4 l<"/ 

-4 .8 7'.L > - 7 

** 1 y *4 « / £ ♦ «/ 1 

8 . nb* ; 

i \>d 

is 

i.iat-o/ 

-o.auottou 

2. 1 Oof t O b 

7«3i>o6t0o 

/ • 3b8E*Ut> 

-1.30cE+6b 

-a.06bt *07 

-a. 1241+07 

-3, 7c Lt + 67 

~1 , r v < i, 1 0 / 

b . r> r c. i 

♦ Oc 

lb 

1 . 1 9L-07 

~9,4OtiE + 06 

1 .oftlEtOS 

S. 1 6fln + 0b 

S. 0996+06 

~l.913l.tOo 

-l.oott +0/ 

-3 , 1 3br +o7 

-2.1-096 + 0/ 

- 1 * Ubrc. r 0 / 

3, 1*541 

i Ou 

1 7 

i .a/t-ov 

-i .Zfl/6+07 

••1 , 3 S 96 <06 

3.4bl6tOo 

3, (12 lfc<0b 

-1.99 71 ,06 

- I ,21136+0 7 

-2,31bi +0/ 

-2, 1259+07 

*bb-.l + vio 

2.1731 

+ C*o 

lb 

1.35E-07 

-1 ,6 /9t + 0 7 

-a,3b9t<i)a 

a. ia9t+oo 

2.1 l‘b| + Qo 

-1 .7306+06 

-9. 3026+06 

-l,o + 7 L + C* 7 

-1 .4991 + 0 7 

- c . S b V i ioU 

1 . 2b’+i 

t Oo 

1 4 

i .oat-o/ 

-2.13ot*07 

-2 . b69E + 06 

1 . Io3t'<06 

i .a/ottoo 

-1 ,245E + u(> 

-6 . 0o8f + Ob 

-1.6286+0/ 

-4 , 4b,;c + 06 

l?U +VO 

i, .« Sot 

tO‘j 

ao 

1 ,506-0 7 

-2.605E+07 

-a.uaac.+06 

a.ouae+06 

b.7/46<05 

-o.aiit' + 'ib 

-2.911+6 + 06 

-4.92+6+06 

- 4 , 6341 +06 

-r.oibt *U6 

2.9201 

'Or* 





2 

l 

l 

2 

L 

Z 

/ 

/ 


7 




3.0C6-07 

3.3OE-07 

3,606-07 

3.90t -07 

4.206-07 

a.b0f-07 

«,bbt~b7 

b » 1 (M -07 

b,^or -0/ 

5, /ot “0 / 

o o 

I 

tn i j 












1 

i.fcbr-09 

1 . 4 11 b C + Ob 

7,0346+0/ 

3,o6cL +07 

. t 1 S ’ i v. 4 0 / 

1 ,U‘,t+)/ 

. iltySk * y 0 


if J« i+E + Cb 

**»:*• i» i‘ v Of. 

l./t’- 1 31 7 


5 

1 . lbf-Ort 

1 .4ftc6 + 0o 

7 . 0 3aE + 0 7 

l.cbbf+Oy 

el * V^‘jC t 7 

c , lu5f.+ n 

6 * dOt L + 1 lib 

tl t vOo 

* 2 , 3 7 \ f t y b 

-t t Ab vl t 00 

“• i . ^ « +»Jf» 

^ M 

? 

i .art -06 

1 . J9nfc + 0H 

o . 7 iol +07 

-,.S/9r. + 0? 

*»‘ic + (+ / 

1 . S«-f 1 J 7 

b « 1 t bf* 

7 . J l 0 «*t r»iO 

v f* % t Ob 

" f » i" ▼ i* c 

•» 5 , 1 3<'i c + 1 ' C» 



2 . o9i -Ob 

1 , 2bOL+ 08 

b . 5 1 tl c + 1+ 7 

5 . - J 46+ 67 

1 tij7 

1, 1 141 +77 

b/r/jf t Or 

/ y U V* ;*< * +. 'J 

- , Ww i £ t vb 

- c 1 « *Jf i r ■* ^ 

" * 'y V t-c.^ 

a* 

u 

3.9b! -01, 

1 .0942* 6b 

b.oO+t +u? 

3.2446*07 

1 > cl / - c. -t 0 / 

1 ,9/Vt til 

S . 7oot tyc 

<* . 6 r > t < t> 

** > » V> » “b 

t + + f. 

" “ . t' / 7 l Y J*. 

v ir 3 

o 

- .231 -Ob 

4, 3C|4E +0/ 

b , 2 3 + f +07 

3,0226+07 

i ,77^1 +07 

1 . CiOf. <07 

b . b 1 4 C 1 tr f» 

2 j J i»<‘l v Oh 

. Ww -»0t> 

-t'.Si'-LMt? 

" '< , 1 0 . J ' « •/ *1 

WT) 

7 

b.i.Oi -CH 

7.7296+07 

4.ob2t +0/ 

6. 7746+07 

1 « oc.<d t 0 7 

9.7.- It -06 

b , k 1 ft t + »f 0 

2,00/+ 1 0 0 

-b.AiiL+yb 

"C.ti/iftityo 

*• •* . i t S* t *.o 


<i 

b . 7 7 1 •• 0 8 

b.330E + C+7 

a, 07H6+07 

2 . S2S6 + 0 7 

1 > U jdL r 07 

9 .(){'> l <0b 

+♦ # + ub 

1 y « > l Y 0 li 

" 0 . U> 1 i *» b b 

"i . /l^utOci 

- J , ‘J f' 1 { t y 1 . 



o . l , «6-Ob 

b. 1146 + 07 

3 , S 3 3 E +0 7 

2.2676,07 

1 

1, , 4 0 3 6 *06 

■i| , SiMlf + (/t* 

1 , {> 4 / l t'7 0 

- / * (r ( + 0*> 

- i . *) ) l l T 0 0 

" f J , a e* -.1 J T \J ( 

f 

to 

7 . 3 1 1' - 0 5 

a , i 2 if + 07 

3.0b7c + 0 7 

2.0136 >07 

i « 7. 7u£»C7 

7 . o 6 4 1 +66 

0 , 1 3or ^ Otv 

1 . <ib 51 + 0 <*• 

** L , b 

-i , o s bi vuc. 

^ il.L * V<> 

; : 

i* . Or c - oft 

3,2 65-.+ 6 7 

a * s +9/ 

l.t (,9c + 07 

loll.uO/ 

6 , 9 .i 1 !. + , 6 

\.rs 7L-VC+ 

3 , st* 

M l> ♦ Of' <• t "/ b 

- J i ! 

"b.f 27c J 

1 7Z 

I : 

+, , ob' -Oo 

6,:98i<07 

2.1 6,69 + 07 

1 . 6 3 1 £ + 0 7 

1 * v 0 7 i t 0 7 

c , lb 76 - 0+. 


1 > 

w i , Ob • C 1 1 >0 

" : . ^ * 2 : 1 y*/o 

• j # ‘> * j J- + ; w 

■ C/5 

* i 

9.621 - 0 4 

2.0326+6/ 

i , Vl- c,t + 0 7 

1.5106+0/ 

o * /OVutOb 

6 .4 i2C**.h 


. .vj /[ v i,b 

- 1 , \ * Cl - L* O 

■* i g ’ r f . T »K. 

k fi,vi *1 d' 


1 ^ 

t .046-0/ 

1 ,So4tt07 

1 ,41.31 +07 

1 . 1 029+0 7 

7.*H7:. + (70 

4 .ofcftl + i. b 

2 « b 1 ■+ 1»*+> 

o.V^bl +' b 

- 1 . it*«c »*jc 

■ j • b >' / 1 'Oo 




'i.iar-c7 

1 , 2 1 ct + 07 

1.1776+47 

9,09 5£ + 0o 

o* H /it » Oo 

3 , 94 46 r 06 

2 , 1 i2L+*0o 

b , i \ ot t uS 

"1 » 

" 3 « t 'jbt v .>0 

• / » i ? •* i ♦ 0 y 


U 

l.SVt-07 

9,lb9t+0o 

9.2481 * 0 6 

7 . 3 0 .J c <06 

S* 0 % c > h +0c» 

’3 , 2 2 6 T < 0 0 

1 . 1c Ah + (*y 

1% . V 6 b ! r 0 b 

M , lblfcU'o 

" i « <' ) i l +l!j) 

" ! *yt> 


17 

1.276-0? 

b . 9 4 1. .• 0 o 

7. (1026 + 1+b 

5 , b 3 c t *0+i 

7 =t + 0<> 

2,6319 + 6+. 

1 , id/r+vo 

>. w 7 : / L ♦ C b 

•• » , ob*i J + 0t» 

"J.il/I +'/o 

*■ 0 , v 1 b 1 1 1 * 0 


ta 

1 .381.-07 

4 , s52E + 06 

b , OOn* tub 

4 , 0 6 7 c + 1+ b 

2,‘v0v£t0o 

1 , 0686 + Ot 

9,bbit >Ct> 

1 ♦ *’t»t Ltu5. 

-0 .h'yit * Ob 

- i . /. •> '* r. + 1 0 

yOi^t. +U7 


1^5 

1 .«2l.-07 

2. L , 126 + 06 

3 ,2a 3l. + Ob 

2,1.411 < i; 6 

t 

1 .2121 ’ 1 Of, 

0 . ^ | ti i“ * v b 

1 . 

- & ► bu Ac tOb 

" t « O - ^ V Y / O 

W 1 , 1 ff 6 ’1/ 


bM 

1 .but-07 

1 , 3Sjl.6<’0t 

1 , 6 , 4 ft 6 + 0 u 

1 .bUuE+Ob 

^•c<! jf-rCb 

S.424C+ 65 

A , 1 t Ob 

« ,b‘»dl tiW) 

-i t0*3 

- i 4 fi Hi + 0 j 

•* t . Sob* + J? 



A5-16 


i'Loiriv'L CHiHuL ICC.'- - . 1 d.'.b 


0 

0 

0 

3 

3 

46 

16 (. 


'lb 

u 

2 

0 

0 

0 

0 

0 

2 

7 

6 

«2 

129 

ft5b 

7Vj 

1 lift 

10 

2 

0 

0 

0 

• 0 

0 

0 

0 

21 

61 

iao 

6911005 

105 

ift 

0 

0 

0 

0 

0 

0 

0 

0 

II 

30 

ifl" 

55o 

892 

103 

15 

0 

c 

0 

0 

0 

0 

c 

0 

c 

17 

7 « 

3o9 

59 S 

1 Oft 

10 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

3ft 

17/ 

‘319 

,77 

ft 

0 

0 

0 

0 

ft 

V 1 

0 

0 

u 

0 

2 

65 

12b 

3/ 

1 

0 

0 

c 

0 

0 

0 

0 

0 

0 

0 

0 

11 

35 

ft 

1) 

0 

0 

6 

0 

a 

0 

0 

0 

0 

0 

0 

0 

2 

1 

1 

0 

0 

(J 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

a 

0 

0 

0 

0 

0 

ft 

0 

0 

0 

0 

0 

0 

0 

0 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

(J 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a 

0 

0 

0 

(1 

0 

0 

c' 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

c 

0 

0 

0 

0 

0 

0 

0 

- 0 

0 

0 

0 

0 

6 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

c 

0 

a 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


iitC.A r I V£. CHAKUE LQCAUUN& 

6 o 0 0 0 0 30 b92 Oft i> ft 0 0 0 0 

0 0 0 0 0 0 1 fj o 1 120 2iu 2 0 ft 0 Cl U ■ 0 

0 0 0 0 0 0 22ol<fH0 iftft 22 0 0 -0 0 0 

0 0 0 0 0 0 1 oO ! 306 320 2o 0 0, 0 0 0 

(i o o o o o ho a<ia no 2000000 
0 0 0 0 0 0 72 fttll) i 30 0 0 0 0 0 0 

0 0 0 0 O 0 30 2lM CO 2 0 0 0 0 0 

000000o52t600Uft0 
(1 000000022 0 0000 

<i l U 0 0 0 0 2 Cl 0 0 0 0 0 0 

0 c 0 0 0 0 ft 0 0 0 0 0 0 0 0 

80000 0 0 00000000 
OOOCOOOOODOOOoO 
wOOOOOOOOOOOOOO 
OOOOOODOOOOOOOO 
0 0 0 0 000 OOOOoOOO 

COOOOOOOOOOOOOO 

0 ' 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

OOC U 0000 V 040000 


0 0 0 0 o 
0 0 0 0 0 
0 0 0 0 0 
0 0 0 0 0 
0 0 0 0 0 
I! 0 0 0 0 
OOOO u 
v 0 0 0 0 
0 0 0 0 0 
0 0 0 0 0 
00000 
0 0 0 0 O' 

00000 
0 0 0 c 0 
0 0 0 0 0 
0 0 0 e u 
0 0 ' 0 0 0 
0 0 0 0 0 
0 0 0 0 0 


00000 
0 0 0 0 0 
0 0 0 0 0 
0 0 0 0 0 
0 0 0 0 0 
00000 
0 0 0 0 0 
0 0 0 u V 

00000 
u y 0 0 0 
00000 
0 (- 0 0 0 
0 0 0 Ci 0 

00000 
0 0 0 0 0 
0 0 ■ 0 0 0 
0 0 u a 0 
0 0 0 0 0 
0 0 0 0 0 



Appendix 6 


DESCRIPTION OF ACORN 

I. PURPOSE 

ACORN is a computer code to calculate streamer 

development initiated by electron avalanche. 

II. 1 START-UP- REQUIREMENT - INPUTS 

The following inputs are needed for an ACORN run: 

CARD1 C3I5) — contains three integer quantities: 

KM-T = 1 — for starting from cards with new data. 

= 0- — for restarting from magnetic tape 
(TAPE 3) at dump number KNR. 

KNR — writing on TAPE 3 commences at dump 

number KNR.+1 . 

KNL =0 — if no enlargement of grid, spacing 

is required. 

=1 — if double spacing is required for 

consecutive run, 

CARD 2 C3I5) — contains three integer quantities: 

NR — number of radial mesh points covering 
discharge volume. 

NZ — number of axial mesh points. 

NRR — number of radial mesh points used in 
calculation of potential.. 
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CARDS 3-9 (1P3E15.8) — contains parameters for- the 

dielectric avalanche. 

o 

P — gas pressure. Newton/m — Vestigal. 

GAP — imaginary gap length, M. 

RAD — maximum radius at which charge density 
will be calculated, M. 

FRA — fraction of distance between initial 
charge cloud origin and left-hand 
boundary to the gap length. 

EPS — dielectric constant. 

VG — voltage across the gap due to external 
field-. 

ALC1 — constants in formula for the primary ionization 

- ALC3 coefficient, a: 

a/P=ALCl*EXP [-(ALC2) / (E/P) ] + ALC3 . • 

DELTC'1 — constant for detachment (for gas- run 
purpose only).. 

WPC1 — constants for the ion drift velocity; 

- W PC 2 Vp = [WPG1 (E/P ) + WPC2] • E/E. 

WMC1 — constants for electron drift velocity: 

- WMC3 V e = [WMC1+WMC2/(E/P)I • E/P for E/P > WMC3 ; 

= const ^ [l-const 2 (E/P) ] • E/P for E/P < WMC3 . 

(The constants are chosen so that the magnitude 
and slope of V e are continuous at E/P = WMC3 . ) 

D —.material specific heat, J/M 3 - °K. 

ZB — initial charge cloud length, M. 

RML — avalanche length at the beginning of the 
calculation, M. 

TOLPOT — fractional accuracy required in iterations 
for potential distortion. 

TOLWM — fractional accuracy required in iterations 
for position of foot of characteristics. 

DTMIN — minimum allowed time step. 
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CARD 10 (415) 


contains four integers. 


NT 

NTR 


I D I V 
J D IV 


number of detailed printouts between 
successive dumps to TAPE 3. 

number of time steps between detailed 
printouts. (NT*NTR is therefore the 
number of time steps between dumps.) 

radial spacing of 'arrays in print-out » 

axial spacing of arrays in print-out 
(must be a power of two) . 


After each dump to TAPE 3, these four quantities 
are read again, as many times as desired. 


LAST CARD — A blank card, or four zeros, in (415) 
terminates the run. 


II. 2 START-UP REQUIREMENTS 


A magnetic tape is required to store the data for 
consecutive runs. For test runs, tape is not required. 

III.1 CODE OUTPUTS - PRINTED OUTPUT 

The code prints out first the input data, followed 
by a listing of the radii of radial mesh points. Next, the 
current values of NT, NTR, ID IV, and JDIV are printed. 

At the end of each time step, a summary of the 
results is given comprising the present values of the time, 
the time step, and the number of iterations. 
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Every NTR time steps the following more detailed 

th 

information is printed f.or every IDIV radial mesh point 

X l. - 

and every JDIv axial mesh point (in addition to the above 
summary ) : 


a) 

b) 

c) 

d) 

e) 

f) 

g) 

h) 

i) 


axial electric field, V/m. 
radial electric field, V/m 

3 

electron charge density, C/M . 

3 

positive ion charge density, C/M . 

3 

net charge density, C/M . 

temperature, °K. 

total electron charge, C, 

'dipole length, M. 
electron charge center, M. 


After NTR*NT time steps the record number of the 
dump to TAPE 3 is printed, followed by the next set of values 
of NT, NTR, IDIV, and JDIV. 

The above sequence of output is repeated until a 
data card is encountered with NT=0. 

III. 2 CODE OUTPUTS - MAGNETIC TAPE 3 

The time, time step, order of enlargement, electron 
density, positive ion density, net charge density, and tempera- 
ture up to the present time are dumped to TAPE 3 in the form of 
an information record. The calculations may then be restarted 
at a later date. 
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FIGURE A6.1. GEOMETRY OF CODE ACORN 


PROGRAM LOGIC 


The logic flow in ACORN proceeds as follows 


SET UP GRIDS AMD READ INPUT' PARAMETER 


READ DATA ON TAPE 3 IF RESTART THE RUN 
AND ADJUST THE MESH SPACING AS NECESSARY 


DEFINE THE INITIAL CHARGE CLOUD 
IF START FROM MEW DATA 


CALCULATE THE BOUNDARY POTENTIAL 


CALCULATE POTENTIAL, FIELD, AND DRIFT VELOCITY 


CALCULATE THE DIFFUSION EFFECT 


INTEGRATE ALONG ELECTRON CHARACTERISTICS 


INTEGRATION OF THE NET-DENSITY EQUATION. 
CALCULATE NEW ELECTRON DENSITY, POSITIVE ION 
DENSITY, AND NET CHARGE DENSITY. 
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TEST RUN OUTPUT 


TEIPORAL IROITH OP AVAL I NCI E- - TIO-IIHENSIOIAL CAICULATIIH 


IHFUT KITft 


KHI « I, KNI * 0, INL * 0 


(* • 2 t, HI • 43 , «Rt * 24 


I.OOIOtOOIE’OO 

3.004040IOE-07 

1 .54I0400IE-07 

2.00100 IOIE-1! 

i.77toioioe-n 

4.00I0040I£-*02 

l.aSI00404E»08 

7.22404OO4£«OS 

0 . 

0. 

0. 

0 . 

I.40I0040IE-04 

0. 

1 . OOICOIOIE*04 
2.3OI04OIOE1O4 
1 .06IQI004E-07 
1.33IOIOOIE-OS 
l.OOIOIOIOE-OI 
1.00I0040IE-I2 
5.00I0400IE-18 


GJS FRESSIAE (TORI) 

GIP LEIGTI tl) 

RIBIIS OF OIJCIARIE (It) 

FIACI10I IF IISIAICE OF IIITIAL CIARIE CLIUB Tl TIE SAP 
BIELECSRIC CON STINT 
VILTiGE ACROSS (HE GIP 
) 

) CINITAITS II FIR4UL4 FIR FRIHAtT IOIISITION COEFFICIENT 

> 

CINSIAIT IN FOIKULA FOrEICITATIOl CIEFFICIENT 
) 

) CINSTAITS IN F8RRULI FIR 101 BIIF1 VELICITT 
) 

) CONSTAITS II F0RBUL4 FIR ELECTION BUTT VELBCIIY 
) 

BEHSITI OF TIE SIUJ NITEII1U XG7H3 ) 

AIIAL IANIE IF INITIAL CNIRGC CLQIB I I I 
AIALINCHE LEIGTH I I ) 

FIACTIINAl ACCIRICT RE1UIIE) II PITEITIAL 

FIACIIQIAl ACCUIACT IEIUIIEB IN SLOPE OF CHAIACTEIISIIC 

HIHIHUI ALLOIEI VILVE OF THE STEF 


1.00I00404E*00 OSIER OF EILAIGEIEIT 


RIBIIL NESH 


I 

i. 

2 

7 .I 004 E-I? 

3 

i .sioiE-ia 

A 

1 . 340 IE- 4 S 

5 

3 . 1 I 0 IE -43 

6 

j.jtoiE-ia 

7 

I.I 80 IE-I 3 

8 

5 . 4 I 04 E-I 3 

9 

4 . 240 IE-I 8 

10 

7 .I 20 IE-IB 

n 

7 . 8 IOIE -43 

12 

1.580 IE- 1 8 

13 

1 . 340 IE -1 a 

14 

1 . 0144 E-I 7 

IS 

1 . 0 I 2 IE-I 7 

16 

1 . 170 IE-I 7 

17 

I.H 84 E-I 7 

ia 

I. 324 IE-I 7 

19 

I. 4044 E -47 

20 

I.I 824 E -47 

21 

1 .S 404 E-I 7 

22 

I. 438 IE-I 7 

23 

I. 7 I 4 IE-I 7 

24 

I. 7744 E-I 7 

23 

I.I 720 E-I 7 

26 

1 .950 IE -47 



5T4RT OF CALCULATION 


IIIT1AL CONIITIONS * 

THE • SEC t SAP VOLTAGE - 4.IO0OOE+O2 VOLTS 

NT * 1 HTR - 5 IBIV * 2 JD1V • « 


TIIE 

SEC 

1.828738-t 4 
UI9726E-H 
2.49740E-H 
4.64U0E-H 
4.641IOE-14 


TIRE 5TE7 
SEC 

1.I2I73E-14 
6.85256E-I6 
4.00340E-15 
2, 14370E-1 4 
5.0COOOE-1I 


NUMBER 

ITEIATIOHS 

2 

2 

2 

2 

2 


DISCHAIGE VARIABLES AT TIIE 4.44I80E-H SEC* 


AIIAL ELECTIIC FIELD HOLTS/ rt) 


I.01000E+08 
8. O0?78E*O8 
8.009T8E+08 
8- 00829008 
8.00724E*08 
I.0*6I5E*Q8 
S.fl051OE*OB 
8.004 ME+O# 
8.003318*08 
I.C0262E+O8 
I.00203E+08 


1.13104E-1I 

7.94961E-20 

9.30737E-23 

t.l9808E-25 

1.73934E-28 

2.T2389E-31 

3.73860E-34 

4.04025E-37 

3.7J438E-40 

3.04380E-43 

9.4I539E-48 


1.0381 98+00 
8.O3573E+08' 
8.02997£*O| 
8.023358+08 
9.O1724E+09 
I.01223E+GS 
S.OO834E+00 
8.00551 E+Ot 
1,003498*08 
I.04210E+08 
0.001 UE+48 


5.91 1 57E-07 
2.9423IE-08 
2.51 407E-I 1 
3.09520E-14 
5.02299E-17 
7.54480E-2O 
9.34302E-23 
1.01 4 JOE-25 
9.2I879E-2* 
7.46432E-32 
2.039418-33 


8.2O483E+08 
8. 180648+08 
I.09970E+08 
S,O4873E+08 
8.Q21528+I8 
8.OO760E+48 
8.OO085E+O* 
7.997128*01 
7.99663E+I8 
7.9963IE+08 
7i99453E+08 


7.741 13E+01 
2.074148+01 
8.J0813E-01 
1 .048428-42 
5.43O8IE-0S 
1.5809TE-O7 
Z.95044E-10 
4.010828-13 
4. 302918-16 
3.I3100E-I9 
t.09672E-20 


7.33I2OE+08 
7. 46359E+0* 
7.690248*08 
7. 8311 2E+C0 
7. 70221 E+OI 
7.94 039E+08 
7.96190E+0I 
7.9745IE+0I 
7.9I223E+08 
7.9I717E+08 
7.99H5E+0I 


2.92127C+05 
2. 081 49E+05 

S. 38072E+04 
6.40I76E+03 
1.54822E*02 
9.12432E+0O 
1 .215678-01 
9.JO4240E-I4 
4.0831 2E-04 

T. I1S43E-OI 
5. 1523 IE- 10 


8.34332E+I8 
I.28535E+08 
8.1541 28+16 
8. 55 1 84E+06 
I. 00490E+08 
7. 99035E+08 
7.187I6E*08 
7.918848*08 
7.99I59E+88 
7.99225E+O0 
7 . 19364E+0I 


0 . - 

1.47027E+I7 
2.221 27E*07 
1.75328E+I7 
1 . 1 II77E+07 
4.483248*04 
3.l82t7E*Q6 
2.422168+06 
1.574448*06 
1 . 061 13E+06 
7.310338+05 


5.447798+04 
4.64USE+04 
2.130168*04 
5.96342E+03 
1.001 60E+03 
1 .823871+02 
4.415538+10 
2. 4732 IE-81 
5.96352E-03 
9.1381SE-05 
1 .085798-85 


8. 853148+01 
I. 151 938+01 
8.0435TE+0I 
0.I33OOE+OI 
I.82293E+0I 
8.114808+01 
8. IO893E+0I 
I.08494E+0I 
0 .002428+01 
8.000838+01 
7. f 99I8E+0I 


0. 

1 .091248+04 
1 .908378+04 
2,299598+04 
2.32052E+04 
2.1 4964E+04 
1 ,800038+04 
I.4790IE+04 
I - 1 9007E+04 
9. 47642E+05 
7.31874E+05 


2. I 69 998-01 
1 .199488-01 
1 . J 126 IE-01 
4.457368-02 
1.23I62E-02 
2.48301E-03 
3.429698-04 
3.I2305E-05 
2. 263488-04 
I .016968-07 
2.164288-01 


I.II428E+OI 
8.013938*01 
1.01294 £+01 
8.11 148E+0I 
I. 80976E+0I 
I.IO798E+0I 
I.I0631E+OI 
I.0O4I3E+OI 
I.0O339E+OI 
I.08258E+OI 
I.08179E+0I 


8. 

1.90055E+05 
3.57849E+05 
4.04912E+O5 
5.700538*05 
4.48527E+05 
4.09877E+05 
5.84 1 49E+03 
5.41 327E+03 
4.I9593E+03 
4. 350 1 5E +05 


4.037058-11 
J.44477C-11 
2.364I3E-1I 
1.I3874E-1J 
4. 10451 E-12 
1.110108-12 
2.259138*13 
3.469I2E-14 
4.03452E-T5 
3.54340£*1 4 
9.I37J3E-17 


t. 085408+08 
1.005538*08 
S.O0531E+O8 
1.004978*08 
1.I0454E*08 
1.004058*08 
1.103538*08 
l.083#2E*08 
1.102538*08 
1.082088*08 
1.001468*08 


I. 

5.545478*04 

I. 077298*05 

J. 531058*05 
I.S96S1E+05 
2. I6320E+0S 
2.331138*05 
2.4084OE*O5 
2-'4fl921E*05 
2.348988*05 
2.244038*05 


4.O701OE-23 
3.844038-23 
2.552448-23 
I .284078-23 
4. 92544E-24 
1.44534E-24 
3.2595TE-25 
5.47440E-26 
7.46606E-27 
I.0I727E-28 
2.43079E-28 


RADIAL ELECTRIC flELl (VOLTS/ M) 

8 . ». 0 . 0 . 

3- -7.47435E+OS -1.6460IE+07 -3.O4233E+07 

-1.25811 E+04 -1.29524E+07 -3.086128*07 

0- -1.49629E*0 6 -9.152458*06 -1.61495E+07 

1 -?I.33092E+04 -4.452548*04 -7.119148*06 

*• -1.44I44E+04 -4.51 555E+04 -3.260348*06 

*• -1.28T84E+04 -3.15918E+04 -1.442458*06 

-1.11726E+86 -2.224 9 4 E+04 -9.130728*05 

0- -9.47771 E+05 -1.38458E*04 -5.33412E*05 

8. -7.934998+05 -1.I4388E+06 -3:334058+05 

0. -4.591 10E+05 -1.375418*05 -2-17838E*05 


ELECTRO I lENSITT (C0DL/I31 



I ,002738*08 
8.002718*08 
8.002448*08 
1.002338*08 
I. 002398*08 
1.002228*08 
1.082038*08 
1.081138*08 
1.001438*08 
8.001438*08 
6.001248*06 


«. 

0, 

0. 

0. 

0. 

I. 

0, 

0. 

0. 

0. 

0. 


4,779518-34 
4.65990E-34 
3.I40I3E-34 
1 .406578-34 
4. 30746E- 35 
1 .9C702E-35 
4.444788-36 
1.1 2494E-37 
1 .1 56I7E-37 
1 .29404E-38 
3.94796E-39 



nsnm ion density (cbul/mji 


-I.22O30E-1I -4.190B2E-10 6.3X69 6E*04 
-4.Q4349E-20 -2.33474E-M 1.I3094E+O* 
-I.98128E-23 -3.AO*33E-14 1.57I70E+I2 
-3.45754E-24 -9.87161E-17 9.49I23E-I2 
-4.00423E-29 -1 .3X1 *66-1? 1.5J377E-06 
-7.80949C-32 -1.3IS6IE-22 -3.175096-*? 
-I.I3749E-35 -1.36717E-25 -3.I9I86E-1Z 
-4.994 18E-38 -I.04I93E-2* -J.87234E-13 
-5.02938E-41 -4.881206-32 -2.0S243E-1Q 
•3. t HUE-44 -4.02980E-33 -I.I7I40E-2T 
-f . 17062E-45 -9.51431 £-37 -2.70027E-23 


4. 628516+43 1.54I57E+I4 2.5914^-02 

1.I8448E+03 1.243976*04 2.272506-0? 

3. 49O40E+04 5.44310E+03 U33334E-02 

I.44377E+04 1.S6312E+03 5.510006-03 

5.34784E+I2 2.64205E+I2 1.591046-01 
I .175296+41 2.77507E+01 3.261056-04 

1.11 4121-01 1.76397£+<?* 4.480C3E-0S 

4.414216-04 4.836356-02 4.72961E-04 

8.216116-07 1.6H32E-43 3.3642(6-07 
5.94279E-1 1 2.354506-05 1.618326-01 

“5. 191 146-14 I.7I179E-06 3. 138516-09 


1. 135166-12 1 .744176-25 -3.724376-37 

1.I3597E-T2 1.163436-25 -1.396*26-37 

6.79MSE-1J 1 . ( 18)86-25 “1. 4 3715c- 3 7 

3. 374461-13 5.74527E-26 -1.263856-37 

1 .24426C-13 2. 248426-26 -3.034516-38 

3.39442E-14 4.903196-27 -I.54444E-38 

7.7*0756-13 1.4267JE-27 -3.712036-39 

I.273I5E-15 7.978576-28 -7.0*0116-40 

1.5971 U-U 4.259176-2? -1.039246-40 
1 .529*8E-1 7 4.711 196-30 -I .2*5*26-41 

4.674476-1 i t .491556-30 -7.149396-42 


NET CHARGE 1EH5ITT (C0UL/I3) 

-2.371 34E-1I -5.91S766-07 6,3371 26+04 
-1.20133E-1? -2.744726-08 1.148I7C+04 
-1 .148846-22 -2.521 47E-1 I ' 1.5704*6+02 
-1.563846-25 -3.I0507E-14 1.482436-02 
-2.339746-28 -S.03440E-17 -5.297226-03 
-3.704746-31 -7.55985E-20 -I.61274E-07 
-4,554346-34 -9.556706-21 -2.917436-10 
-4.73947E-37 -I.01714E-25 -4.141336-13 
-4.25732E-40 -9. 293476-2? -4 .323416-16 
-3.37391E-43 -7.44I35E-32 -3.86271 E-M 
-9.180246-45 -2.040546-33 -1 .0994IE-20 


1.70731E+«5 -4.09123E+44 -1.91*646-01 
9. 979136+04 -3. 417486*04 -1.672236-01 
1.1 0?6iE+*4 -1 .585636+04 -9.69979E-02 
3.828946+03 -4.4*0506*03, -3,90t3iE-02 
1.299436+02 -2. 353946+42 -1.098816-02 
2.42659E+OS -7.463+4E+01 -2.15701E-O1 
-9.94973E-03 -4.641546+40 -2.961696-04 
-4.428J8E-0* -I.7I9SIE-01 -2.85009E-05 
-3.174506-14 -4.344206-03 -1.T26IS6-04 
-1.154216-01 -6.783456-05 181 326-01 

-5.152976-10 -8.2741 36-16 154426-01 


-3.923546-11 -4.452456-23 -4. 713236-34 
-3.54317E-I1 -1.827396-23 -4.64354E-34 
-2.296I4E-1 1 r2.541 25E-23 r3. 143276-34 
-T. 10303E-1 1 -1 .278336-21 -1.407S4E-34 
-3.979I8E-12 -4.9*2766-24 -6.312716-35 
-1.474I5E-12 -l.43843E 1 24 -1.909396-35 
-2.18J52E-U -3.24332E-25 -4.468546-36 
-3.342326-14 -3.64461 E“26 -8.112006-37 
-3.174*1 E-15 -7 .623476-27 -1.I5711E-37 
3. 41 459E-T 4 -I .*3?45E-28 -I.29524E-38 
9.349I8E-17 -2. 423I8E-28 -3.755UE-3? 


TEHPCRATU1E (KI 


3.010*06+02 
3.0000QE+02 
3.000006+02 
3.00000E+02 
3.00000E+02 
3.00000E+02 
J.COOOOE+O? 
3.00040E+02 
3. 00040E+02 
3. OOOOOE+02 
3.000006+02 


1.0*00*6+02 

3.00000E*02 

3.00000E+«2 

3.9COOOE+02 

3.0t004£*02 

3.0000IE+02 

3.04000E+02 

3.10000E+02 

3.00*00£+02 

3.40OC0€*02 

3.05000E+02 


3.00000E*02 

3.0I004E+I2 

3.0IOOIE+I2 

3.009OIE-M2 

3.00C04E+42 

3.00009E+42 

3.0«0«0E*02 

3.040I0E4O2 

3.A0010E+02 

3.0J040E+02 

3.90000E+I2 


3.0046IE+I2 
3. 00339 E+ 02 
3.00I9SE+02 
3.4O412E+0I 
3.00101 £+42 
3.40400E+92 
3.00004E+42 
3.00040E+02 
3.*OOIOE*02 
1.0000QE+02 
3.40IOOE+O2 


3.40I84E+92 
3.00469E+42 
3,00*31 E+02 
3.00I09E+02 
3.40101 E+*2 
3. 40040 £+02 
3.0IOOOE+I2 
3.0l000£*02 
4.04090€«02 
3.04000E+42 
3.0OOIOE+O2 


3.I0400E+02 
3.40900E+02 
3.40000E+02 
3. IOOflOG+02 
3, 10900E+02 
3. IOOCOE+02 
3-00900E+0? 
3.OI0OIE+O2 
T.0900IE+02 
3.0040CC+02 
3 .000006+02 


3. I0040E+02 
3.0*0*06+02 
3.0*0*06*02 
3.0*0106*02 
3»0*0*0E+02 
3.00000E+02 
3.*0O*OE+O2 
3.0*0*06*02 
3.0*0*06*02 
3. IC0*0E+02 
3.I0900E+02 


3,0*0*06*02 
3.0*00*6*02 
3.000016*02 
3.0*0*06*02 
3.0*0*0£*02 
J. 100*06*02 
3.000*06*02 
3.0*0*06*02 
3.0*0*06*02 
3.000006*02 
3.900106*02 


3.0*0*06*02 
3»0*00*E*02 
3.0*00*E+02 
3.0*00*6+02 
3.0*00*6+02 
3. *00*0£*02 
3.0*0*06*02 
3.010006*02 
3.010006*02 
3 .0*0*06*02 
3.«O0*OE*O2 


T8TAI NET CHARGE * 4.535992156-1 9 (C3UU 

TCTAt EiECTRON CH*R5E « 4 . I I222240E-1 7 (COIL) 

DIPOLE LEIGTH « 1 . 39206442E-OI W> 

CENTER OF ELECTROI CL08B « 7.11947742E-07 (N) 


CENTER OF DIPOIE « 2.*50*74I9E-07 (H) 

VARIABLES COIRESPDNDING T* AIOVE TINE iTEP DittPED TO TAPE 3 AT R6C0RI T 



HI - 


t ITR 


S ID1V « 


2 JPIV • 


I 


II« TIKE STE7 , WISER 

SEC SEC IIEIAT10NJ 

1.1S7SSC-U 1 .0513 TE-M 2 

2.0fSME-tJ 1.013901-14 2 

2.18974E-U 9.807B3E-I5 2 

2.284I4E-I1 9.49234E-I5 2 

2.374181-U 9.20I78E-15 2 

DISCIAIGE VARI4Bi.ES AT TIIE 2.17418E-13 SEC 


AIIAl ELECTIIC FIEU) (I0L1S/ HI 

1.034551*0* I. 104721*08 8.5I771E«08 7.771311+08 4.78I44E+08 8.3375CE+OI l.l*4*0e*0* ».03?03E*08 8. 014801+08 

1.O33T8E+08 8. 1OQ35E+08 8.05277E+0B 7.948I5E+0* 4.877411*08 1.471101+0* 8.103091+0* I. *34(21*01 *,I14441*0B 

8.032351*08 8.(87321*01 1.3)3321*01 *.000471*08 7.131341*18 8.3725SE*0i I.49339E+0* 1.034471*08 8.0I4I4£*08 

8.027871*06 6.075141*01 l.17SSTE«08 7.884031*08 7.417011*18 8.2419(1+01 8.171431*01 1.(31181*01 ». 115(01*08 

B.I2483E+I0 8. I404SE+0I *.118291*01 7.84757E*08 7.434501*08 8.13843E*Ot I.I4434E+OI S.*2843£*08 1.114(11*08 

8.023531*06 8,046871*08 8.043911+78 7.84488E*0» 7.74158E*0I 1.071081+01 1.047591*01 1.024741*08 8.013251*08 

8.020171*01 1.035211*06 8.03]4(E*08 7.873351*08 7.141421*08 8.031521+01 I.03447E+OI 1.021191*08 8.011781*08 

8.017021*08 8.02373E+OI 8.014311*01 7.118101*08 7.874071*08 I. 009871+01 S.I2414E+OI 1.017591*08 1.010471*08 

8.014111*08 8.018301+01 8.(03591*01 7.73821E*0» 7.925901*08 7.9784J1+0I 1.017131*01 I. 014301*08 8.0(9341*08 

8.011531*08 8.012441+08 7.997431+18 7.(52571+0* 7.74544E+08 7.193181+OS 8.11178E+0I 1.11(401*08 1. 108101*08 

8.009341*08 8.00848E+0S 7.79440E+0I 7.74312E+OI 7.938751*01 7.990851+01 S.S0730E+0I 1.008921*08 1.0(4731*08 


R4DIH. ELECTIK FIEU (VOLTS/ H) 


»- 7 . o. 0. 

#- -1.49I44E+04 -2.34020E+07 -1.374*41*08 
•- -3.011471*04 -2.394351+17 -1 .52481 E+(8 
»* -3.841101*04 -2.13583E+47 -7.S759SE+S7 
»- -4.22372E+04 -1.772211*07 -5.13052E+I7 
0- -4.25923E+04 -1.434951*07 -2.794171*77 
0. -(.040111+04 -1.120911+07 -1.404531+07 
0. -3.22(31 E+Oi -8.42194E+04 -7.719511*04 
0. -3.32575E+04 -4.59878E*04 -4.157341+04 
0. -2.91474E+04 -S.053I01«04 -4.040421*04 
0. -2.5218(1+04 -3. 831441+14 -2.77I45E+74 


«. 0 . 0 . 0. I. 

(.474721*07 1.3241(1+07 1.42473E+04 4.(57741*05 7. 

8.04724E+07 2.134781+07 3.IOOI1E+04 7.792431*05 7. 

7.48264E+I7 2.302871+07 3.954351+04 1.073101*04 0. 

4.793441+07 2.050241+07 4.444301+04 1.332711*04 0. 

2.5474OE+07 1.449071+07 4.537111+04 1.(91(91*04 (. 

1.38(941+07 1.241071+07 4. 338571+04 1.572511+04 0. 

1.034(91*04 (.439511+04 3.944S9E+04 1.517741*04 0. 

(.950571*04 7 .028071+04 3. 514271+04. 1.551551*04 (. 

3.203841+04 5.247471+04 3.051401+04 1.177*41*04 I. 

2.I53/4E+04 3.750711*04 2.410701*04 t.3S0I1£*04 0. 


1LEC1R0H I1ISITY ICSUL/H3! 

7. 575431-14 3.17091E-05 2.1321(1*02 
2.788351-14 1.310711-03 7.57371E+H 
I.275T2E-13 4.337451—07 Z.7544IE+O0 
1.71598E-17 1.0I4S41-08 5.713171-02 

1.447871-1? 7.350501-11 1.033941-03 

4. 498411-72 3.025951- 13 1.I1051E-OS 
2.18630E-24 2.020391- IS 2.522991-07 
5.495411-27 7. 22*131-18 2.40314E-79 
1.J2121E-29 2.34441E-20 2.009411-11 

2.925341-32 B. 4308(1-23 I. 193441-13 
2.138721-33' 5.063511-24 8.73I25E-I5 


3.422011+05 9.4571 9E+05 1.713301+03 
1.144441+05 8.45S3BE+05 1.424731+03 
1.523011+04 5.071191*05 8.311111*02 

4. 171331+03 1.48048E+7S 2.82I34E+02 

3.C7257E+02 2.781241+14 4. 391241+01 
1.34070!+0t 2.573181*03 1.045421+01 

3.7443*1-01 1.514431+72 1.201921+00 

4.18I18E-I3 5.75821 E+OO 7.524171-02 
7. 171421-15 1.513151-71 5.3S5741-03 

8.94224E-07 2.83392E-03 2.147891-04 
8.41 I20E-0* 1.82279E-04 4.410801-05 


1.143471-0* 1.815541-14 1.7(8(71-24 

I . 74302E-0I I. 474491-14 1.(44221-24 

4.34420E-05 1.084331-14 1.119441-24 

2.79S03E-OS 3.2(4301-15 3.544046-25 

4.728001-04 1.90141E-1S 2.1254(1-25 

2.011171-04 3. 211741-14 4.2I46IE-24 

3.544*2E-0? I. 079721-14 1.377431-24 
4.32374E-OS 1.710(31-17- 2.431551-27 

4,275(21-07 2.0II3JE-18 3.2732IE-2B 

3. I4573E-T • 1.749421-1? 1.514*51-29 

1.3151 -41—1 1 4.224451-20 1.578711-29 



POSITIVE 10* PENSJTT ( CCXft/MI > 


-2. 779 I 3E- 13 
-t.43A3?E-14 
-2.97701E-15 
-3.42079E-17 
-2.0A141E-T? 
-7.9I8J8E-22 
-2.22728E-24 
-3.23273E-27 
•I , I3227E-2? 
-2.37001E-32 
-1.7A230E-33 


-S.47089E-06 
-1.2937AE-04 
-3.731 77E-OI 
-6.3I301E-10 
-5.I214U-12 
-3.2I653E-J4 
-I.41519E-14 
-5.57T19E-I? 
-2.13230E-21 
-7.40252E-24 
-4.39600E-25 


6,43264E*04 
t.36315E*04 
1 .4061 1E+02 
1.2S27IE-IS 
2.6A3I2E-04 
2.1 3981E-04 
1.ZB918E-0I 
2.7702U-11 
-2.S3360C-13 
-3.69046E-I3 
-3.32H1E-16 


1.001 ltE+06 
i,44777E*05 
1.606I4E4O5 
1.90313E+O4 
1.0*333Et03 
2,0Bt4IE+#1 
I.2440SE-0I 
4.30541 £-03 
3.464I2E-05 
2.245T3E-07 
1.78785E-0I 


7.38M 3E+I5 
t.4#77«+*5 
3.3I583E+45 
1.01J24C+O3 
1.348J7E404 
1.5311 4E+I3 
1.01377E*01 
3.4I3IOE400 
8. 4454BE-02 
1.44H0E-O3 
1.11 2S5E-04 


4.41401E+OJ 
3.91 204E+02 
2.0I797E+02 
6. J8373E+01 
1.&242JE+0I 
2.46183E+OI 
3.07J04E-01 
2.4I427E-03 
1 .4244IE-03 
5.84634E-05 
I . I96I1E-03 


1.74952E-05 
1.S8727E-03 
P.71178C-04 
4.J8400E-0* 
1.37544E-OA 
3.229131-07 
5.S8504E-0I 
7.1 5776E-0f 
A.B43P7E-1 • 
4.93830E-1I 
1.J5398E-1I 


2 . 23417 E -15 -1 
2 .I 6241 E -15 -I 
I.I 3726 E -15 -i 
4 . 47578 E -14 
2 . 33272 E-I 6 - 
4 . 4 S 4 I 3 E -17 - 
U 34312 E -17 - 
2.1 36 ME -18 
2 . 61130 E -19 
2 . 4 I 233 E -20 
7 . 151 45 E -21 * 


NET CHARGE OENSITT (C0«L/13> 

— J.5193IE-1 3 -3.7J801E-05 6.41T33E+04 6.581O1E+05 
-I.24247E-13 -I.4403IC-O5 1.3 5 757 £+04 4.A6333E+05 
-4.24713E-T5 -4.71283 E-07 1.57453E+02 1. 33314 £+03 
-5.33677E-17 -1.14837E-08 6.I1473E-42 I.36400E+04 
-3.50T26E-1 9 ~9.*T3244E-n -7.6458ZE-04 7.71073E+02 
-1.44043E-21 -3.35460E-13 -t. 5 9 4531-05 1.S20S8E*01 
-4.41558E-24 -2.U191E-I5 -2.3M07E-47 4.99473E-02 
-I,0?2a3E-26 -7.78405E-1 « -Z.37344C-0P -2.48257E-03 
-2.45J48E-29 -2.75784E-Z0 -Z.03797E-T1 -5.72461 E-05 
-S.21545E-32 -J.37905E-23 -1.2J031E-13 -6.6743IE-07 
-3.701 03E-33 -5.523HE-24 -9.27040E-15 -6.83A35E-08 


-2.270O6E+I5 -1.43161E+01 -P.84713E-03 -I.5f2T5E‘14 - 
-2.247A2E+05 -1.233S2E+03 -1.042128-05 -1.44813E-14 - 
-1.70647E+05 -4.3I314E+02 -S.39S02E-03 -I.S24ME-15 - 
-4.A7434E+04 -2.13977E+02 -2.34943E-05 -4.3I677E-15 - 
“I.13t8 TE+04 -4.9670IE+01 -7.55256E-04 -I.646I3E-15 - 
-1.»35#4E*03 -7.9I542E+04 -1.75118E-04 -4.55656E-14 - 
-3.95249E+01 -8.M913E-01 -3.40I32E-O7 -1.45324E-17 - 
-2.35520E+00 -7.041 9IE-02 -3.J08I8E-0I -1.19674E-1? - 
-4.4S518E-02 -3.J3126E-05 -3.591 I2E-Q1 -1.120Z5E-18 “ 
1.39224E-I3 -I.58323E-04 -2.55210E-tO -I.72USE-I9 - 
2.40T95E-04 -3.213J8E-03 -i.fo2loE-U -5.43948E-20 - 


TENPERATIRE <K) 


3.0I0OCE402 
3 . OOQOOE +02 
3.00000E+02 
3.00000E+02 
3.O0OOOE+O2 
3.00000E+02 
3.00000E+02 
3.00000E+02 
3.00000E+02 
3.00000E+02 
3.00000E+02 


3.00004E+02 
3.00001E+02 
3.00400C+02 
3.00000E402 
3.00000E*02 
3.00000E+02 
3.00000E+02 
3.00000E+02 
3.00001E+02 
3. 10000E+02 
3.00000E+02 


3.00002E+42 

3.50001E+92 

3.00CCOE+02 

3.01001E+02 

3.O0OIOE4O2 

3.00001E402 

3.IOOIOE+12 

3.00010E+02 

3.00000E+02 

3.00004E+02 

3.00000E+02 


3.02725E+I2 
3.01 143 E+02 
3.10448E+02 
3.00054E+02 
3.0000JE*02 
3.00000E+02 
3.IOO0OE+O2 
3.OCOOSE+02 
3.00400E+02 
J.OCOOOE+02 
J.00090E+02 


3.43391E+I2 
1.02I53E4O2 
3.0I549E4O2 
3.00473E-+I2 
3.00I77E+42 
3.00I07E+02 
3.0IO0OE+O2 
3.0IOOOE+I2 
3.00600E+02 
3.00I00E+42 
I . OOOIOE+02 


3-I0002E+02 
3. 00002E+02 
3.04001E+02 
3. IOOOOt+02 
3.00001E+02 
3.IOO0OE+O2 
3.000CCE+02 
3.00000C+02 
3.IOOOCE+02 
3.00000E+0? 
3.00000E+02 


3. 40000E+02 
3.O0OIOE4O2 
3.00000E+02 
3.00000E+02 
3.00000E+02 
3.00000E402 
3.IO0OOE+O2 
3. OOOOOE+02 
3.00000E+02 
3.00000E+02 
3.00000E+02 


J.*00*0£402 

3.00090E402 

3.00040£t02 

3.0I010E+02 

3.00040E«02 

3.40040E402 

3.00040E+02 

3.40000E+02 

3.00000E402 

3 .ocoaoE*o 2 

3.04000E402 


TOTAL *ET CHARSE ■ 1 . I2205654E-18 (C3U.) 

TOTAL CL£CT20f4 CHARGE - 2.8U490I8E-H <COW_) 

DIPOLE LEIGH s t .445 10A77E-08 (H) 

CENTER OF ELECTROI aOOD « 2.37335770E-07 {«» 

CENTER Of DIPOLE • 2.300J0236E-I7 (HI 

VARIABLES CORRESPONDING TC AIO?E TINE STEP 50NPED TO TAPE 3 AT RECCRI 4 


NI * 0 NTR * 1 IBIV * 1 JDTV * 



. 1 (54 4E* 25 
.04641E-25 
.34573E-26 
. 460 P 5 E -24 
.31732E-26 
1. 079f 7E- 27 
.I58I6E-28 
I.5I967E-28 
2.M654E-29 
2.21183E-30 
1.45303E-3O 


1.01 50 IE-24 
1.77383E-24 
1 . IB270E-24 
5.9101 4E-25 
2.24340E-25 
4.42280E-26 
1 . 48922£’26 
2.590S2E-2 7 
3. 54713c- 28 
3 . 745 53c~ 2? 

I.7H21E-29 


3.40000E»02 
3.04000E*02 
3.0»00#E*02 
3.0I040E402 
3.04000E+02 
3 .10000E*02 
3.S0050E402 
3.0l000£*02 
3.04000E402 
1.O0O4OE+O2 
3.0I040E402 



Appendix 7 


ANALYTIC APPROACH TO THE CALCULATION 
OF IONIZATION COEFFICIENTS 


The interest of calculating ionization coefficients 
in semiconductors dates back to about three decades ago. 

McKay (1954) derived from experimental data a curve of o^ 
versus applied electric field. Wolff (1954) has developed a 
method to compute a^‘ by expanding the electron distribution 
function in terms of Legendre polynomials and keeping the first 
two terms in solving the Boltzmann transport equation in steady- 
state. This approximation is justified at high electric fields. 
Shockley (1961) has conjectured that the steady-state electron 
distribution function must have a spike in the direction of 
the electric field and, by neglecting all other electrons not 
in this spike, he was able to obtain an a whose logarithm is 
inversely proportional to the electric field. This approach 
is reasonable for low electric fields. Baraff (1962) solved 
numerically the Boltzmann equation and obtained an important 
plot of ctX versus E^/e£\ for various ntaf values. Here, X 
is the optical mean free path', E^ , the energy gap between 
the valance band and the conduction band, rito the energy of the 
optical phonon, and <? the external electric field. Further 
publication by Baraff (1964) concentrated on the high field 
case in which spherical harmonics expansion was made and a 
truncation procedure was used to obtain a closed system of 
equations. Temperature consideration was also given later on 
(Okuto and Crowell, 1972; Crowell and Sze , 1966). Finally, 
macroscopic ionization rates from theory were studied and found 
to be generally different due to the possible field dependence 
of the carrier's drift velocities in the avalanche regime 
(Beni and Capasso, 1979). 



In this appendix, we will approach the problem of 
compuxing ionization coefficients - from a point of view different 
from that of solving Boltzmann equation numerically (Baraff, 

1962, 1964). Specifically, we treat the collision processes 
as Markov processes. If X is the mean free path of a collision, 
then the probability of no collision for an electron ' traveling a 
distance x is exp(-x/X); the probability of one collision in 
dx is dx/X. The mean distance of a collision is therefore: 

CO 

x « j xe~ x/1 ~ = X (A7.1) 

o 

In the following, we will concentrate on the applica- 
tion of the concept of Markov processes to the calculation of 
the ionization coefficient in the case that scattering effect 
is neglected. 

An electron is released at E=0, gaining energy from 
the constant electric field S, losing energy n<u by emitting 
an optical phonon with a constant mean free path X. What is 
the mean distance in which the electron will reach the ioniza- 
tion energy E^? This is the question we would like to address 
analytically . 

Shockley's (1961) conjecture about the aforementioned 
spike is based on the Markov concept in a very trivial way, as 
follows. The number of electrons, per unit volume, which have 
survived the transport to energy E^ without collisions is 
proportional to expC-Ek/e^X) where E i /ei? is the distance traveled 
by the electron to reach energy E^. Those are the electrons 
which can ionize electrons in the valance band. The number of 
ionizations in a unit length is therefore proportional to 
exp (-E^/ edTX) : 
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a 


(A7.2) 


= c expC-E^/e^X) 


Results from our analysis, on the other hand, 
possess 'more complicated dependence on (E^/ed^X). 

We choose E^ = (n-rl)Tuo for simplification of the 
analysis. Let XZ^ he the mean distance in which an electron 
released with E = at origin will reach Ej_ . Define 
R = X = E^/e£\, Y - aX , a = exp(-RX), 8 = RX, and 

AZ s -fua/eSx = 8. Apparently, Z^ = Z q - AZ. We define K(i) 
by: 


Z ± = Z - K(i) AZ 


CA7.3) 


Obviously, K(l) = 1. The recursive formula for K(i-tl) can be 
obtained by considering the probabilities of possible final 
states the electron with energy i tzw at origin can end up with. 

As shown in Figure (A7-1) , it has probability a to go up to 
(i+l>?fto and then reaches E^ at (Z i+1 + AZ). It has ctP(i,j) 
probability to go to energy j7!co( j^l) , and then reaches E. at 

i_ -i-i 

(Z. + AZ) where P('i,j) =8 /(i-j+l)I Finally, the remaining 

probability is to drop down to below Tzco and reaches E^ at 
(Z x + AZ). 

i 


Z. 

l 


' a < Z i +1 


+ AZ) + a 


j=2 


P(i, 3 ) 


(Z . 
3 


+ AZ) 


+ ( 1 — a 


E 

n — O 


P(i , j ) ) (Z 1 + AZ) ‘ . (A7.4) 
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By substituting Z, = Z -K(k)AZ into the above equation, one 

K O 

gets : 

i 


aK(i+l) 


1 + K(i) - 


a 


'y 1 P(i, j)K(j ) 
3=2 


• - (1 - a 


a P(i. 0)K(1) 
3=2 


(A7.5) 


This equation and K(l)" = 1 thus determine all K(k)'s. It is 
noted that eq. (A7.5) depends only on the two dimensionless 
quantities R = and X = E^/ef^A. 

At energies greater than E^ = (n+l) 7 zw, the electron 
is assumed to be capable of ionizing electrons from the valance 
band with a mean free path X^. We also assume infinite mean- 
free path for optical phonon emission when the electron energy 
is greater than E^. Under these circumstances, = A^/A. 

Eqs . (A7. 3) and (A7.5) then determine Z q as a function of R 
and X . . 


by 


The ionization coefficient per unit length is given 


1 

a ” Z~T 

. o 

aA = 1/Z = f (R’,X) 

Y = f(R,X) (A7.6) 
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We have carried out the calculation for a number of 
R's. Fig. (A7-2) shows the results with X^ = X. The curve for 
R=0 (n=<» in this case) is obtained differently as follows. In 
the limit of zero energy for the optical phonon, the electron 
does not lose energy. Z Q X = + X^, i.e. , 

Y = iix for x i = x - 

The curves in Fig. (A7-2) do not have the simple 
behavior (a straight line in semi-log plot) of eq. (A7.2), 
except for R n* 0.5. This is anticipated because there are 
terms of the type e mRX in eq. (A7.5) for m = 1, . . . ,n+l. It 
is only when m=n+l does e = e , as eq. (A7.2). 

Compared with Baraff's (1962) work in which the 
isotropic scattering effect is included, our results agree 
wixh his only for the small X where the electric field is so 
high that the direction of the electron after an. isotropic 
scattering almost immediately goes back to that of the applied 
field. For large X (small electric field), the scattering 
effect should be very important and our results do deviate 
greatly from Baraff's (1962) work. 

o 

For the case of Teflon, we have put X = 26 A, 
independent of the electron energy Ceq. (A. 2. 8) does depend 
on energy, however!, for the optical phonon mean free path. 

= 6.5 eV, fzco = 0.11 eV. The calculated ionization length 

by this model is plotted in Fig. (A7.3) as a function of the 

applied electric field. The region for which the electric 

9 

field is above 10 V /m is characterized as high field region 

and the present model is more reliable. Comparison with 

Fig. 2.6 which is obtained by the Monte Carlo approach with 

energy-dependent X and with scattering, shows reasonable 

agreement between the two approaches; the ionization lengths 

9 

being within a factor of two of each other for E > 10 V/m. 
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FIGURE A7.1. ALL POSSIBLE STATES OF AM ELECTRON 
RELEASED WITH E=iR(ii AT ORIGIN 
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X = E_j /esA 


FIGURE A7.2. IONIZATION COEFFICIENT AS A FUNCTION 
OF ELECTRIC FIELD FOR VARIOUS R 


A7-7 



IONIZATION LENGTH (m) 


10“ f 






ELECTRIC FIELD CVOLTS/m) 



FIGURE A7.3. IONIZATION LENGTH AS A FUNCTION OF 
EXTERNAL ELECTRIC FIELD FOR CF ? 
(SCATTERING IS NEGLECTED) 



Appendix 8 

ELECTRIC FIELD PENETRATION INTO AN ISOTROPIC PLASMA 


In one dimension the equations that describe 
the steady state penetration of an electric field into 
an isotropic plasma with fixed ions are: 

O = u n(x)E(x) - D (A. 8.1) 

Ji (n( „) - n(x)> ■ (A. 8. 2) 

dx £ o 

where n(x) is the electron number density, E(x) is the 
electric field, y is the charge mobility for electrons ; 

D is the diffusion coefficient for electrons, e is the 
charge on an electron (e<0), and is the permittivity 
of free space. 

Equation (A. 8.1) is a statement of the balance 
of particle current due to diffusion and due to E field 
drifting. Equation (A. 8.2) is Poisson's equation. Here, 
we will assume that u and D are constant, although in 
general this is not a restriction to the numerical solu- 
tions of Eqs . (A. 8.1) and (A. 8. 2). 

We are interested in the problem of a plasma 
with a density n and electron temperature kT, existing 

CO 

from x=o to x=«=, where at x=0 the electric field has a 
known value. Equations (A. 8.1) and (A. 8.2) are easily 
nondimensionalized by the following substitutions: 


A8-1 • 



n(x) = nji(x) 


E(x> = E q E(x) 

x = x x 
o 


(A. 8. 3) 


where : 


n = n(<») 

E q = E (0 ) 

E kT 
o 

' 2 
n e 

CO 

In the definition of x we have used the Einstein relation 

o 

that states: 



U = _e 
D kT 

x q is the standard definition for d Debye length. 

In nondimens ional "hat M variables Eqs . (A. 8.1) 
and (A, 8. 2) become: 



3E 

A 

3x 


"I * A 
^(1-n) 


(A. 8. 4) 


(A. 8 .5) 
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where : 






n kT 

CO 


1/2 


The parameter 5 is the square root of the ratio of the 
electric field energy density at x=0 to the electron 
thermal energy density at x=«. In what follows we drop 
the "hat" notation and assume all variables to be non- 
dimensional . 


Refering to Eq. (A. 8. 3 ) the boundary condi- 
tions for the problem we wish to solve are: 


n(co) .= 1 

E ( 0 ) = ±1 


(A. 8. 6) 


Before we solve Eqs. (A. 8. 4 ) and (A. 8 - 5 ) with 
boundary conditions Eq. (A.S.6) it is profitable to look 
at the first integral. If we multiply Eq. (A. 8. 5 ) by E 
and replace 11E and E by derivatives of n using Eq. (A. 8. 4 ) 
it is simple to show: 


_d 

dx 



£n(n) 


n 


0 


Integrating this expression we immediately get: 
r 2 2 

E + £n(n) - n = const. 


Applying the boundary conditions at x=“ (we tacitly 
assume for the moment E(“)= 0 , that is the plasma shields 
the E field) we get: 
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-g— E 2 + £n(n) + 1 - n = 0 

which is valid for all x. In particular for x=0 we have: 

2 

r + £n(n ) + 1 - n =0 (A. 8-7) 

2 O O 

where : 

n Q = n(0) 

Figure (A. 8.1) is a plot of left hand side of Eq . (A. 8. 7) 

2 

versus n . It is obvious from the plot that for c, >0 there 
are always two solutions which make the left hand side of 
Eq. (A.S.7) zero. One solution is always less than 1 cor- 
responding to the case where the E field is repeling elec- 
trons from the surface. The other solution is always 
greater than 1 corresponding to the situation where the E 
field attracts electrons to the surface. 


If £>>1, from Eq. (A. 8. 7) we can estimate the 


two solutions for n to be: 

o 


n = 
o 


e - \ + 2 

i + s! 

2 


For £<<1 the corresponding solutions would be: 

1 - C 
n 

o 

1 + C 
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Figures (A. 8.2a) and (A. 8.2b) are plots of n Q versus 5 
for both branches. The dotted lines correspond to the 
estimates for large and small 5 . 

Knowing n Q makes Eqs . (A. 8. 4) and (A. 8. 5)' 
straight forward to integrate, since by the first inte- 
gral we know two boundary conditions at x=0 . If we 
divide Eq. (A. 8 . 4 ) by n and differentiate we obtain: 


— ^ 3 - iln(n) 
3x 2 



(A. 8 . 8 )' 


Substituting Eq. (A. 8. 5 ) into Eq. (A. 8. 8) and letting 
a(x) = £n(n) Eq. (A.S.8) becomes: 


■2-f = (e a - 1) 


(A. 8 .9) 


with boundary conditions 


a(0) = £n(n Q ) 


d a 
dx 



= ±C 


(A .8 .10) 


Eq. (A. 8. 9) with boundary condition Eq. (A. 8. 10) is in 
standard form for numerical integration. Once' the solu- 
tion to a(x) is obtained the profiles n and E are simply 
computed by: 
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n (x) 


= e a(x) 


E(x) 


1 dg(x) 
5 dx 


Figure (A. 8. 4) shows a numerical solution to Eq. (A. 8.. 9) 
for 5=2, for both attracting and repelling -E fields. The 
profile E(x) r.eally does fail away to zero at x=«>, veri- 
fying our --as sumption about the first integral. 


As a measure of the penetration of' the E field 
we define an E field penetration deprh as the first mo- 
ment of the E field, profile-: 

r°xE(x)dx 
x = 1 ° 

■ E /“E(x)dx 


Likewise we define the distance over which the density 
varies appreciably as: 


j“x]n(x) r l.| dx 
’ J"|n(x) - 1 1 dx 

Figures (A. 8. 4) and (A. 8.. 5) are a plot of A E and A^ versus 
5 for the cases where the E field attracts and repels the 
electrons from the surface, respectively. A-g and A^ were 
computed from the numerical solution of Eq. (-A.8.10). 

The x integration was terminated when both A £ and A N were 
changing by less than a hundredth of a percent. 


In Figure (A. 8.4), the ca.se where the E field 
attracts electrons to the surface, the penetration depth 
decreases with increasing E (recall 5 ~E for constant 
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thermal energy). This reflects the fact -that as E in- 
creases more and mor.e electrons migrate to the surface, 
more effectively shielding the E field. For small E 
(no matter how small) the penetration depth is fixed at 
one Debye length. The density damps more quickly than 
the E field. 

In Figure (A. 8. 5), the case where electrons 
are repelled from the surface, the penetration depth 
increases with increasing E. This happens because as 
E increases the density can only decrease to zero ef- 
fectively letting the E field penetrate undamped in 
this very low density region. For this case the E 
field is damped more effectively than the density. 

Finally, we present an actual physical case 
of interest in this paper. For rhe following parameters 

E - 2 x 10 8 volts/m 
o ' 

kT = -1 ev. 

n Q = 1 x 10 18 /cm 3 

Figure (A. 8. 6) is a plot of the E field and density vs. 
distance for both the attracting and repelling cases. 

The penetration lengths are computed as: 

= 52.9 a. ') 

\ attracting 
X E = 59.1 8 ) 
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repelling 


X. T = 97.7 a 
N- 

X E ■= 82.9 a 
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FIGURE A. 8.1. FIRST INTEGRAL OF EQS. (A. 3.4) AND (A. 0.5) Vs. n Q . THE 
VALUES OF n Q FOR WHICH THE FIRST INTEGRAL IS ZERO GIVE 
THE BOUNDARY CONDITION AT x=o NEEDED TO SOLVE EQS. (A. 8. 4) 
AND (A. 8. 5) NUMERICALLY . 


FIGURE A.8.2A. n Q vs. 5 FOR THE CASE WHERE THE E FIELD IS 

attracting ELECTRONS TO THE SURFACE. 

DOTTED LIMES ARE ANALYTIC APPROXIMATIONS. 



C 


FIGURE A. 8.2S. n Q vs. 5 FOR THE CASE WHERE THE E FIELD 

REPELS ELECTRONS FROM THE SURFACE. 

DOTTED LINES ARE ANALYTIC APPROXIMATIONS. 
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FIGURE A. 8. 5. 


E FIELD PENETRATION LENGTH’ (X E ) AMD DENSITY 
SCALE LENGTH (A fJ ) FOR REPELLING SOLUTION 
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ANGSTROMS 


FIGURE A. 8. 6 . THE ELECTRIC FIELD AND DENSITY vs. 

DISTANCE FOR A TYPICAL PHYSICAL 
CASE OF INTEREST IN THIS PAPER 
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Appendix 9 

chapman-enskog’ theory 


The Boltzmann equation can be written in the 
following form: 


3f 

3t 


+ v 


V x f 


eE 

m 


7-if 

v 



K(v,v o )f 


(A. 9.1) 


Let e be the expansion parameter according to the standard 
Chapman-Enskog procedure (Fersiger and Kaper, 1972), and 
expand f 


= f + ef i + 
o 1 


(A. 9. 2) 


In order to include the effects of the electric field, as 
well as the scattering process in the lowest order solution 
f , we formally order these two terms of 0(l/e), The equa- 
tion of order 1/e is, therefore: 


eE 

m 


V-if 
v o 


■ / 


dv K(v, v )f 
o v ’ o o 


(A. 9. 3) 


The solution of this equation determines f Q . The equation 
of order one is the following equation: 



i 

v 


V-^f 

X o 


eE 

m 




/' 


dv o 


K(v,v o )f 1 


(A. 9. 4) 
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According to the Chapman-Enskog procedure, the time depend- 
ence of the distribution function is determined through the 
dependence of f on the density, drift velocity, and tempera- 
ture, and this is determined by the solution to the lowest 
order equation. If we -restrict our attention to the case 
of constant drift velocity and temperature, we can write: 


8f o 

3t 


3f 

o 

3n 


3n 

3t 


(A. 9. 5a) 


3f 
o 

3x 


3n 

—A 

3x 


3f 

o 

3n 


(A. 9. 5b) 





(A. 9. 5c) 


To lowest order in e: 


3n 

3t 







?^>n 

x 


0 


(A. 9. 6) 


where 



1 f 

5 J fvdv 


.(A. 9. 7) 


Therefore, eq. (A. 9.4) can be written-: 


* V-nl = J dv o K(v,v o )f 1 + ^ • V i ± 1 . (A. 9. 8) 
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We can obtain a particular solution of this equation by taking 
the first velocity moment: 


C<v . v .’> 
i 3 


f -a 

<v.xv.>1^ = / dv dv v.Kf- 
1 3 3 x i Jo 3 1 


(A. 9. 9) 


where 


<v.v .> 
i 3 


<v. > 
i 




v .dv 
J 


fg V. 

J to O 1 


dv 


(A. 9. 10) 


(A. 9. 11) 


Define: 


y -A -A _A _i _i 

dv dv v .K(v,v )f (v ) 
o 3 o 1 o 

= f dv v. K-. (v^) f- (v ) 
J o jo 1 o 1 o 


(A. 9. 12) 


A solution of eq. (A. 9. 9) is the following: 


f 


1 



Cv - 


-A 



1 


V-»n 

x 


(A. 9. 13) 
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A knowledge of f^ allows us to compute the average electron 
velocity to order 1: 

<v> = v^ + ^ y*dvf^(v)vdv (A. 9. 14) 

and we can write: 

<v> - V D “ D • Vn (A. 9. 15) 


where D is the diffusion coefficient and can be written: 


D 



CvCv-v D )l 


(A. 9. 16) 


Because of the symmetry of the scattering process, D is 
diagonal, but with different values parallel and perpendicular 
to the E field. 
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